Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,314,683 T→A 7.5% intergenic (+25/+35) ompW → / ← yciE outer membrane protein W/putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,6830TA 7.5% 333.4 / 12.6 124intergenic (+25/+35)ompW/yciEouter membrane protein W/putative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (70/43);  new base A (8/2);  total (79/45)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 127 total aligned reads displayed.

CATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGC  >  NC_000913/1314535‑1314831
                                                                                                                                                    |                                                                                                                                                    
cATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAatt                                                                                                                                                     >  1:1034119/1‑149 (MQ=255)
 aTGTCAGTGGGGTCCATGGATATCGATACCACCGCCAATTATAAGCGGGCCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAatat                                                                                                                                                    <  2:2038058/150‑1 (MQ=255)
   gTCAGTGTGGTACTTGGATATCGATACATACGCAATTTATATGCTGGGCGGTGGACAGCAAGACGATAGCGTAATCCTCGATCCTTGGGTGTTTATGTTCTCTGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATAc                                                                                                                                                  <  2:2002286/150‑1 (MQ=255)
   gTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATAc                                                                                                                                                  >  1:1294918/1‑150 (MQ=255)
     cAGTGTGGTACATGGATATCGCTACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGCTGTTTATGTTCTCAGCAGGATATCGATTTTAATTCCGCAAAAAAACGACCCCGTAAAAAACgg                                                                                                                                                >  1:107399/1‑150 (MQ=255)
     cAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACgg                                                                                                                                                >  2:601558/1‑150 (MQ=255)
      aGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTATTATACggg                                                                                                                                               >  2:798209/1‑150 (MQ=255)
      aGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACggg                                                                                                                                               >  2:1249173/1‑150 (MQ=255)
       ttgtgGTACATGGAAATCGATACCACCGCCAATTATAAGCTGGTCGTTGCACAGCAACAAGATAGCGTACGCCTCGATCCGGGGGTGTTATTGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACgggg                                                                                                                                              <  1:1590934/149‑1 (MQ=255)
       gtgtGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACgggg                                                                                                                                              <  1:1095054/150‑1 (MQ=255)
        tgtgGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGt                                                                                                                                             <  1:1202048/150‑1 (MQ=255)
           ggTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCaa                                                                                                                                          >  2:1938327/1‑150 (MQ=255)
           ggTACATGGATATCGAAACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTACGCACAAAAACGACCCAGTAATATACGGGGTCaa                                                                                                                                          >  1:1837233/1‑150 (MQ=255)
            gTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGCCCCCGTAATATACGGGGTCAAt                                                                                                                                         >  2:341689/1‑150 (MQ=255)
            ctACATGGATATCGACACCACCGCCAATTAAAAGCTGGGCGGTGAACAGCAACACCAAAGCGTACGGCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTATATTACGGGGTcgtt                                                                                                                                         <  2:387400/149‑4 (MQ=255)
              aCATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATaa                                                                                                                                       >  1:2102160/1‑150 (MQ=255)
               cATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGTCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAg                                                                                                                                      >  2:1685372/1‑150 (MQ=255)
                aTGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTATATTACGGGGTcgtttttg                                                                                                                                     <  1:1229980/150‑8 (MQ=255)
                     aatCGATACCACCGCAAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCAGTGGGTGGCTATGTTCTCAGCAGGATATCGTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATgg                                                                                                                               <  2:76489/149‑1 (MQ=255)
                      atCGATACCACCGCCAATTATAAGCTGTGCGGTGCACATCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTGTAAGAAGGATATCTTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATgg                                                                                                                               <  2:1739190/150‑1 (MQ=255)
                      atCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATgg                                                                                                                               <  1:695854/150‑1 (MQ=255)
                       tCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGt                                                                                                                              >  2:843429/1‑150 (MQ=255)
                         gATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGAGATCGTTTTTAAATCCGCACAAAAACGACCCCGTAATATACTGGGTCTATAATGACATGGtat                                                                                                                            >  1:2056363/1‑150 (MQ=255)
                          aTACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGtata                                                                                                                           >  2:393616/1‑150 (MQ=255)
                             ccaccGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAg                                                                                                                        >  1:941470/1‑150 (MQ=255)
                             ccaccGCCAATTATAAGATGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTAAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAg                                                                                                                        >  2:315541/1‑150 (MQ=255)
                               accGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTCCGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAggg                                                                                                                      <  2:96871/150‑1 (MQ=255)
                                ccGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGt                                                                                                                     >  1:511457/1‑150 (MQ=255)
                                ccGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTAAATTCCTCACAAAAACGACCCCGTAATATACTCGGTCAATAAGGACATGGTATAAAGGGt                                                                                                                     >  1:1154155/1‑150 (MQ=255)
                                   ccAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCATCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAATGACATGGTATAAAGGGtttt                                                                                                                  >  2:296066/1‑150 (MQ=255)
                                      aTTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTattatt                                                                                                               >  2:616813/1‑150 (MQ=255)
                                           aaGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATttcttg                                                                                                          >  2:395219/1‑149 (MQ=255)
                                               tGGGCGGTGCACAGAAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACATAAACGACCCCGTAATATACGTGGTGAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCtt                                                                                                      <  1:47369/150‑1 (MQ=255)
                                                gggCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGACAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTc                                                                                                     >  2:2079594/1‑150 (MQ=255)
                                                   cGGTGCACAGCAACACGATAGCGTACGCCTCGCTCCGTGGGTGTTTATGTTCGCAGCAGGATATCGTTTTTAATTCCGCACAAAAAAGACCCCGTAATATACGGGGGCAATAAGGACATGGTAAAAAGGGGATTATTTCTTCGATTCTAc                                                                                                  >  1:833030/1‑150 (MQ=255)
                                                       gCACAGCAACACGATAGCGGACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCa                                                                                              <  1:185394/150‑1 (MQ=255)
                                                        cacaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCAt                                                                                             >  2:2102424/1‑150 (MQ=255)
                                                         acaGCAACACGAAAGCGTACGCCTCGATCCGTGGGTGGTTATGTGCTCAGCAGGATATCGTTTTTATTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTCCGCTTCTACGCCATc                                                                                            <  1:125167/150‑1 (MQ=255)
                                                          caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATTCGGGGTCAATAATGACATGGTTTAAAGGGTATTATTTCTTCGCTTCTACGCCATCa                                                                                           >  2:141487/1‑150 (MQ=255)
                                                          caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATcg                                                                                           >  2:236282/1‑149 (MQ=255)
                                                          caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCa                                                                                           <  1:1738012/150‑1 (MQ=255)
                                                          caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAAGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCa                                                                                           >  2:338192/1‑150 (MQ=255)
                                                           aGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAg                                                                                          <  2:1541847/150‑1 (MQ=255)
                                                             cAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGACAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCATCTACGCCATAAGtt                                                                                        >  1:1180372/1‑150 (MQ=255)
                                                                 acGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCag                                                                                    <  1:721721/150‑1 (MQ=255)
                                                                    aTAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGc                                                                                 <  2:189810/150‑1 (MQ=255)
                                                                         gTACGCCTCGATGCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAAtt                                                                            <  2:1092542/150‑1 (MQ=255)
                                                                          tACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTa                                                                           >  2:1503323/1‑150 (MQ=255)
                                                                          tACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTa                                                                           >  2:1155544/1‑150 (MQ=255)
                                                                             gCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATTTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaaa                                                                        >  1:1840968/1‑150 (MQ=255)
                                                                             gCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGGAATATACGGGGTCAATAAGGACATGGTATAAAGGGTAGTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaaa                                                                        >  2:1960436/1‑150 (MQ=255)
                                                                              ccTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaaaa                                                                       >  1:2101880/1‑150 (MQ=255)
                                                                                  gATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                                                                   <  2:68415/150‑1 (MQ=255)
                                                                                  gATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                                                                   <  1:328547/150‑1 (MQ=255)
                                                                                            tgtTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGAGAATTAAAAATTTCACAGTTGttt                                                         >  2:2065983/1‑150 (MQ=255)
                                                                                               ttATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCg                                                      >  1:2019461/1‑150 (MQ=255)
                                                                                               ttATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCg                                                      >  2:1608470/1‑150 (MQ=255)
                                                                                               ttATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCg                                                      >  1:1861124/1‑150 (MQ=255)
                                                                                               ttATGTTCTCAGAAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGGATCAGAGCGAATTAAAAATTTCTCACTTGTTTGCg                                                      <  2:1885667/150‑1 (MQ=255)
                                                                                                tttgttgtCATCAGTATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGTGTCAATTAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTTAAAATTTCTCAGTTGTTTGCgg                                                     <  2:983925/148‑1 (MQ=255)
                                                                                                  tGTTCTCAGCAGGATATCGTTTTTAATTACGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGAGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGATTTTTGCGGaa                                                   >  1:1308799/1‑150 (MQ=255)
                                                                                                   gTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGATCGAATTAAAAATTTCTCAGTTGTTTGCGGAat                                                  >  2:45988/1‑150 (MQ=255)
                                                                                                   gTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGaaa                                                  >  2:798829/1‑149 (MQ=255)
                                                                                                        cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCt                                             <  2:1206490/150‑1 (MQ=255)
                                                                                                        cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCt                                             >  1:1472138/1‑150 (MQ=255)
                                                                                                        cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGAGGTCAATAAGGACATGGTATAAAGTGTATTATTTCTTCGCTTCTAAGCCATCAGTTTCAGAGAGAATTAAAAATTTCTAATTTGTTTGCGAAATATTCt                                             >  2:2050242/1‑150 (MQ=255)
                                                                                                         tgtagGTTATTGTTTTTATTTCCGCAATTAAACGACCCCGTATTATACGGGGTCAATAAGGACATGGTATATAGGGTATTATTTCTTCGCTTCTTCGCCATCAGTTTCAGAGGGATTTAAAAATTTATCAGTTGTTTGCGGAATATTCTg                                            <  2:1180372/147‑1 (MQ=255)
                                                                                                            agGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACTGGGTCAATAATGAAATGGTATAAAGAGTATTATTTCTTCGCTTATACCCCAACAGTTTCAGAGCGAAAAAAAAATTTCTAATTTGTTAGCGGAATATTCTGAAt                                         >  1:112954/1‑150 (MQ=255)
                                                                                                            agGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAAt                                         >  1:332906/1‑150 (MQ=255)
                                                                                                             ggATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTATGAATc                                        >  2:1211615/1‑150 (MQ=255)
                                                                                                               atatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAg                                      <  1:689454/150‑1 (MQ=255)
                                                                                                                 atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                    >  2:532074/1‑150 (MQ=255)
                                                                                                                 atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                    <  2:1619059/150‑1 (MQ=255)
                                                                                                                 atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                    <  1:315541/150‑1 (MQ=255)
                                                                                                                    gTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGAATTATTTCTTCGCTTCAACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGt                                 >  2:544471/1‑150 (MQ=255)
                                                                                                                    gTTTTTAATTCCGCACAAAAACGACCACGTAATATACGGGGTCAATAAGTACATGGTATAAAGGGTATTATTTCTACGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCTGAATATTCTGAATCAGCCatt                                 >  1:94611/1‑148 (MQ=255)
                                                                                                                          aaTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTGAGGTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCa                           >  1:260586/1‑150 (MQ=255)
                                                                                                                          aaTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCa                           <  1:1249173/150‑1 (MQ=255)
                                                                                                                           aTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCAt                          >  1:507027/1‑150 (MQ=255)
                                                                                                                           aTTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCAt                          >  2:967873/1‑150 (MQ=255)
                                                                                                                                gCACAAAAACGACCCCTTAATATACGTGGTAAATAAGGACATGGTATAAAGGGAATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAACTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCt                     <  1:44252/150‑1 (MQ=255)
                                                                                                                                  acaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCttt                   >  1:152928/1‑150 (MQ=255)
                                                                                                                                   ccttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTGGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtttt                  >  1:1814747/8‑150 (MQ=255)
                                                                                                                                   caAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtttt                  <  1:1938327/150‑1 (MQ=255)
                                                                                                                                   caAAAACGACCCCGTAATATACGGGGTCAATAAGGACAAGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTACTCAGTTGTTTGAGGAATATTCTTAATCAGCCAGTCATCCATTTGCtttt                  >  1:391471/1‑150 (MQ=255)
                                                                                                                                    cttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTc                 >  1:866070/7‑150 (MQ=255)
                                                                                                                                     ttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                >  1:63142/6‑150 (MQ=255)
                                                                                                                                     aaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTCCGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                >  2:1912127/1‑150 (MQ=255)
                                                                                                                                     aaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                >  1:1274951/1‑150 (MQ=255)
                                                                                                                                      tattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCt               >  2:497000/5‑150 (MQ=255)
                                                                                                                                         tgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCAt            >  1:1677113/2‑150 (MQ=255)
                                                                                                                                         tgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTATCAGTTGTTTGCGGAATATTCTGAATCAGCCAGACGGCCATTTGCTTTTCCTCAt            >  2:1412778/2‑150 (MQ=255)
                                                                                                                                            ccccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTa         <  2:283175/150‑1 (MQ=255)
                                                                                                                                            ccccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTa         <  2:1067373/150‑1 (MQ=255)
                                                                                                                                             cccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAATAATTTCTCAGTTGGTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaa        <  1:1041323/150‑1 (MQ=255)
                                                                                                                                              ccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaa       <  1:616813/150‑1 (MQ=255)
                                                                                                                                              ccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaa       >  1:744613/1‑150 (MQ=255)
                                                                                                                                               cGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAAATTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCAGTTGCTTTTCCTCATTTaaaa      <  2:1111581/150‑1 (MQ=255)
                                                                                                                                                 tAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATc    <  2:1740338/150‑1 (MQ=255)
                                                                                                                                                   atatACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGc  <  1:541155/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGC  >  NC_000913/1314535‑1314831

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: