Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 15.1% intergenic (+24/+79) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3110CG15.1% 424.5 / 60.8 152intergenic (+24/+79)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (68/61);  new base G (15/8);  total (83/69)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 154 total aligned reads displayed.

ACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACC  >  NC_000913/2286215‑2286410
                                                                                                |                                                                                                   
aCGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGACGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTa                                                                                                  <  2:763193/100‑1 (MQ=255)
 cGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                 >  1:682130/1‑100 (MQ=255)
   tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgg                                                                                               <  1:487051/100‑1 (MQ=255)
      cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                            <  1:3059582/100‑1 (MQ=255)
       gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGGTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTATggttgg                                                                                           >  2:2636116/1‑100 (MQ=255)
        cAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggt                                                                                          <  1:2109889/100‑1 (MQ=255)
          gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttt                                                                                        >  1:1439399/1‑100 (MQ=255)
           ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                       <  1:16518/100‑1 (MQ=39)
            cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                      >  1:1772439/1‑100 (MQ=255)
             tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttttt                                                                                     <  2:2166632/100‑1 (MQ=255)
             tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                     >  1:2023417/1‑100 (MQ=255)
              ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    <  1:1230070/100‑1 (MQ=255)
              ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    >  1:3381950/1‑100 (MQ=255)
              ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    <  1:250881/100‑1 (MQ=255)
              ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGGTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    >  2:2095871/1‑100 (MQ=255)
              ggTAGCGCACCGTCATGGGATGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                    >  2:1840527/1‑100 (MQ=255)
               gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGATCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                   <  2:2675236/100‑3 (MQ=38)
               gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                   <  1:928837/100‑3 (MQ=38)
               ccgatctaCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                   <  1:240882/93‑1 (MQ=255)
                  gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAACAAACCAACCGTAAGGGTTGGttttttcttag                                                                                <  1:2655190/100‑6 (MQ=37)
                  gcgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTatat                                                                                >  2:2237394/1‑100 (MQ=255)
                  gcgaACCGTCATGGGGTGTCGGGGGTCAGAGGTTCAAATCCGATCGTGCTGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                 <  2:753240/99‑5 (MQ=25)
                   cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                               >  1:3545287/1‑100 (MQ=255)
                    gcACCGTCATGGGGTGTCGGGGGTTGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                              <  1:3124538/100‑1 (MQ=255)
                    gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGCAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                              >  2:3129933/1‑100 (MQ=255)
                      aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtaagatcgg                                                                            >  2:240882/1‑93 (MQ=255)
                       ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                           <  1:1779928/100‑1 (MQ=255)
                          tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                        <  1:1305078/100‑1 (MQ=255)
                             tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                     >  2:844961/1‑100 (MQ=255)
                              ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                    >  1:1367935/1‑100 (MQ=255)
                               nggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                   >  1:2628895/2‑100 (MQ=255)
                                ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                   >  1:434317/1‑99 (MQ=255)
                                ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                  <  1:955823/100‑1 (MQ=255)
                                 gtgtAGGCGGTCGCAGATTCAAATCCTCCCGTTCCGAACAAAAATCCCAAGCAAAACCCAACCCTTCCCGTTTGTTTTCTTATATCTGCAATTAATTc                                                                   <  2:2330998/98‑1 (MQ=255)
                                   gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGaag                                                               >  2:3447234/1‑98 (MQ=255)
                                   gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                               >  1:514190/1‑100 (MQ=255)
                                     cGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                             >  1:2030814/1‑100 (MQ=255)
                                       ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                           <  1:2210826/100‑1 (MQ=255)
                                       ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCGTACTGTTGGTGTTTTTATATCTGCAATTAATTCGATAAACa                                                           >  2:325583/1‑100 (MQ=255)
                                       gaggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                           <  1:84742/98‑1 (MQ=255)
                                        gggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAg                                                          >  1:1318808/1‑100 (MQ=255)
                                          gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAc                                                        >  2:3910096/1‑100 (MQ=255)
                                            cGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                                      >  1:1633601/1‑100 (MQ=255)
                                             ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                     >  2:3821799/1‑100 (MQ=255)
                                              gAGGTTCAAATCCTCTCGTGCCGACCAAACATCTCAGAAAAAACCACCCCTTACGCTTAGTCTTCTTATAGCTGCCATTAATACGATACACCGACCCTg                                                    <  2:1163732/99‑1 (MQ=255)
                                               aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                   <  2:1850410/100‑1 (MQ=255)
                                                ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                  <  1:48993/100‑1 (MQ=255)
                                                 gTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCCACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                  <  1:3280334/99‑1 (MQ=255)
                                                   tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                               <  1:359483/100‑1 (MQ=255)
                                                   tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                               <  2:167255/100‑1 (MQ=255)
                                                    cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                              >  1:1745113/1‑100 (MQ=255)
                                                    cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                              >  2:3057092/1‑100 (MQ=255)
                                                     aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                             <  1:870750/100‑1 (MQ=255)
                                                         ccTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                         <  2:200762/100‑1 (MQ=255)
                                                          ctctCGTGCCGACCAAAAATACCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcacatc                                        >  2:2401179/1‑98 (MQ=255)
                                                           tctcGTGCCGACCAAAAATCCCAATAAAAAACCAACCCTTACGGTTGGTTTTGTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                       <  2:652927/100‑1 (MQ=255)
                                                           tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGCTAAACAGACCGTGACACATCACAGcc                                       <  1:2408075/100‑1 (MQ=255)
                                                           tctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                       <  1:3029838/100‑1 (MQ=255)
                                                             tcGTGCCGATCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                       <  2:1876223/98‑1 (MQ=255)
                                                                  ccGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                >  1:2752609/1‑100 (MQ=255)
                                                                   cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                               >  1:2217360/1‑100 (MQ=255)
                                                                    gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                              >  2:1968460/1‑100 (MQ=255)
                                                                    gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGAAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                              >  2:1433594/1‑100 (MQ=255)
                                                                     aCCAAAAATACCAAGAAAAAACCAACCCTTGCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                             <  2:25879/100‑1 (MQ=255)
                                                                       cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                           >  1:505943/1‑100 (MQ=255)
                                                                        aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                          <  2:461221/100‑1 (MQ=255)
                                                                        aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGGCTGTGATGTGTCACGGTCTGTTTAtcgaat                           <  2:1962754/99‑6 (MQ=18)
                                                                        aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                           <  2:2197721/99‑1 (MQ=255)
                                                                         aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                         <  1:2195337/100‑1 (MQ=255)
                                                                              cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                    <  1:3529817/100‑1 (MQ=255)
                                                                               ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   >  2:3769745/1‑100 (MQ=255)
                                                                               ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   >  1:3472233/1‑100 (MQ=255)
                                                                               ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                   <  2:3939221/100‑1 (MQ=255)
                                                                               ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCAATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                    >  2:180337/1‑99 (MQ=255)
                                                                                tataacaaaaCCCACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGGTAAACAGACCGTGCAGAATCACCGCCTGTTTATTGTCTGTTATc                     >  2:881920/7‑97 (MQ=19)
                                                                                tataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  2:2928759/5‑100 (MQ=38)
                                                                                cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                  >  1:3628245/1‑100 (MQ=255)
                                                                                 ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                 >  1:3255982/4‑100 (MQ=38)
                                                                                  taaaaaaaCCAACCGTACGGGTTGGATTTTTTATATCTGCAATTAATTCAATAAACAGACAGTGACACATGACAGCATGTTTATTTTCTGTTATCAGaa                 >  2:2861430/3‑99 (MQ=37)
                                                                                  taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  1:499919/3‑100 (MQ=38)
                                                                                  taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  1:1039645/3‑100 (MQ=38)
                                                                                  taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                >  2:2742692/3‑100 (MQ=38)
                                                                                   gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  1:1042042/1‑100 (MQ=255)
                                                                                   gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  1:2737444/1‑100 (MQ=255)
                                                                                   aaaaaaaCCAACCGTAGGGGTTGGTTTTTTTATATCTGCAATTAACTCGATAAACAGACCGTGCCACATCAAAACCTGTTTATTCTCAGTTATCAGAAc                >  2:1745326/2‑99 (MQ=25)
                                                                                   aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  1:1168872/2‑100 (MQ=38)
                                                                                   aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg               >  1:3698783/2‑100 (MQ=38)
                                                                                    aaacaaCCAACCGTAAGGGTTGGTTTTTTTATATCTACAATTAATTCGATAAACAGACCGTGACACATCCCAGCCTGTTTATTTTCTGTTATCAAAAc                >  2:166611/1‑98 (MQ=37)
                                                                                     aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc             >  1:2082438/1‑100 (MQ=39)
                                                                                     aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc             >  1:3744676/1‑100 (MQ=255)
                                                                                      aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc            >  1:2199013/1‑100 (MQ=255)
                                                                                       aaaCCAACCCTTACGGTTGGTTTTTTTTTTTTTTTTTTTAATTCGATAAACAGACGGTGACACATCACATCTTGTTTATTTTCTGTTATCATAACCTcaa           >  2:961774/1‑98 (MQ=255)
                                                                                       aaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATTACAGCCTCTTTATTTTCTGTTATCACAACGTCCa           >  2:2226638/1‑100 (MQ=255)
                                                                                         accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg          <  1:3413681/99‑1 (MQ=255)
                                                                                          ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg          <  2:3338528/98‑1 (MQ=255)
                                                                                          ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc        >  1:2549124/1‑100 (MQ=255)
                                                                                          ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc        >  1:3816076/1‑100 (MQ=255)
                                                                                          ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc        <  1:742553/100‑1 (MQ=255)
                                                                                             agccTTACGGTTGATTTTTTTATATCTGCAATTAATTCGATAGACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac     >  1:2671551/3‑100 (MQ=255)
                                                                                             accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac     <  1:2888484/100‑1 (MQ=255)
                                                                                                |                                                                                                   
ACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACC  >  NC_000913/2286215‑2286410

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