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

GTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  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)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttc                                                                                   <  1:1382482/100‑2 (MQ=38)
            ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                   <  1:2025092/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                  <  1:1020034/100‑1 (MQ=255)
             gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                  <  1:566385/100‑1 (MQ=255)
              tAGCGCACCGTCATGGGGTGTCGGGGATCGGAGGTTAAAATCCTCTCGTGACGACCAAAAATCCCAAAAAAAAACCAACCGTAAGggttggtt                                                                                         <  2:949953/93‑1 (MQ=37)
              tAGCGCACCGTCATGGGGTGTCGGAGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                 <  2:606070/100‑1 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                <  1:779166/100‑1 (MQ=255)
               aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                <  1:1386726/100‑1 (MQ=255)
                   cACCGTAATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCAACCAAAAATCCCAATAAAAAACCAACCCTTACGGTTGGttttttt                                                                                   <  2:789352/93‑1 (MQ=255)
                   cACCGAAATGAGGTATAGGGGGTCGGAGGTTCAAATCCTAAAAAAAAAAAAAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                               <  2:1009228/97‑6 (MQ=2)
                     ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                          <  1:1614039/100‑1 (MQ=255)
                      cGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttggg                                                                              <  2:1196938/95‑7 (MQ=37)
                       gTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                        <  1:1012242/100‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                       <  1:1338251/100‑1 (MQ=255)
                        tCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                       <  1:538801/100‑1 (MQ=255)
                             gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                   <  2:1631304/99‑1 (MQ=255)
                                tgtCGAGGGTCGGAGGTTAAAATCCTCTCGTGCCGGCCGAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                               <  2:1391328/100‑1 (MQ=255)
                                 gtaggaggTCGGAGGTTAAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                              <  2:127819/94‑1 (MQ=255)
                                 gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACGAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                              <  2:74186/100‑1 (MQ=255)
                                 gtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATa                                                              >  1:1171850/1‑100 (MQ=255)
                                        gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAc                                                       >  2:169606/1‑100 (MQ=255)
                                          cGGAGGTTCCAATCCTCTCGTAACGACCAAAAATCCCAAGAAAAAACCAACTCTTACGGTTGATATTTTTATATCTGCAATTAATTCGATAAACAGAcag                                                     >  2:2086788/1‑98 (MQ=255)
                                          cGGAGGTTCAAATCCTCTCGTGCCGACCAAAAACCCCAAGAAAAAACAAACCATTACGGCTGGTTTTTTTATATAAGCAAATAATTCGATAAACAGAc                                                       >  2:478507/1‑98 (MQ=255)
                                           ggAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCTTTACGGTTGGTTTTTTTTATATCTGCAATTAATTCGCTAAACAGACCg                                                     >  1:1251517/1‑100 (MQ=255)
                                            gAGGTTCAAATCCTCTCGTGTCGACGAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTCTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                   <  2:1055596/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                 <  1:1695908/100‑1 (MQ=255)
                                              ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                 <  1:1675981/100‑1 (MQ=255)
                                                ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                               >  2:654657/1‑100 (MQ=255)
                                                 tCAAATCCTCTCGTGCCGACCAAATATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATAACTGCAATTAATTCGATAAACAGACCGTGAcaca                                              <  2:1669168/100‑1 (MQ=255)
                                                 tCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGAcaca                                              <  1:1956449/100‑1 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAACACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                             >  1:1806512/1‑100 (MQ=255)
                                                  cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAt                                             >  1:1574163/1‑100 (MQ=255)
                                                    aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATca                                           >  2:1036329/1‑100 (MQ=255)
                                                       ccTCTCGTGCCGACCAAAAATCCCAAGAAAATACCAACCCTTACGGCTGGTTTTTTTATATCTGCAATGAATTCGATAAACACACCGTAACACATca                                           >  2:1198354/1‑97 (MQ=255)
                                                           tcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                    <  1:303187/100‑1 (MQ=255)
                                                             gTGCCGACCCAAAAACCCAGGAAAAAACCAACCCTCACGCTTGGATTCTTTACACCCGCAATTAATTCGATAAACAGACCGTGACCCATcccatc                                       >  2:2148084/1‑93 (MQ=255)
                                                              tGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttt                                 >  2:1175958/1‑100 (MQ=255)
                                                               gCCGACCAAAAATCCCAAGAAAAAACCAACCTTTACGGTTGGTTTCTATTATATCTACAATTCATTCGATAAACAGACCGTGACACATCACAGc                                       >  2:2092871/1‑94 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                             >  1:2344072/1‑100 (MQ=255)
                                                                  gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTCTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                             <  1:1910681/100‑1 (MQ=255)
                                                                     cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                          <  1:2054233/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTCCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                         <  2:39348/100‑1 (MQ=255)
                                                                      aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                         <  1:2451441/100‑1 (MQ=255)
                                                                       aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                        <  1:2351391/100‑1 (MQ=255)
                                                                        aaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGtt                       <  2:548231/100‑1 (MQ=255)
                                                                         attcacAAGAAAAAACAATCCCTTACGGTTGGTTTTTGTAGATCTGCAATTAATTCGATAAACCGACCGTGACACATCACAGCCTGTTTATTTTCTGt                        <  2:2340342/93‑1 (MQ=255)
                                                                          aTCCCAAGAAAAAACCAACTCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACCTCACAGCCTGTTTATTTTCTGTTa                      <  2:30080/99‑1 (MQ=255)
                                                                            gatataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                   >  1:2307535/7‑100 (MQ=37)
                                                                            cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCa                   <  1:2526951/100‑1 (MQ=255)
                                                                             ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAg                  >  1:1047687/1‑100 (MQ=255)
                                                                              cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                 <  1:318482/100‑1 (MQ=255)
                                                                               ataaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                >  1:2309336/4‑100 (MQ=38)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACGGACTGTGCCATATGCCCGCCTGTTTATTTTCTGTTATCAGAAc               >  2:511241/3‑100 (MQ=25)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTTAAACATCACAGCCTTTTTATTTTCTGTTATCAGAAc               >  2:392539/3‑100 (MQ=37)
                                                                                taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  2:147808/3‑100 (MQ=38)
                                                                                aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  1:300296/1‑100 (MQ=255)
                                                                                aGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTAACACATCACAGCCTGTTTATTTTCTGTTATCAGAAc               >  1:1803475/1‑100 (MQ=255)
                                                                                 gAAAAAGCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGGTAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:123283/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:328503/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:2258883/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:1471388/100‑1 (MQ=255)
                                                                                 gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1519182/1‑100 (MQ=255)
                                                                                 aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:2411843/2‑100 (MQ=38)
                                                                                 aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1216158/2‑100 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACCGCCTGTATATTTTCTGTTATCAGAACGt             >  2:1197240/1‑99 (MQ=38)
                                                                                   aaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:1718281/1‑100 (MQ=39)
                                                                                   aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  2:2082645/1‑100 (MQ=255)
                                                                                   aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:2317250/1‑100 (MQ=255)
                                                                                   aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:965766/1‑100 (MQ=255)
                                                                                    aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           >  2:362885/1‑100 (MQ=39)
                                                                                    aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc           <  1:1909703/100‑1 (MQ=255)
                                                                                       accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCCGAACGTCCa          >  2:1781155/1‑98 (MQ=255)
                                                                                        ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAc       <  1:448727/100‑1 (MQ=255)
                                                                                         caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc      <  2:922130/100‑1 (MQ=255)
                                                                                         caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc      <  1:1674829/100‑1 (MQ=255)
                                                                                          aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca     <  1:1331277/100‑1 (MQ=255)
                                                                                          aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAAGTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTAACAGAACGTCCAGACca     <  2:1016651/100‑1 (MQ=255)
                                                                                             ccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcacan  <  2:1683491/100‑2 (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: