Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 18.0% 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,3150CG18.0% 298.5 / 60.9 122intergenic (+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 (44/56);  new base G (15/7);  total (59/63)
Fisher's exact test for biased strand distribution p-value = 5.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 124 total aligned reads displayed.

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCC  >  NC_000913/2286217‑2286414
                                                                                                  |                                                                                                   
  aGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggt                                                                                                  >  2:1407696/1‑100 (MQ=255)
   gcgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAAGTTCAAATCCTCTCGTGCCGACCAAAAAAAAAAAAAAAAAACCAACCCTTAc                                                                                                     <  1:2061364/96‑1 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCATAAGTCCCAAGAAAAAACCAACCCTTACggttgg                                                                                               <  1:1920821/100‑1 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                               <  2:1767293/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)
        gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttt                                                                                            <  1:2017038/100‑1 (MQ=39)
         ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                           <  2:2907873/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                        <  2:1235425/100‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttn                                                                                        <  2:1179249/100‑2 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                       <  2:892438/100‑3 (MQ=38)
              tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                      <  2:2639395/100‑1 (MQ=255)
               aacaaatagTCAAGGGGTGTCGGTGATAGGAGGTTGAAATCCTGTCGTGTCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                           <  1:1123517/86‑1 (MQ=255)
               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)
                    aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                                <  2:2588539/100‑1 (MQ=255)
                    aCCGTCATGGGGTGTAGGGGGTCGGAGGTTCAAATCGTCTCGTGCCGGCCAAATATCCCAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTtata                                                                                     <  1:446233/95‑1 (MQ=255)
                     ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                              >  2:1687353/1‑100 (MQ=255)
                       gTCATGGGGTGTCGGGAATCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                              <  1:1084610/99‑1 (MQ=255)
                         cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatta                                                                           >  1:40364/1‑100 (MQ=255)
                          aTGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaa                                                                          >  2:1074508/1‑100 (MQ=255)
                           tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                         >  1:2125256/1‑100 (MQ=255)
                             gggTGTAGGAGGTCGGAGGTTCACATCCTCGCGTGCCGTCTAAAAATCCCAAGAAATAACCACCCCTTGCGGTTGGTTTTTTTAGATCTGCAATTAATTc                                                                       <  1:828085/100‑1 (MQ=255)
                              ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                      <  1:2122744/100‑1 (MQ=255)
                              ggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                      >  2:2122757/1‑100 (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)
                                    atctgTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                <  1:1433950/96‑1 (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)
                                            gAGGTTCAAATACTCTCGTGCCGACCAAAATTCCAAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                <  1:110300/92‑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)
                                                  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:1911092/1‑100 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                   >  2:364523/1‑100 (MQ=255)
                                                                 cGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtt                                   <  2:541974/100‑1 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                  <  2:2325649/100‑1 (MQ=255)
                                                                    ccAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                <  2:66107/100‑1 (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)
                                                                           tCCCAAGAAAAAACCAACCCTTTCGGTTTGTTTTTTTATATCTACAATTATTTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                         <  1:1368115/100‑1 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCTATTAATTCGATAAACAGACCGTGACACATAACAGCCTGTTTATTTTCTGTTATCa                        >  2:84327/1‑100 (MQ=255)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                        <  1:1857963/100‑1 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                      >  2:2347298/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                      >  2:822216/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                      <  1:2654833/100‑1 (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)
                                                                                aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAAAAc                    >  1:1646304/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)
                                                                                   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)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa               >  1:1307664/1‑99 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCa               >  1:2293761/1‑99 (MQ=39)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg              >  2:2188630/1‑100 (MQ=255)
                                                                                      aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg              >  1:2892456/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)
                                                                                                 aCGGTCGGTTTTTTTAAATCCGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAGCAGAACGTCCAGACCACAcc      <  1:2320380/97‑1 (MQ=255)
                                                                                                  cGGTTGGTTTTTTTATATCTGAAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc  >  1:1429231/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCC  >  NC_000913/2286217‑2286414

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