Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 16.5% intergenic (+24/+79) proL → / ← yejO tRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,3110CG16.5% 221.6 / 37.2 87intergenic (+24/+79)proL/yejOtRNA‑Pro/pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (25/45);  new base G (10/4);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  NC_000913/2286212‑2286409
                                                                                                   |                                                                                                    
ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAAccc                                                                                                      >  2:1922162/1‑100 (MQ=255)
     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:1386726/100‑1 (MQ=255)
                    aGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                <  1:779166/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:1471388/100‑1 (MQ=255)
                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:328503/100‑1 (MQ=255)
                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1519182/1‑100 (MQ=255)
                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              <  1:2258883/100‑1 (MQ=255)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:1216158/2‑100 (MQ=38)
                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg              >  1:2411843/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:965766/1‑100 (MQ=255)
                                                                                        aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc            >  1:2317250/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      <  1:1674829/100‑1 (MQ=255)
                                                                                              caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAcc      <  2:922130/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)
                                                                                                   |                                                                                                    
GGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAC  >  NC_000913/2286212‑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: