Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9780AC22.5% 540.3 / 197.4 257intergenic (‑188/+120)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (99/87);  new base C (54/0);  total (170/87)
Fisher's exact test for biased strand distribution p-value = 5.20e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-02

Only 100 of 260 total aligned reads displayed.

GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076
                                                                                                   |                                                                                                       
 aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                        >  2:1748445/1‑100 (MQ=255)
   gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                      >  2:1400367/1‑100 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    <  2:3697246/100‑1 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    <  2:1868168/100‑1 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    <  2:1203372/100‑1 (MQ=255)
     aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                    <  2:1583895/100‑1 (MQ=255)
      ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                   >  2:723952/1‑100 (MQ=255)
        gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                  >  1:2123468/1‑99 (MQ=255)
           gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                              >  1:1841179/1‑100 (MQ=255)
             cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                            >  1:1621528/1‑100 (MQ=255)
              ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                           >  2:2463373/1‑100 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:958081/100‑1 (MQ=255)
               gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                          <  2:818949/100‑1 (MQ=255)
                 tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                        >  1:2555284/1‑100 (MQ=255)
                  ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                       >  1:4095968/1‑100 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGCCTAC‑a                                                                                     >  1:343967/1‑100 (MQ=255)
                   gcgcTGCACGCAGAGGCACAGCGCCACAGGGTATACAGATTGTTGACAAGGGGCGCCTTGCTTATGCCGGAGGCGGCGTCAACGCCTTATCCGCCCTCC‑a                                                                                     >  2:3850843/1‑100 (MQ=37)
                    cgcTGCACGCAGAGGCACAGCGCCACCGGGTATATAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAGCCGCCCTAC‑aa                                                                                    <  2:2198763/100‑1 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                >  1:3633464/1‑100 (MQ=255)
                        gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                >  2:2566778/1‑100 (MQ=255)
                         cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                               <  2:3602570/100‑1 (MQ=255)
                           cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                             <  2:795566/100‑1 (MQ=255)
                             cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                           >  1:3344190/1‑100 (MQ=255)
                               gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                         <  1:57648/100‑1 (MQ=255)
                               gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                         >  2:2842144/1‑100 (MQ=255)
                                agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                        <  1:74454/100‑1 (MQ=255)
                                 ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                       >  1:2482844/1‑100 (MQ=255)
                                  ncacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                      >  1:970644/2‑100 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     <  2:3380946/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     <  2:181838/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     <  1:1625364/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     <  1:3485065/100‑1 (MQ=255)
                                   cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATACGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                     >  2:2677384/1‑100 (MQ=255)
                                     caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                   <  2:1037196/100‑1 (MQ=255)
                                        cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                <  1:44975/100‑1 (MQ=255)
                                         gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                               >  1:3750925/1‑100 (MQ=255)
                                          ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                              >  2:3973068/1‑100 (MQ=255)
                                            agcGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                            >  1:2836025/3‑100 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                         <  1:4132410/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                         <  2:3898464/100‑1 (MQ=255)
                                              cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                         <  1:1561023/100‑1 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                      >  1:3845921/1‑100 (MQ=255)
                                                ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                      <  1:2891323/100‑1 (MQ=255)
                                                   ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑Attt                                                   <  2:3959336/100‑1 (MQ=255)
                                                     aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTc                                                 <  1:3302885/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑g                                              <  2:698052/100‑1 (MQ=255)
                                                      cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑g                                              <  1:1996007/100‑1 (MQ=255)
                                                        gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTa                                            <  2:3468785/100‑1 (MQ=255)
                                                         aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAg                                           >  2:3096364/1‑100 (MQ=255)
                                                          ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAgg                                          >  1:276384/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGa                                     >  2:3418752/1‑100 (MQ=255)
                                                               gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGa                                     >  1:1495035/1‑100 (MQ=255)
                                                                 cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATa                                   <  1:3056097/100‑1 (MQ=255)
                                                                  aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATaa                                  <  2:2514453/100‑1 (MQ=255)
                                                                      tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCg                               <  2:1099332/99‑1 (MQ=255)
                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                       >  1:1268481/8‑100 (MQ=18)
                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATCTCCGTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:538629/8‑100 (MQ=18)
                                                                             tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑ACTAT‑‑‑GTAGGCCTGATAAGCATAGCGCAt                       >  1:2056668/8‑100 (MQ=18)
                                                                              gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                      <  2:3614070/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                     <  1:635868/100‑1 (MQ=255)
                                                                               tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                     >  1:3918225/1‑100 (MQ=255)
                                                                               tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GACTAT‑‑GTAGGCCTGATAAGCATAGCGCATCa                     >  1:3775130/6‑100 (MQ=18)
                                                                                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                    <  1:1400378/100‑1 (MQ=255)
                                                                                ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                    <  2:1841165/100‑1 (MQ=255)
                                                                                gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg                    >  1:3268704/5‑100 (MQ=18)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                   >  2:1429436/1‑100 (MQ=255)
                                                                                 tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                   >  2:420042/1‑100 (MQ=255)
                                                                                  aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  1:3984123/3‑100 (MQ=18)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcg                 >  1:3600999/3‑99 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  1:3571019/3‑100 (MQ=21)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGATCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  2:2802287/3‑100 (MQ=17)
                                                                                  aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AAACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc                  >  2:2604278/3‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GCCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:2678399/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:3112094/2‑100 (MQ=17)
                                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCG‑AG‑CTTCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCa                 >  1:2938883/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  2:1461546/2‑100 (MQ=18)
                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                 >  1:2089800/2‑100 (MQ=21)
                                                                                    gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                >  2:3145687/1‑100 (MQ=255)
                                                                                    gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                >  1:1207656/1‑99 (MQ=17)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat               >  1:973996/9‑98 (MQ=25)
                                                                                     ccggataaggcgTTCACGCCGCATCCGGCCTA‑‑TAAATCATGCAAATCCAAAATATTGCA‑GG‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt              >  2:2003735/9‑100 (MQ=21)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt               >  1:4021250/1‑100 (MQ=255)
                                                                                     ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt               >  1:2313246/1‑100 (MQ=255)
                                                                                     ccGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                 >  1:2346837/1‑98 (MQ=17)
                                                                                      cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GCCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatt              >  1:3198393/8‑100 (MQ=18)
                                                                                      cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatt              >  1:241319/8‑100 (MQ=21)
                                                                                      cGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt              <  1:2146669/100‑1 (MQ=255)
                                                                                       ggataaggcgTTTACGCCGCATCTGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt             >  2:2318665/7‑100 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CTCCAC‑‑GTAGGCCTGATAAGCGTAACGCATCAGGCa                 >  1:1913010/7‑96 (MQ=18)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaattc             >  1:834315/7‑96 (MQ=21)
                                                                                       ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGATTAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt             >  1:1007299/7‑100 (MQ=17)
                                                                                       ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑ACGAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta             >  1:2949985/7‑96 (MQ=25)
                                                                                        gATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt            >  1:3813783/1‑100 (MQ=255)
                                                                                         aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg           <  2:584437/100‑1 (MQ=255)
                                                                                           gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc         <  1:2574081/100‑1 (MQ=255)
                                                                                              gcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt      >  1:1061649/1‑100 (MQ=255)
                                                                                               cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg     <  2:45364/100‑1 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc    >  1:2240378/1‑100 (MQ=255)
                                                                                                gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc    <  2:3123626/100‑1 (MQ=255)
                                                                                                 ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc   >  1:605600/1‑100 (MQ=255)
                                                                                                   |                                                                                                       
GAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCA  >  NC_000913/3705879‑3706076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑