Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 5.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,3110CG 5.5% 386.8 / 12.7 127intergenic (+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 (57/62);  new base G (4/3);  total (62/65)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 129 total aligned reads displayed.

TCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGC  >  NC_000913/2286163‑2286457
                                                                                                                                                    |                                                                                                                                                  
tCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTTGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCt                                                                                                                                                   >  1:660510/1‑150 (MQ=255)
  ttACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGAACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGCGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTa                                                                                                                                                 >  2:681750/1‑150 (MQ=255)
   tACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                <  2:1603564/150‑1 (MQ=255)
   tACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                <  1:433145/150‑1 (MQ=255)
    aCTTGCAATCGGGGTCGAAAACGGTAGAATTCGCAGCCCCGAGGCGGCACGTAGCGCATCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAACGCCTCACGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                                                                               <  2:1314976/150‑1 (MQ=255)
       tGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                            >  1:902567/1‑150 (MQ=255)
       tGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGGGCGGACCAATAATCCCAAGAAAAAACCAACCCTTACggtt                                                                                                                                            >  2:272109/1‑150 (MQ=255)
            tCGGTGTGGAAAACTGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGAAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttt                                                                                                                                       <  2:1209986/150‑1 (MQ=255)
             cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTTGTAGCTCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                      <  2:1524750/150‑1 (MQ=255)
             cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttt                                                                                                                                      <  2:813868/150‑1 (MQ=255)
              ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtgttt                                                                                                                                     >  1:1775090/1‑150 (MQ=255)
              ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTAGGttttt                                                                                                                                     >  1:733699/1‑150 (MQ=255)
                tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                                                                   <  1:1316491/150‑1 (MQ=255)
                tgtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                                                                   <  2:1454324/150‑1 (MQ=255)
                tgtgGAAAACGGTAGTATAAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                                                                   <  1:1622676/150‑1 (MQ=255)
                  tgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                                                                 >  1:621200/1‑150 (MQ=255)
                   ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                <  1:870404/150‑1 (MQ=255)
                    gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                                                                               <  1:883337/150‑6 (MQ=255)
                      aaaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCGTACGGTTGGTTTTTTTATATCt                                                                                                                             >  2:1777867/1‑150 (MQ=255)
                       aaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTg                                                                                                                            <  2:693898/150‑1 (MQ=255)
                       aaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGGCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTg                                                                                                                            <  2:287096/150‑1 (MQ=255)
                          cctagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGATTTTTTATATCTGCaa                                                                                                                         >  1:1216711/3‑150 (MQ=255)
                           gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                                                                        <  1:1263434/150‑1 (MQ=255)
                           gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                                                                        <  1:901968/150‑1 (MQ=255)
                           gtagtaTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTATCTGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaat                                                                                                                        >  1:776515/1‑150 (MQ=255)
                            ttttaTTAGCAGCCACGAGTCGGCTCGTAGCGAAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaatt                                                                                                                       <  1:335240/147‑1 (MQ=255)
                               taTTAGCAGCCACGAGTCGGCACGTAGCGCATCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaat                                                                                                                    <  1:42913/150‑1 (MQ=255)
                                 ttAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTc                                                                                                                  <  1:567470/150‑1 (MQ=255)
                                  tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTTTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                 <  1:1271392/150‑1 (MQ=255)
                                  aagcagcCACGAGTCGGCACGTAGCGCATCCTGGTAGCGCAACGTCATGGGTTGTCGGGGGTCGGAGGATCAAATACTCTCGTGACGACCAAAAATCCCAAGAAAAAACCAACCCTAACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                 <  1:149697/149‑1 (MQ=255)
                                   agcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTAGTGCCGACCAAAATTCCCAAGAAAAAACCAACCATTTCGGATGGTTTTTTTATATATGCAATTAATTAGa                                                                                                                >  2:1243681/1‑150 (MQ=255)
                                    gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAt                                                                                                               >  1:1533971/1‑150 (MQ=255)
                                        ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                           <  2:1557844/150‑1 (MQ=255)
                                        ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                           <  1:1734516/150‑1 (MQ=255)
                                        ccACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAc                                                                                                           >  2:1568280/1‑150 (MQ=255)
                                           cGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGa                                                                                                        <  2:1831316/150‑1 (MQ=255)
                                                cGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                                                                   <  2:830771/150‑1 (MQ=255)
                                                cGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTg                                                                                                   <  2:114924/150‑1 (MQ=255)
                                                cGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTGACGGGTGGGTGTTTTATATCTGCAATTAATTCGATAAACCGACCGTg                                                                                                   >  1:1121582/1‑150 (MQ=255)
                                                 ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGa                                                                                                  >  2:864109/1‑150 (MQ=255)
                                                 ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATATGCAATTAATTCGATAAACAGACCGtca                                                                                                  >  2:792893/1‑148 (MQ=255)
                                                 ggCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGAGTCGGGGGTCGGAGGGTGAAGTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCATACGGATGGTTTGTTTATATCTGCAATAAAGTCGATAAACAGAACGTgc                                                                                                  >  1:265077/1‑149 (MQ=255)
                                                    aCGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGacac                                                                                               <  2:1337332/150‑1 (MQ=255)
                                                       tAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAAGCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATc                                                                                            >  1:1318276/1‑150 (MQ=255)
                                                        aGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATca                                                                                           >  2:1695880/1‑150 (MQ=255)
                                                             aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                                                                      >  2:427895/1‑150 (MQ=255)
                                                             aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                                                                      <  1:564031/150‑1 (MQ=255)
                                                             aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGcc                                                                                      <  1:556741/150‑1 (MQ=255)
                                                              gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                                     >  1:706285/1‑150 (MQ=255)
                                                               ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                                                                    <  2:863116/150‑1 (MQ=255)
                                                                   gTAGCGCACCGTCATGGGGTGTCGGGGGTCTGAGGTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAATAACCAACACTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                                                                <  1:1100548/150‑1 (MQ=255)
                                                                    tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                                                               >  2:559639/1‑150 (MQ=255)
                                                                      gcgcACCGTCATGGGGTGGCGGGGGTCGGAGGTTAAAATCATATCGTGACGACCAAAAATCCCAAGAAAAAACCAACCCTTACTGTTGGTTTTTTTATATCTTCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAttt                                                                             <  2:117513/150‑1 (MQ=255)
                                                                       cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                                                            <  1:777708/150‑1 (MQ=255)
                                                                       cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAtttt                                                                            >  1:688536/1‑150 (MQ=255)
                                                                        gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                           <  2:70782/150‑1 (MQ=255)
                                                                        gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                           >  1:409365/1‑150 (MQ=255)
                                                                         cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                                                                          <  1:1051908/150‑1 (MQ=255)
                                                                          aCCGTCATGTGGTGTCGGTGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                                                                         <  2:98485/150‑1 (MQ=255)
                                                                          aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTg                                                                         <  1:657289/150‑1 (MQ=255)
                                                                               cATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATc                                                                    <  2:1332205/150‑1 (MQ=255)
                                                                                  ggggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATACGATAAAAAGACAGTGACACATAACAGCCTGTTTATTTTCTGTTTTCAGa                                                                 >  2:1631468/1‑150 (MQ=255)
                                                                                     gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg                                                              <  2:852035/150‑1 (MQ=255)
                                                                                     gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGGGCCGACCAACACTCCCAAGAAAAAACCAACCCTTACGGGTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGCCCGTGCCACCGCACAGCGTGTTTATTTTCTGTTATCAGAACg                                                              >  2:1842055/1‑150 (MQ=255)
                                                                                     gggtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg                                                              <  2:1501784/148‑1 (MQ=255)
                                                                                        tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATGTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc                                                           <  1:1265433/150‑1 (MQ=255)
                                                                                              gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca                                                     >  2:1567919/1‑150 (MQ=255)
                                                                                                      ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTTTTATCAGAACGTCCAGACCACACCCGcc                                             <  1:1695880/150‑1 (MQ=255)
                                                                                                      ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                             <  2:508235/150‑1 (MQ=255)
                                                                                                      ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGCTTTATTTATATCTGCAATTAATTCGATAAACAGACCGGGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                             >  1:96496/1‑150 (MQ=255)
                                                                                                        cAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTg                                           >  2:837691/1‑150 (MQ=255)
                                                                                                         aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTTGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa                                          <  2:1839706/150‑1 (MQ=255)
                                                                                                         aaaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa                                          >  1:1052711/1‑150 (MQ=255)
                                                                                                          aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                         <  1:828291/150‑1 (MQ=255)
                                                                                                          aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                         <  1:1313504/150‑1 (MQ=255)
                                                                                                            tCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTAAATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTa                                       <  2:349151/150‑1 (MQ=255)
                                                                                                                ctcGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTcc                                   <  1:1553602/150‑1 (MQ=255)
                                                                                                                  cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCAc                                 <  2:1640017/150‑1 (MQ=255)
                                                                                                                        gACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAC‑‑CACCCGCCTGAGTATTCCACGGCGATTc                         >  1:148821/1‑150 (MQ=255)
                                                                                                                           cAAAAATCCCACGAAAAAACCAACCCTTACGGTGGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                        <  2:1216711/150‑1 (MQ=255)
                                                                                                                             aaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAc                      >  1:1190839/1‑150 (MQ=255)
                                                                                                                               caTCCCAAGAAAAAACCACCCCTTAGGGGTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCGGGCCACACCCGCCTGAGTATTCCACGGCGATTCTAccc                    >  2:100823/2‑150 (MQ=255)
                                                                                                                                 tCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGCCCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCt                  >  2:1780095/1‑150 (MQ=255)
                                                                                                                                   ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                >  1:434972/1‑150 (MQ=255)
                                                                                                                                   ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACACGCCTGAGTATTCATCGGCGATTCTACAACTgc                >  1:119380/1‑150 (MQ=255)
                                                                                                                                   ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACAGTGACACATAACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCGCACCCGCCTGAGAATTCCACGGAGATTCTACCCCTgc                >  1:946434/1‑150 (MQ=255)
                                                                                                                                   ccAAGAAAAAACCAACCCTTACGGTTGGTATTTTTATATCTGCAATTAAATCGATACACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGAACACACCCGCCTGAGTATTCCACGCGGATTATACCCCTgc                >  1:71203/1‑150 (MQ=255)
                                                                                                                                    cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg               >  2:443162/1‑150 (MQ=255)
                                                                                                                                      taaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGcc             >  1:1036147/3‑150 (MQ=255)
                                                                                                                                       gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa            >  1:437878/1‑150 (MQ=255)
                                                                                                                                       gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa            >  1:1013462/1‑150 (MQ=255)
                                                                                                                                       gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa            >  1:1108156/1‑150 (MQ=255)
                                                                                                                                       aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATTTCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGACTGTTTATTTTATGTTATCATAACGTCCAGCCCACACCCGCCTGAGTATTCCACGGCGATTCTACCACTGCGCCa            >  1:994181/2‑150 (MQ=255)
                                                                                                                                         aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt          >  2:1788747/1‑150 (MQ=255)
                                                                                                                                         aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCATGCGCCAtt          >  2:98234/1‑150 (MQ=255)
                                                                                                                                          aaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTAATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTg         >  2:1007158/1‑150 (MQ=255)
                                                                                                                                          aaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGTGCCATTg         <  1:1181226/150‑1 (MQ=255)
                                                                                                                                            aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc       <  2:187361/150‑1 (MQ=255)
                                                                                                                                             accaaccCTTACGGTTGGTTTTTTTATATCTGGAATTAATTAGATAAACAGACCGGGACACATGACAGCCTGTTTATTTTCTGTTTTCAGAACGTCCAGCCCACACCCGCCTGAGTATTCCACGGCGAATAGACCCCTGCGGCAtatcca      >  2:915473/1‑145 (MQ=255)
                                                                                                                                                 accCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTTTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGc  <  2:1318276/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGC  >  NC_000913/2286163‑2286457

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