Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,311 C→G 6.8% 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 6.8% 346.1 / 14.0 117intergenic (+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 (54/54);  new base G (2/6);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 118 total aligned reads displayed.

CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCT  >  NC_000913/2286164‑2286458
                                                                                                                                                   |                                                                                                                                                    
cTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGCGTCGGCACGTCGCGCAGCCTGGTAGCGCACCGTCATGCTGTTTCGGGCGTCGTATTTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCCACCCtt                                                                                                                                                    >  2:1091388/1‑150 (MQ=255)
  aaCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCGGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTAc                                                                                                                                                  <  1:1151860/149‑1 (MQ=255)
   aCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACg                                                                                                                                                 <  2:1154210/150‑1 (MQ=255)
           tCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttt                                                                                                                                         <  2:486656/150‑1 (MQ=255)
            cGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGtttt                                                                                                                                        <  1:1160219/150‑1 (MQ=255)
             ggggtgGAAACCGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttt                                                                                                                                       <  2:199165/147‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAGGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGTTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAATCCAACCCTTACGGTTGGttttt                                                                                                                                       <  2:1305722/150‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttt                                                                                                                                       <  1:1109725/150‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                       <  2:1641008/150‑1 (MQ=255)
             ggTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATACTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttt                                                                                                                                       >  1:1360871/1‑150 (MQ=255)
              gtgtGGAAAACGGTAGTATTAGCAGGCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGTGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                      <  2:1488728/150‑1 (MQ=255)
              gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                      >  1:1184411/1‑150 (MQ=255)
              gtgtGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAACATCCCAAGAAAAAACCAACCCTTACGGTTGGtttttt                                                                                                                                      >  1:554240/1‑150 (MQ=255)
               ggtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCAAGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGttttttt                                                                                                                                     <  2:787217/149‑1 (MQ=255)
                gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttct                                                                                                                                    <  1:1232387/150‑3 (MQ=255)
                gtgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTta                                                                                                                                    >  1:1506508/1‑150 (MQ=255)
                 tgGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTtat                                                                                                                                   >  1:1014850/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)
                   gAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCGTAAGGGTTGGttttttcttgg                                                                                                                                 <  1:370733/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)
                               aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTTTTTATATCTGCAATTAATTc                                                                                                                    <  2:1186363/150‑1 (MQ=255)
                               aTTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGG‑TTTTTTATATCTGCAATTAATTc                                                                                                                    <  2:193163/150‑1 (MQ=255)
                                ttAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGGC‑GGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                   <  2:500527/150‑1 (MQ=255)
                                 tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                   <  1:1718120/150‑1 (MQ=255)
                                 tAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCg                                                                                                                   <  2:1642586/150‑1 (MQ=255)
                                  agcagcCACGAGTCGGCACGTAGCGCCGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGGTGGGTTTTTTATATCTGCAATTAATTCGa                                                                                                                  >  1:122292/1‑150 (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:730096/1‑150 (MQ=255)
                                      gcCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATaaa                                                                                                              >  2:1352339/1‑150 (MQ=255)
                                        cACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                                                                            >  1:1471314/1‑150 (MQ=255)
                                         aCGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGCCCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCATTTAATTCGATAAACAg                                                                                                           >  1:1430588/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)
                                                            aGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAAAAGACCGTGACACATC‑ACAGcc                                                                                       >  2:586808/1‑150 (MQ=255)
                                                             gCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCt                                                                                      <  2:173197/150‑1 (MQ=255)
                                                             gCCTGGTAGCGCACCGTCATGGGGTGGCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCt                                                                                      <  2:561981/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTTTCGGGGGTCGGAGTTTAAAATACTCTCGTGCCGAACAAAAATCCCAAGAAAAAACCACCCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTg                                                                                     <  2:310842/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAATAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTg                                                                                     <  2:1282386/150‑1 (MQ=255)
                                                              ccTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTg                                                                                     >  1:1494035/1‑150 (MQ=255)
                                                               cTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGt                                                                                    <  2:1430271/150‑1 (MQ=255)
                                                                  gTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTa                                                                                 <  1:730096/150‑1 (MQ=255)
                                                                      cgcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAAACCTCTCTTGCCG‑CCCAAAATCCCAAGAAAAAACCAACCCATACGGTTGGTTTTTTTATATCTTCAATTAATTCGATAAACATACCGTGACACATG‑CCAGCGGGTTTATTTTc                                                                            >  2:227599/1‑150 (MQ=255)
                                                                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTc                                                                            <  1:1718950/150‑1 (MQ=255)
                                                                       gcACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTACCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTc                                                                            <  2:343643/150‑1 (MQ=255)
                                                                        cACCGTCATGGGTTGTCGGGGGGCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCCATAAAAAACACACCCTTACGGTTGGTTTTTTTATATCTGCAATTACTTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCt                                                                           <  2:853000/150‑1 (MQ=255)
                                                                        cACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGCTTATTTTCt                                                                           <  2:630355/150‑1 (MQ=255)
                                                                         aCCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACAATTACGGTTTGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTg                                                                          <  2:197207/150‑1 (MQ=255)
                                                                          ccGTCATTGGGTGTCGTGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGt                                                                         <  1:182959/150‑1 (MQ=255)
                                                                                tGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTACCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTATTTATATCTGCAATTAATTCGATAAACAGACCGTGACACAGC‑ACAGCCTGTTTATTTTCTGTTATCAg                                                                   >  2:352987/1‑150 (MQ=255)
                                                                                    gtgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCACAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTTCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACg                                                               <  1:431416/150‑1 (MQ=255)
                                                                                     tgtCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGt                                                              <  2:1506508/150‑1 (MQ=255)
                                                                                       tCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTcc                                                            <  1:1417421/150‑1 (MQ=255)
                                                                                       tCGGGGGTCGGAGGTTCAAATCCTCTCGTGACGACAAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGTTTTTTTTATATATGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTATATTTTCTGTTTTCAGAACGTcc                                                            >  1:1274539/1‑150 (MQ=255)
                                                                                         gggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAg                                                          <  2:895874/150‑1 (MQ=255)
                                                                                          ggggTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGa                                                         >  2:680331/1‑150 (MQ=255)
                                                                                             gTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACca                                                      >  2:118536/1‑150 (MQ=255)
                                                                                              tCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcac                                                     <  2:472592/150‑1 (MQ=255)
                                                                                                ggaggaTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTAAAGGTAGGTTTTTTTATATCTGCAATTAATTCGAAAAACAGACCGTCCCACAAC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACcaacc                                                   >  1:1422132/1‑147 (MQ=255)
                                                                                                ggAGGGTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGCCCGTGACACATG‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCacac                                                   >  1:88556/1‑150 (MQ=255)
                                                                                                 gAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACAcc                                                  <  2:1669712/150‑1 (MQ=255)
                                                                                                   ggTTCAAATCCTCTCGTGCCGACCCAAACTCCCAAGAAAAAACCAACCCTTACGGTTGGTGTGTTTAGGGCTGCAATTTATTCGATAAACAGACCGTGACCCATCGACTG‑CTGTTTATTTTGTGTTATCAGAACGTCCAGCCAACACCCg                                                >  2:1194720/1‑150 (MQ=255)
                                                                                                   ggTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCg                                                <  2:1014850/150‑1 (MQ=255)
                                                                                                     ttCAAATCCTGTCGTGCCGACCAAAAATCCCTAGAAAAAACCAACGCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTAATCAGAACGTCCAGACCACACCCGcc                                              <  1:1091388/150‑1 (MQ=255)
                                                                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCCGACCACACCCGcc                                              >  2:1265317/1‑150 (MQ=255)
                                                                                                     ttCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGcc                                              <  2:591612/150‑1 (MQ=255)
                                                                                                                 cGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCAc                                  <  1:1568113/150‑1 (MQ=255)
                                                                                                                      cGACCAAAAATCCCAAGAAAAAACCAACCCTTTCGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGAACGTGACACATC‑ACAGCCTGTTTATTTACTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTACACGGCGa                             >  1:1150715/1‑150 (MQ=255)
                                                                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGCCACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                         >  1:50572/1‑150 (MQ=255)
                                                                                                                          cAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCt                         >  1:1296966/1‑150 (MQ=255)
                                                                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCCCGGCGATTCtg                        >  2:276970/1‑149 (MQ=255)
                                                                                                                           aaaaaTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTa                        >  1:1028760/1‑150 (MQ=255)
                                                                                                                               aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACTGCGATTCTAcccc                    >  2:1674251/1‑150 (MQ=255)
                                                                                                                               aTCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTAcccc                    <  1:586808/150‑1 (MQ=255)
                                                                                                                                 cccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTg                  >  1:1114654/1‑150 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGt                                                                         <  2:410682/94‑1 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGt                                                                         >  1:410682/1‑94 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                 >  2:1526034/1‑150 (MQ=255)
                                                                                                                                  ccAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgc                 >  2:569670/1‑150 (MQ=255)
                                                                                                                                   cAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTgcg                >  1:1567473/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGAATACCACGGCGATTCTACCCCTGCGCCa             >  2:92132/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCATGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa             >  1:1018139/1‑150 (MQ=255)
                                                                                                                                      aaaaaaaCCAACCGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCa             >  2:1149262/2‑150 (MQ=255)
                                                                                                                                       aaaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGCACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAt            >  1:787609/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCAtt           >  1:792167/1‑150 (MQ=255)
                                                                                                                                           aaccaaccGTAAGGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc        >  2:1713128/1‑150 (MQ=255)
                                                                                                                                           aaccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGcc        >  1:939956/1‑150 (MQ=255)
                                                                                                                                            accaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCa       <  1:562072/150‑1 (MQ=255)
                                                                                                                                             ccaaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCAt      <  2:1273540/150‑1 (MQ=255)
                                                                                                                                              caaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATa     <  2:493796/150‑1 (MQ=255)
                                                                                                                                               aaccCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAg    >  2:1426434/1‑150 (MQ=255)
                                                                                                                                                 cccTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGct  <  2:1271312/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATC‑ACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGAGTATTCCACGGCGATTCTACCCCTGCGCCATTGCCATAGCT  >  NC_000913/2286164‑2286458

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: ‑