Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,313 T→A 7.4% intergenic (+26/+77) 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,3130TA 7.4% 314.9 / 11.8 117intergenic (+26/+77)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 T (52/56);  new base A (3/6);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

Only 100 of 120 total aligned reads displayed.

CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACT  >  NC_000913/2286164‑2286462
                                                                                                                                                     |                                                                                                                                                     
cTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGCGTCGGCACGTCGCGCAGCCTGGTAGCGCACCGTCATGCTGTTTCGGGCGTCGTATTTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCCACCCtt                                                                                                                                                       >  2:1091388/1‑150 (MQ=255)
  tACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                     <  2:1528928/150‑1 (MQ=255)
  aaCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCGGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                     <  1:1151860/149‑1 (MQ=255)
   aCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                                                                                    <  2:1154210/150‑1 (MQ=255)
       gCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACggttg                                                                                                                                                >  1:1273540/1‑150 (MQ=255)
         aaTCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACggttggt                                                                                                                                              >  1:493796/1‑150 (MQ=255)
            cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                           <  1:1160219/150‑1 (MQ=255)
             ggggtgGAAACCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttt                                                                                                                                          <  2:199165/147‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttt                                                                                                                                          <  1:1109725/150‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                          <  2:1641008/150‑1 (MQ=255)
              gtgtGGAAAACGGTAGTATTAGCAGGCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGTGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         <  2:1488728/150‑1 (MQ=255)
              gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                         >  1:1184411/1‑150 (MQ=255)
                gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                                                                       <  1:1232387/150‑3 (MQ=255)
                gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                                                                       >  1:1506508/1‑150 (MQ=255)
                 tgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                                                                      >  1:1014850/1‑150 (MQ=255)
                  ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTtatata                                                                                                                                     >  2:562072/1‑150 (MQ=255)
                  ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                     <  1:1374692/150‑1 (MQ=255)
                  ggAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTAAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtata                                                                                                                                     >  2:1665872/1‑150 (MQ=255)
                   gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                                                                                    <  2:396264/150‑6 (MQ=255)
                   gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTatat                                                                                                                                    <  2:935948/150‑1 (MQ=255)
                     aaaCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAACCCAACCCTTACGGTTGGTTTTTTTATATCt                                                                                                                                  <  1:983323/150‑1 (MQ=255)
                       aCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGACTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCAGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGc                                                                                                                                <  1:1274680/150‑1 (MQ=255)
                        cGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCa                                                                                                                               <  1:906221/150‑1 (MQ=255)
                               aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaattaatt                                                                                                                        >  1:591612/1‑150 (MQ=255)
                               aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTc                                                                                                                       <  2:193163/150‑1 (MQ=255)
                               aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTATATCTGCAATTAATTc                                                                                                                       <  2:1186363/150‑1 (MQ=255)
                                ttAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGGCGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                      <  2:500527/150‑1 (MQ=255)
                                 tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                      <  2:1642586/150‑1 (MQ=255)
                                  agcagcCACGAGTCGGCACGTAGCGCAGACTGGTAGCTCACCGTCATGGGGTGTCGGGGGTCGGAGGTTAAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGATGGTTTTTTTATATCTGCAATTAATTCGa                                                                                                                     <  1:310355/150‑1 (MQ=255)
                                   gcagcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGaa                                                                                                                    >  1:38940/1‑149 (MQ=255)
                                     agcCACGAGTCGGCACGAAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGGTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTAGATCTGCAATTAATTCGATaa                                                                                                                  >  1:27778/1‑150 (MQ=255)
                                      gcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                 >  2:1352339/1‑150 (MQ=255)
                                      gcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                                 >  2:730096/1‑150 (MQ=255)
                                         aCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCATTTAATTCGATAAACAg                                                                                                              >  1:1430588/1‑150 (MQ=255)
                                            aGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAcc                                                                                                           >  2:1729884/1‑150 (MQ=255)
                                            aGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTGTATATCTGCAATTAATTCGATAAACAGAcc                                                                                                           >  2:1175363/1‑150 (MQ=255)
                                                 gCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAATAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                                                                                      <  1:654118/150‑1 (MQ=255)
                                                    cGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCATTAAATTCGATAAACAGACCGTGAcaca                                                                                                   >  1:717485/1‑150 (MQ=255)
                                                          gcAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGGAATTAATTCGATAAACAGCCCGTGACACATCACAg                                                                                             >  1:1514037/1‑150 (MQ=255)
                                                             gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                                          <  2:173197/150‑1 (MQ=255)
                                                             gCCTGGTAGCGCACCGTCATGGGGTGGCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCt                                                                                          <  2:561981/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTTTCGGGGGTCGGAGTTTAAAATACTCTCGTGCCGAACAAAAATCCCAAGAAAAAACCACCCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                                                                         <  2:310842/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAATAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                                                                         <  2:1282386/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTg                                                                                         >  1:1494035/1‑150 (MQ=255)
                                                               cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                                                                        <  2:1430271/150‑1 (MQ=255)
                                                                tGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGtt                                                                                       <  1:1415329/150‑1 (MQ=255)
                                                                 ggTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGttt                                                                                      <  2:423112/150‑1 (MQ=255)
                                                                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTa                                                                                     <  1:730096/150‑1 (MQ=255)
                                                                   tAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTAt                                                                                    <  2:750524/150‑1 (MQ=255)
                                                                      cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAACCTCTCTTGCCGCCCAAAATCCCAAGAAAAAACCAACCCATACGGTTGGTTTTTTTATATCTTCAATTAATTCGATAAACATACCGTGACACATGCCAGCGGGTTTATTTTc                                                                                >  2:227599/1‑150 (MQ=255)
                                                                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTc                                                                                <  1:1718950/150‑1 (MQ=255)
                                                                        cACCGTCATGGGTTGTCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCCATAAAAAACACACCCTTACGGTTGGTTTTTTTATATCTGCAATTACTTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCt                                                                               <  2:853000/150‑1 (MQ=255)
                                                                        cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGCTTATTTTCt                                                                               <  2:630355/150‑1 (MQ=255)
                                                                          ccGTCATTGGGTGTCGTGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                             <  1:182959/150‑1 (MQ=255)
                                                                          ccGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                             <  1:651457/150‑1 (MQ=255)
                                                                                tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTACCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTATTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAGCACAGCCTGTTTATTTTCTGTTATCAg                                                                       >  2:352987/1‑150 (MQ=255)
                                                                                  gggTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGaa                                                                     <  2:240671/150‑1 (MQ=255)
                                                                                    gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCACAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTTCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACg                                                                   <  1:431416/150‑1 (MQ=255)
                                                                                       tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTcc                                                                <  1:1417421/150‑1 (MQ=255)
                                                                                       tCGGGGGTCGGAGGTTCAAATCCTCTCGTGACGACAAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGTTTTTTTTATATATGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTATATTTTCTGTTTTCAGAACGTcc                                                                >  1:1274539/1‑150 (MQ=255)
                                                                                         gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg                                                              <  2:895874/150‑1 (MQ=255)
                                                                                          ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGa                                                             >  2:680331/1‑150 (MQ=255)
                                                                                             gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca                                                          >  2:118536/1‑150 (MQ=255)
                                                                                             gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTATGTTATCAGAACgaccagacca                                                          >  1:1147443/1‑150 (MQ=255)
                                                                                              tCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac                                                         <  2:472592/150‑1 (MQ=255)
                                                                                                ggaggaTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTAAAGGTAGGTTTTTTTATATCTGCAATTAATTCGAAAAACAGACCGTCCCACAACACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaacc                                                       >  1:1422132/1‑147 (MQ=255)
                                                                                                ggAGGGTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGCCCGTGACACATGACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCacac                                                       >  1:88556/1‑150 (MQ=255)
                                                                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                    <  2:1014850/150‑1 (MQ=255)
                                                                                                     ttCAAATCCTGTCGTGCCGACCAAAAATCCCTAGAAAAAACCAACGCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTAATCAGAACGTCCAGACCACACCCGcc                                                  <  1:1091388/150‑1 (MQ=255)
                                                                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCCGACCACACCCGcc                                                  >  2:1265317/1‑150 (MQ=255)
                                                                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                                  <  2:591612/150‑1 (MQ=255)
                                                                                                         aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCATTTAATTCGATAAACAGACCGTGACCCATCACGGCCTGTGTATTTTGTGTTTTCTGACCGTCCAGCCGGGACCCGGCTGAg                                              >  2:1225230/1‑150 (MQ=255)
                                                                                                         aaTCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAg                                              <  2:1360871/150‑1 (MQ=255)
                                                                                                                 cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCAc                                      <  1:1568113/150‑1 (MQ=255)
                                                                                                                      cGACCAAAAATCCCAAGAAAAAACCAACCCTTTCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAACGTGACACATCACAGCCTGTTTATTTACTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTACACGGCGa                                 >  1:1150715/1‑150 (MQ=255)
                                                                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGCCACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                             >  1:50572/1‑150 (MQ=255)
                                                                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                             >  1:1296966/1‑150 (MQ=255)
                                                                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCCCGGCGATTCtg                            >  2:276970/1‑149 (MQ=255)
                                                                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTa                            >  1:1028760/1‑150 (MQ=255)
                                                                                                                               aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACTGCGATTCTAcccc                        >  2:1674251/1‑150 (MQ=255)
                                                                                                                               aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAcccc                        <  1:586808/150‑1 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                             >  1:410682/1‑94 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGt                                                                             <  2:410682/94‑1 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                     >  2:1526034/1‑150 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                     >  2:569670/1‑150 (MQ=255)
                                                                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                    >  1:23357/1‑150 (MQ=255)
                                                                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                    >  1:1567473/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGAATACCACGGCGATTCTACCCCTGCGCCa                 >  2:92132/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCATGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                 >  1:1018139/1‑150 (MQ=255)
                                                                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa                 >  2:1149262/2‑150 (MQ=255)
                                                                                                                                       aaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGCACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAt                >  1:787609/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt               >  1:792167/1‑150 (MQ=255)
                                                                                                                                           aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc            >  2:1713128/1‑150 (MQ=255)
                                                                                                                                           aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc            >  1:939956/1‑150 (MQ=255)
                                                                                                                                            accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCa           <  1:562072/150‑1 (MQ=255)
                                                                                                                                             ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt          <  2:1273540/150‑1 (MQ=255)
                                                                                                                                              caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATa         <  2:493796/150‑1 (MQ=255)
                                                                                                                                               aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAg        >  2:1426434/1‑150 (MQ=255)
                                                                                                                                                 cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGct      <  2:1271312/150‑1 (MQ=255)
                                                                                                                                                     tACGGTTGGTTTTTTTATATCTGCAATTAATTCTATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGctgact  <  2:824696/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCTGACT  >  NC_000913/2286164‑2286462

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: