Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 19.8% 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,3110CG19.8% 265.2 / 63.3 116intergenic (+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 (43/50);  new base G (15/8);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 119 total aligned reads displayed.

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  NC_000913/2286217‑2286409
                                                                                              |                                                                                                  
gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                               >  2:451410/1‑100 (MQ=255)
  aGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggt                                                                                             >  2:1407696/1‑100 (MQ=255)
   gcgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAAGTTCAAATCCTCTCGTGCCGACCAAAAAAAAAAAAAAAAAACCAACCCTTAc                                                                                                <  1:2061364/96‑1 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggt                                                                                             <  1:1244895/97‑1 (MQ=255)
      cAGCCTGGTAGCGCACCGTCATGGTGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggt                                                                                         <  2:1028405/100‑1 (MQ=255)
      cAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttggt                                                                                         >  2:469461/1‑100 (MQ=255)
       aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGggttggtt                                                                                        <  2:1017851/100‑1 (MQ=255)
        gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                       <  1:2017038/100‑1 (MQ=39)
         ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                      <  2:703817/100‑1 (MQ=39)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttn                                                                                   <  2:1179249/100‑2 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                  <  2:892438/100‑3 (MQ=38)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAAACAACCCTTACGGTAGGTTTTTTta                                                                                  >  1:2699824/1‑100 (MQ=255)
               aacaaatagTCAAGGGGTGTCGGTGATAGGAGGTTGAAATCCTGTCGTGTCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                      <  1:1123517/86‑1 (MQ=255)
               aGCGCACCGTCATTGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                <  2:2376162/100‑5 (MQ=37)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                >  1:261238/1‑100 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                <  2:867079/100‑1 (MQ=255)
                 cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                              <  2:2572477/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTCGGGTGTCGGAGGTTCAAATCCTCTCGTGCCGATCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                <  1:1125490/95‑5 (MQ=37)
                    aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                           <  2:1432470/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                           <  2:2588539/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTAGGGGGTCGGAGGTTCAAATCGTCTCGTGCCGGCCAAATATCCCAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTtata                                                                                <  1:446233/95‑1 (MQ=255)
                     ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                         >  2:1687353/1‑100 (MQ=255)
                       gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACTAAAAATCCCAAGAAAAAACCAACCCGTACGGTTGGTTTTTTTATATCTGCaa                                                                         >  1:1530941/1‑99 (MQ=255)
                       gTCATGGGGTGTCGGGAATCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                         <  1:1084610/99‑1 (MQ=255)
                         cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                      >  2:2886674/1‑100 (MQ=255)
                         cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                      >  1:40364/1‑100 (MQ=255)
                          aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaa                                                                     >  2:1074508/1‑100 (MQ=255)
                           tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                    >  1:2125256/1‑100 (MQ=255)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                  <  2:1732130/100‑1 (MQ=255)
                             gggTGTAGGAGGTCGGAGGTTCACATCCTCGCGTGCCGTCTAAAAATCCCAAGAAATAACCACCCCTTGCGGTTGGTTTTTTTAGATCTGCAATTAATTc                                                                  <  1:828085/100‑1 (MQ=255)
                              ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                 <  1:2122744/100‑1 (MQ=255)
                                tgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                               <  2:1025308/100‑1 (MQ=255)
                                tgtCGGGAGTCGTATGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                               <  1:1046528/100‑1 (MQ=255)
                                  tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaa                                                             <  2:941986/100‑1 (MQ=255)
                                  tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaa                                                             >  2:214162/1‑100 (MQ=255)
                                    gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                           >  2:2207149/1‑100 (MQ=255)
                                       ggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                        >  1:1969658/1‑100 (MQ=255)
                                        gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGat                                                       >  2:1433955/1‑99 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                    <  2:799851/100‑1 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGt                                                    >  1:497468/1‑100 (MQ=255)
                                            naGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                   >  2:1914527/2‑100 (MQ=255)
                                            gAGGTTCAAATACTCTCGTGCCGACCAAAATTCCAAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                           <  1:110300/92‑1 (MQ=255)
                                            aaatttCAAACCCTCTCGTGCCGACAAAAAATCTCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                    <  1:2110499/79‑1 (MQ=255)
                                             aGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                  <  1:371228/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                 <  2:2880346/100‑1 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACTAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                  <  1:2353150/97‑1 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                               <  2:170840/100‑1 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACAAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                  <  1:1178265/97‑1 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                             <  2:2095182/100‑1 (MQ=255)
                                                     aTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcac                                          <  2:2783203/100‑1 (MQ=255)
                                                      tCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATcaca                                         >  2:1201352/1‑100 (MQ=255)
                                                       ccTCTCGTGCCGAACAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAg                                        <  2:2746942/100‑1 (MQ=255)
                                                          ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                     >  1:1495971/1‑100 (MQ=255)
                                                          ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                     >  1:66106/1‑100 (MQ=255)
                                                            cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                   <  2:1604426/100‑1 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                              >  2:1911092/1‑100 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                              >  2:364523/1‑100 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                             <  2:2325649/100‑1 (MQ=255)
                                                                    ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                           >  2:792875/1‑100 (MQ=255)
                                                                    ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                           <  2:66107/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                         >  2:1233479/1‑100 (MQ=255)
                                                                       aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                        <  2:1312577/100‑1 (MQ=255)
                                                                          cagatataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                       >  1:1634901/9‑98 (MQ=25)
                                                                           tCCCAAGAAAAAACCAACCCTTTCGGTTTGTTTTTTTATATCTACAATTATTTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                    <  1:1368115/100‑1 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                   <  1:1857963/100‑1 (MQ=255)
                                                                              tataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                 >  2:1098468/5‑100 (MQ=38)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                 >  2:822216/1‑100 (MQ=255)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                >  2:2509407/4‑100 (MQ=38)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                >  2:1752250/4‑100 (MQ=38)
                                                                               aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                >  2:2496934/1‑100 (MQ=255)
                                                                               aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                <  2:1511325/100‑1 (MQ=255)
                                                                                aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAAAAc               >  1:1646304/1‑100 (MQ=255)
                                                                                 gAAAAACCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:2162603/1‑100 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  2:2118110/1‑100 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGGTGGGGTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGCTTATTTTCTGTTAGCAGAACg              >  1:2702539/1‑100 (MQ=255)
                                                                                  aaagaaCCAGGCGTTACGGTTGGTTTTGTGATATCGGAAATTAATTTGATAATCAGAGCGTGGCACATCACCGCCTGTTTATTTTGTGGTATCAgg                 <  1:2401940/96‑2 (MQ=255)
                                                                                  aaacaaCCAACCCTCACGGTTGGTTTTTTTCTATCCGCAATTCATTCGATAAACAGACCGTGACACATCCCAGCCTGTTTATTTTCCGCTATccgca                >  1:2007782/1‑93 (MQ=255)
                                                                                  aaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGt             >  2:2877177/1‑100 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGGTTTTTTTATTTTTTCTATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:744939/1‑100 (MQ=37)
                                                                                   aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATGACAGCCTGTTTATTTTCTGTTATCAGAACGt             >  1:13450/1‑99 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  2:2674544/1‑100 (MQ=39)
                                                                                   aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  2:910806/1‑100 (MQ=255)
                                                                                    aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           >  2:872651/1‑100 (MQ=39)
                                                                                    aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACCTCACAGCGTGTTTATTTTCTGTTATCAGAACGTcc           >  1:1335840/1‑100 (MQ=255)
                                                                                     aaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           >  1:245682/1‑99 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa          >  1:2293761/1‑99 (MQ=39)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa          >  1:1307664/1‑99 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg         >  1:2892456/1‑100 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg         >  2:2188630/1‑100 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg         >  2:2188479/1‑100 (MQ=255)
                                                                                      aaccaaccCTTACGGTTGGTTTTTTTTTTTTTTTTTTTAATTCGATAAGCAGACCATGATATATCATAGCCTGTTTTTTTTCTGTTGTccttccgt             >  1:2887955/1‑89 (MQ=255)
                                                                                      aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa          >  1:2497690/1‑99 (MQ=255)
                                                                                      aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg         <  2:2914433/100‑1 (MQ=255)
                                                                                       accaaacCTGACGATTGGTTTTTTTATATATGCAATTAAGTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGa        <  1:2161850/100‑1 (MQ=255)
                                                                                        ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc       >  2:352207/1‑100 (MQ=255)
                                                                                        ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCAATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc       >  1:2810198/1‑100 (MQ=255)
                                                                                          aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACTTCCAGACca     >  1:978312/1‑100 (MQ=255)
                                                                                            cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACGGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaca   <  1:2080568/100‑1 (MQ=255)
                                                                                            cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaca   >  2:1817766/1‑100 (MQ=255)
                                                                                             ccTTACGGTAAGCTTCTTCAGATCTGCAATGAATTCTATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCacac  <  2:1697367/100‑1 (MQ=255)
                                                                                              |                                                                                                  
GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  NC_000913/2286217‑2286409

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: