Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 15.1% intergenic (+28/+75) 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,3150CG15.1% 423.0 / 58.7 152intergenic (+28/+75)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/60);  new base G (15/8);  total (84/68)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 155 total aligned reads displayed.

CGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCC  >  NC_000913/2286216‑2286414
                                                                                                   |                                                                                                   
cGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                     >  1:682130/1‑100 (MQ=255)
  tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgg                                                                                                   <  1:487051/100‑1 (MQ=255)
     cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                                <  1:3059582/100‑1 (MQ=255)
      gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                               >  2:3887768/1‑100 (MQ=255)
        aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                             <  1:2490066/100‑1 (MQ=255)
        aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAAAAAAAAAAAAAAACCAACCCTTACggttg                                                                                                <  2:2822194/97‑1 (MQ=255)
           cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                          >  1:1772439/1‑100 (MQ=255)
            tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttttt                                                                                         <  2:2166632/100‑1 (MQ=255)
            tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                         >  1:2023417/1‑100 (MQ=255)
             ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                        >  2:3383134/1‑100 (MQ=255)
             ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                        >  1:3381950/1‑100 (MQ=255)
             ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                         <  1:1338835/99‑1 (MQ=255)
             ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                         <  2:2966991/99‑1 (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)
                aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                     <  1:1302151/100‑1 (MQ=255)
                 gcgaACCGTCATGGGGTGTCGGGGGTCAGAGGTTCAAATCCGATCGTGCTGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                     <  2:753240/99‑5 (MQ=25)
                  cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                                   >  1:3545287/1‑100 (MQ=255)
                  cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                                   >  1:1549261/1‑100 (MQ=255)
                   gcACCGTCATGGGGTGTCGGGGGTTGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCt                                                                                  <  1:3124538/100‑1 (MQ=255)
                     aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtaagatcgg                                                                                >  2:240882/1‑93 (MQ=255)
                       cGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                              <  1:177358/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:1840536/100‑1 (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)
                                      ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCGTACTGTTGGTGTTTTTATATCTGCAATTAATTCGATAAACa                                                               >  2:325583/1‑100 (MQ=255)
                                      gaggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                               <  1:84742/98‑1 (MQ=255)
                                         gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAc                                                            >  2:3910096/1‑100 (MQ=255)
                                           cGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                                          >  1:1633601/1‑100 (MQ=255)
                                           cGGAGGTTCAAATACTCTCGTGCCGACCAAAGATCCCAAAAAAAAAACAACCCTAACGGTTGGTTGTTTTATATCTGCAATTAATTCGATAAAc                                                                <  2:2199069/94‑1 (MQ=255)
                                             gAGGTTCAAATCCTCTCGTGCCGACCAAACATCTCAGAAAAAACCACCCCTTACGCTTAGTCTTCTTATAGCTGCCATTAATACGATACACCGACCCTg                                                        <  2:1163732/99‑1 (MQ=255)
                                              aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                       >  1:2166966/1‑100 (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)
                                                   cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                                  >  2:3057092/1‑100 (MQ=255)
                                                   cAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGaca                                                     <  2:3162096/97‑1 (MQ=255)
                                                    aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                                 <  1:870750/100‑1 (MQ=255)
                                                       tCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                              >  1:2201042/1‑100 (MQ=255)
                                                        ccTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                             <  2:200762/100‑1 (MQ=255)
                                                        ccTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGATTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                             <  2:1725488/100‑1 (MQ=255)
                                                         ctctCGTGCCGACCAAAAATACCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcacatc                                            >  2:2401179/1‑98 (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)
                                                                    aCCAAAAATACCAAGAAAAAACCAACCCTTGCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                 <  2:25879/100‑1 (MQ=255)
                                                                      cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                               >  1:505943/1‑100 (MQ=255)
                                                                       aaaaaTCCCAATAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCTATTATTTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                              <  2:395529/100‑1 (MQ=255)
                                                                       aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                              <  2:2283934/100‑1 (MQ=255)
                                                                       aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                              <  2:461221/100‑1 (MQ=255)
                                                                       aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                              >  2:661915/1‑100 (MQ=255)
                                                                       aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATCCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                               <  2:2158406/99‑1 (MQ=255)
                                                                       aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGGTTTTTTATATCTGCAATTAATTCGATAAACAGGCTGTGATGTGTCACGGTCTGTTTAtcgaat                               <  2:1962754/99‑6 (MQ=18)
                                                                        aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                             <  1:222243/100‑1 (MQ=255)
                                                                        aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                             <  1:2195337/100‑1 (MQ=255)
                                                                          aaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTa                           <  1:1676442/100‑1 (MQ=255)
                                                                           aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGCTAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAt                          >  2:1560894/1‑100 (MQ=255)
                                                                             cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                         <  1:1133679/99‑1 (MQ=255)
                                                                             cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                        <  1:3529817/100‑1 (MQ=255)
                                                                              ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                       >  2:3769745/1‑100 (MQ=255)
                                                                              ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                       >  1:3472233/1‑100 (MQ=255)
                                                                              ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCAATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                        >  2:180337/1‑99 (MQ=255)
                                                                               tataacaaaaCCCACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGGTAAACAGACCGTGCAGAATCACCGCCTGTTTATTGTCTGTTATc                         >  2:881920/7‑97 (MQ=19)
                                                                               cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                      >  1:2034752/1‑100 (MQ=255)
                                                                               cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                      >  1:3628245/1‑100 (MQ=255)
                                                                                ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                     >  1:3255982/4‑100 (MQ=38)
                                                                                aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                     >  1:2638924/1‑100 (MQ=255)
                                                                                 taaaaaaaCCAACCGTACGGGTTGGATTTTTTATATCTGCAATTAATTCAATAAACAGACAGTGACACATGACAGCATGTTTATTTTCTGTTATCAGaa                     >  2:2861430/3‑99 (MQ=37)
                                                                                 taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                    >  1:499919/3‑100 (MQ=38)
                                                                                 taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc                    >  1:1039645/3‑100 (MQ=38)
                                                                                  gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg                   >  1:1042042/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)
                                                                                      aaaCCAACCCTTACGGTTGGTTTTTTTTTTTTTTTTTTTAATTCGATAAACAGACGGTGACACATCACATCTTGTTTATTTTCTGTTATCATAACCTcaa               >  2:961774/1‑98 (MQ=255)
                                                                                      aaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATTACAGCCTCTTTATTTTCTGTTATCACAACGTCCa               >  2:2226638/1‑100 (MQ=255)
                                                                                        accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg              <  1:3413681/99‑1 (MQ=255)
                                                                                         ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc            >  1:2549124/1‑100 (MQ=255)
                                                                                         ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc            >  1:3816076/1‑100 (MQ=255)
                                                                                         ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc            <  1:742553/100‑1 (MQ=255)
                                                                                          caaccCTTACGGTTGGTTTTTTTTTTTTTTTTTTTAATTAAATAAACATACCGTGTCACATTTTAGCCTGTTTTTTTTCTGTTGTcttttcg                   >  2:2654107/1‑86 (MQ=255)
                                                                                            agccTTACGGTTGATTTTTTTATATCTGCAATTAATTCGATAGACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac         >  1:2671551/3‑100 (MQ=255)
                                                                                            accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac         <  1:2888484/100‑1 (MQ=255)
                                                                                               cTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAcc      >  1:2724392/1‑100 (MQ=255)
                                                                                                   cGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc  <  1:2544981/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCC  >  NC_000913/2286216‑2286414

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: