Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,313 T→A 18.8% intergenic (+26/+77) 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,3130TA18.8% 252.3 / 54.9 119intergenic (+26/+77)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 T (43/52);  new base A (15/7);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 6.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

Only 100 of 122 total aligned reads displayed.

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCC  >  NC_000913/2286217‑2286411
                                                                                                |                                                                                                  
gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                                 >  2:451410/1‑100 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCATAAGTCCCAAGAAAAAACCAACCCTTACggttgg                                                                                            <  1:1920821/100‑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)
         ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                        <  2:2907873/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                     <  2:1235425/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                    <  2:892438/100‑3 (MQ=38)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAAACAACCCTTACGGTAGGTTTTTTta                                                                                    >  1:2699824/1‑100 (MQ=255)
              tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                   <  2:2639395/100‑1 (MQ=255)
               aacaaatagTCAAGGGGTGTCGGTGATAGGAGGTTGAAATCCTGTCGTGTCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                        <  1:1123517/86‑1 (MQ=255)
               aGCGCACCGTCATTGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                  <  2:2376162/100‑5 (MQ=37)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                  >  1:261238/1‑100 (MQ=255)
                 cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATc                                                                                <  2:2572477/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTCGGGTGTCGGAGGTTCAAATCCTCTCGTGCCGATCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttg                                                                                  <  1:1125490/95‑5 (MQ=37)
                    aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                             <  2:1432470/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTAGGGGGTCGGAGGTTCAAATCGTCTCGTGCCGGCCAAATATCCCAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTtata                                                                                  <  1:446233/95‑1 (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)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                    <  2:1732130/100‑1 (MQ=255)
                             gggTGTAGGAGGTCGGAGGTTCACATCCTCGCGTGCCGTCTAAAAATCCCAAGAAATAACCACCCCTTGCGGTTGGTTTTTTTAGATCTGCAATTAATTc                                                                    <  1:828085/100‑1 (MQ=255)
                              ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                   <  1:2122744/100‑1 (MQ=255)
                                tgtCGGGAGTCGTATGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                                 <  1:1046528/100‑1 (MQ=255)
                                  tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaa                                                               >  2:214162/1‑100 (MQ=255)
                                  tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaa                                                               <  2:941986/100‑1 (MQ=255)
                                    gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                             >  2:2207149/1‑100 (MQ=255)
                                    atctgTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                             <  1:1433950/96‑1 (MQ=255)
                                       ggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                          >  1:1969658/1‑100 (MQ=255)
                                        gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGat                                                         >  2:1433955/1‑99 (MQ=255)
                                        gTCGGAGGTTAAAATACTCTCGTGCCGACCAAAAATCCCAAGAAAAAAACAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                            <  1:2321623/97‑1 (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)
                                                 tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                                <  2:1595976/100‑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)
                                                            cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                     <  2:1604426/100‑1 (MQ=255)
                                                             ttGCAGATCAAAAATCCCAAGAAAAAACTAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAGACAGACCGTGACACATCACAGCCTGtt                                    <  1:1218753/99‑1 (MQ=255)
                                                              tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttt                                   <  2:103722/100‑1 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                <  2:541974/100‑1 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                >  2:1911092/1‑100 (MQ=255)
                                                                    ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                             >  2:792875/1‑100 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                           >  2:1233479/1‑100 (MQ=255)
                                                                       aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                          <  2:1312577/100‑1 (MQ=255)
                                                                          cagatataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                         >  1:1634901/9‑98 (MQ=25)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCTATTAATTCGATAAACAGACCGTGACACATAACAGCCTGTTTATTTTCTGTTATCa                     >  2:84327/1‑100 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                     <  1:1857963/100‑1 (MQ=255)
                                                                              tataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                   >  2:1098468/5‑100 (MQ=38)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                   >  2:822216/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                   >  2:2347298/1‑100 (MQ=255)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                  >  2:1752250/4‑100 (MQ=38)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                  >  2:2509407/4‑100 (MQ=38)
                                                                               aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                  <  2:1511325/100‑1 (MQ=255)
                                                                               aaGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                  >  2:2496934/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)
                                                                                  aaacaaCCAACCCTCACGGTTGGTTTTTTTCTATCCGCAATTCATTCGATAAACAGACCGTGACACATCCCAGCCTGTTTATTTTCCGCTATccgca                  >  1:2007782/1‑93 (MQ=255)
                                                                                   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)
                                                                                    aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACCCATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc             >  1:821259/1‑100 (MQ=38)
                                                                                    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           >  2:2188479/1‑100 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg           >  1:2892456/1‑100 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg           >  2:2188630/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)
                                                                                          aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca       <  1:2118097/100‑1 (MQ=255)
                                                                                            cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACGGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaca     <  1:2080568/100‑1 (MQ=255)
                                                                                            cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaca     >  2:1817766/1‑100 (MQ=255)
                                                                                             ccTTACGGTAAGCTTCTTCAGATCTGCAATGAATTCTATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCacac    <  2:1697367/100‑1 (MQ=255)
                                                                                               ttACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac      <  1:2633672/96‑1 (MQ=255)
                                                                                               ttACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc  <  2:545725/100‑1 (MQ=255)
                                                                                               ttACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAccc  <  2:136233/100‑1 (MQ=255)
                                                                                                |                                                                                                  
GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCC  >  NC_000913/2286217‑2286411

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: