Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,315 C→G 15.9% 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.9% 232.0 / 34.3 87intergenic (+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 (27/46);  new base G (10/4);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  NC_000913/2286217‑2286409
                                                                                                  |                                                                                                
gTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                  >  1:2479867/1‑99 (MQ=255)
 ttgcgcAGCCTGGAAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgg                                                                                                >  1:1737972/3‑100 (MQ=255)
 tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACgg                                                                                                >  2:1338242/1‑100 (MQ=255)
 tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACAAAAAATCCCAAGACAAAACCAACCCTTACgg                                                                                                >  1:712726/1‑100 (MQ=255)
   gcgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                              <  2:2396926/100‑1 (MQ=255)
    cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                             >  1:922138/1‑100 (MQ=255)
    cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                             <  1:2444024/100‑1 (MQ=255)
    cgcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                             >  1:1807026/1‑100 (MQ=255)
     gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttgg                                                                                            <  1:2071758/100‑1 (MQ=255)
       aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                             <  2:2046735/97‑1 (MQ=255)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGG‑ttt‑tttc                                                                                   <  1:1382482/100‑2 (MQ=38)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑ttt‑tttt                                                                                   <  1:2025092/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTta                                                                                  <  1:1020034/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTta                                                                                  <  1:566385/100‑1 (MQ=255)
              tAGCGCACCGTCATGGGGTGTCGGGGATCGGAGGTTAAAATCCTCTCGTGACGACCAAAAATCCCAAAAAAAAACCAACCGTAAGggttgg‑tt                                                                                         <  2:949953/93‑1 (MQ=37)
              tAGCGCACCGTCATGGGGTGTCGGAGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTtat                                                                                 <  2:606070/100‑1 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTtata                                                                                <  1:779166/100‑1 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTtata                                                                                <  1:1386726/100‑1 (MQ=255)
                   cACCGTAATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCAACCAAAAATCCCAATAAAAAACCAACCCTTACGGTTGG‑ttt‑tttt                                                                                   <  2:789352/93‑1 (MQ=255)
                   cACCGAAATGAGGTATAGGGGGTCGGAGGTTCAAATCCTAAAAAAAAAAAAAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGG‑ttt‑tttcttgg                                                                               <  2:1009228/97‑6 (MQ=2)
                     ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCa                                                                          <  1:1614039/100‑1 (MQ=255)
                      cGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGG‑ttt‑tttcttggg                                                                              <  2:1196938/95‑7 (MQ=37)
                       gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCaat                                                                        <  1:1012242/100‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCaatt                                                                       <  1:1338251/100‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCaatt                                                                       <  1:538801/100‑1 (MQ=255)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCaattaatt                                                                   <  2:1631304/99‑1 (MQ=255)
                                tgtCGAGGGTCGGAGGTTAAAATCCTCTCGTGCCGGCCGAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGAt                                                               <  2:1391328/100‑1 (MQ=255)
                                 gtaggaggTCGGAGGTTAAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATa                                                              <  2:127819/94‑1 (MQ=255)
                                 gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACGAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATa                                                              <  2:74186/100‑1 (MQ=255)
                                 gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATa                                                              >  1:1171850/1‑100 (MQ=255)
                                        gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGAc                                                       >  2:169606/1‑100 (MQ=255)
                                          cGGAGGTTCCAATCCTCTCGTAACGACCAAAAATCCCAAGAAAAAACCAACTCTTACGGTTGA‑TAT‑TTTTATATCTGCAATTAATTCGATAAACAGAcag                                                     >  2:2086788/1‑98 (MQ=255)
                                          cGGAGGTTCAAATCCTCTCGTGCCGACCAAAAACCCCAAGAAAAAACAAACCATTACGGCTGG‑TTT‑TTTTATATAAGCAAATAATTCGATAAACAGAc                                                       >  2:478507/1‑98 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCTTTACGGTTGGTTTT‑TTTTATATCTGCAATTAATTCGCTAAACAGACCg                                                     >  1:1251517/1‑100 (MQ=255)
                                            gAGGTTCAAATCCTCTCGTGTCGACGAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TCT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                   <  2:1055596/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                 <  1:1695908/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                 <  1:1675981/100‑1 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                               >  2:654657/1‑100 (MQ=255)
                                                 tCAAATCCTCTCGTGCCGACCAAATATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATAACTGCAATTAATTCGATAAACAGACCGTGAcaca                                              <  2:1669168/100‑1 (MQ=255)
                                                 tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                              <  1:1956449/100‑1 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAACACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                             >  1:1806512/1‑100 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                             >  1:1574163/1‑100 (MQ=255)
                                                    aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATca                                           >  2:1036329/1‑100 (MQ=255)
                                                       ccTCTCGTGCCGACCAAAAATCCCAAGAAAATACCAACCCTTACGGCTGG‑TTT‑TTTTATATCTGCAATGAATTCGATAAACACACCGTAACACATca                                           >  2:1198354/1‑97 (MQ=255)
                                                           tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                    <  1:303187/100‑1 (MQ=255)
                                                             gTGCCGACCCAAAAACCCAGGAAAAAACCAACCCTCACGCTTGG‑ATT‑CTTTACACCCGCAATTAATTCGATAAACAGACCGTGACCCATcccatc                                       >  2:2148084/1‑93 (MQ=255)
                                                              tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttt                                 >  2:1175958/1‑100 (MQ=255)
                                                               gCCGACCAAAAATCCCAAGAAAAAACCAACCTTTACGGTTGG‑TTTCTATTATATCTACAATTCATTCGATAAACAGACCGTGACACATCACAGc                                       >  2:2092871/1‑94 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                             >  1:2344072/1‑100 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑CTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                             <  1:1910681/100‑1 (MQ=255)
                                                                     cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                          <  1:2054233/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTCCGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                         <  2:39348/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                         <  1:2451441/100‑1 (MQ=255)
                                                                       aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                        <  1:2351391/100‑1 (MQ=255)
                                                                        aaaTCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                       <  2:548231/100‑1 (MQ=255)
                                                                         attcacAAGAAAAAACAATCCCTTACGGTTGG‑TTT‑TTGTAGATCTGCAATTAATTCGATAAACCGACCGTGACACATCACAGCCTGTTTATTTTCTGt                        <  2:2340342/93‑1 (MQ=255)
                                                                          aTCCCAAGAAAAAACCAACTCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACCTCACAGCCTGTTTATTTTCTGTTa                      <  2:30080/99‑1 (MQ=255)
                                                                            gatataaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                   >  1:2307535/7‑100 (MQ=37)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                   <  1:2526951/100‑1 (MQ=255)
                                                                             ccAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                  >  1:1047687/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                 <  1:318482/100‑1 (MQ=255)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                >  1:2309336/4‑100 (MQ=38)
                                                                                taaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACGGACTGTGCCATATGCCCGCCTGTTTATTTTCTGTTATCAGAAc               >  2:511241/3‑100 (MQ=25)
                                                                                taaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTTAAACATCACAGCCTTTTTATTTTCTGTTATCAGAAc               >  2:392539/3‑100 (MQ=37)
                                                                                taaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  2:147808/3‑100 (MQ=38)
                                                                                aGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  1:300296/1‑100 (MQ=255)
                                                                                aGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTAACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  1:1803475/1‑100 (MQ=255)
                                                                                 gAAAAAGCCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGGTAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:123283/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:328503/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:2258883/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:1471388/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1519182/1‑100 (MQ=255)
                                                                                 aaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:2411843/2‑100 (MQ=38)
                                                                                 aaaaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1216158/2‑100 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACCGCCTGTATATTTTCTGTTATCAGAACGt             >  2:1197240/1‑99 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:1718281/1‑100 (MQ=39)
                                                                                   aaaaaCCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  2:2082645/1‑100 (MQ=255)
                                                                                   aaaaaCCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:2317250/1‑100 (MQ=255)
                                                                                   aaaaaCCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:965766/1‑100 (MQ=255)
                                                                                    aaaaCCAACCGTAAGGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           >  2:362885/1‑100 (MQ=39)
                                                                                    aaaaCCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           <  1:1909703/100‑1 (MQ=255)
                                                                                       accaaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCCGAACGTCCa          >  2:1781155/1‑98 (MQ=255)
                                                                                        ccaaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc       <  1:448727/100‑1 (MQ=255)
                                                                                         caaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc      <  2:922130/100‑1 (MQ=255)
                                                                                         caaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc      <  1:1674829/100‑1 (MQ=255)
                                                                                          aaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca     <  1:1331277/100‑1 (MQ=255)
                                                                                          aaccCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAAGTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACAGAACGTCCAGACca     <  2:1016651/100‑1 (MQ=255)
                                                                                             ccTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcacan  <  2:1683491/100‑2 (MQ=255)
                                                                                                  |                                                                                                
GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTT‑TTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  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: ‑