Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,854,024 C→A 6.3% intergenic (‑45/+72) ydjE ← / ← ydjF putative transporter YdjE/putative DNA‑binding transcriptional regulator YdjF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,854,0240CA 6.3% 389.0 / 10.0 117intergenic (‑45/+72)ydjE/ydjFputative transporter YdjE/putative DNA‑binding transcriptional regulator YdjF
Reads supporting (aligned to +/- strand):  ref base C (54/54);  new base A (5/3);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 120 total aligned reads displayed.

CCCTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGTA  >  NC_000913/1853875‑1854173
                                                                                                                                                     |                                                                                                                                                     
cccTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTc                                                                                                                                                       >  1:215650/1‑150 (MQ=255)
   tGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACt                                                                                                                                                    >  2:1295928/1‑150 (MQ=255)
    gTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACtt                                                                                                                                                   >  2:408961/1‑150 (MQ=255)
      taataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGc                                                                                                                                                 >  2:45867/1‑150 (MQ=255)
      taataaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGc                                                                                                                                                 >  1:1303360/1‑150 (MQ=255)
         taaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGc                                                                                                                                              >  2:1424916/1‑150 (MQ=255)
         taaCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGc                                                                                                                                              <  1:45867/150‑1 (MQ=255)
           aCAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaa                                                                                                                                            <  2:1132795/150‑2 (MQ=255)
            cAGACTACAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCCTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAGTGAATTTGCCGCaat                                                                                                                                           <  1:754424/150‑3 (MQ=255)
              gACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCtt                                                                                                                                         <  2:1246850/150‑1 (MQ=255)
              gACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCtt                                                                                                                                         >  2:573748/1‑150 (MQ=255)
                cTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                       <  1:1555645/150‑1 (MQ=255)
                       ttATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATcatcat                                                                                                                                >  2:819199/1‑150 (MQ=255)
                        tataATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCAtt                                                                                                                               <  1:793778/150‑1 (MQ=255)
                            ataCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAt                                                                                                                           >  2:1078029/1‑150 (MQ=255)
                             taCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                       >  2:1075950/1‑137 (MQ=255)
                             taCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTa                                                                                                                                       <  1:1075950/137‑1 (MQ=255)
                              aCCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttt                                                                                                                         >  1:599640/1‑150 (MQ=255)
                               ccAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTAttta                                                                                                                        >  1:53598/1‑149 (MQ=255)
                               ccAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTc                                                                                                                        <  2:901193/150‑1 (MQ=255)
                                aaaaTATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTTATCATATTGTTCCATTAGGTGCCTATGATTTGCTAAAAAGGTAGCAATTTCGTCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCa                                                                                                                       <  2:666337/149‑1 (MQ=255)
                                    tataCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGgc                                                                                                                   <  1:1485243/150‑1 (MQ=255)
                                       aCGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCa                                                                                                                >  2:225203/1‑150 (MQ=255)
                                        cGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAg                                                                                                               >  1:471470/1‑150 (MQ=255)
                                        cGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAg                                                                                                               >  1:1001112/1‑150 (MQ=255)
                                          aTAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGa                                                                                                             >  1:1396953/1‑150 (MQ=255)
                                          aTAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGa                                                                                                             <  2:150863/150‑1 (MQ=255)
                                            aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAAc                                                                                                           <  1:708274/150‑1 (MQ=255)
                                            aaTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAAc                                                                                                           <  2:215650/150‑1 (MQ=255)
                                             aTGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACa                                                                                                          >  1:1639351/1‑150 (MQ=255)
                                              tGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACaa                                                                                                         >  1:1312501/1‑150 (MQ=255)
                                              tGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTACCATCATTTTATTTCAAGGCGCAGGAACaa                                                                                                         >  2:1351998/1‑150 (MQ=255)
                                                aaaCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaat                                                                                                       <  2:801474/150‑1 (MQ=255)
                                                 aaCCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaata                                                                                                      >  1:250528/1‑150 (MQ=255)
                                                   ccGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAaataat                                                                                                    <  2:741289/150‑1 (MQ=255)
                                                        ccAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAgg                                                                                                                    <  2:421513/129‑1 (MQ=255)
                                                             ggCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACg                                                                                          >  1:617245/1‑150 (MQ=255)
                                                               cAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACgtg                                                                                        <  2:275180/150‑1 (MQ=255)
                                                                 aGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                      >  2:1551386/1‑150 (MQ=255)
                                                                 aGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGtgtg                                                                                      >  1:1247718/1‑150 (MQ=255)
                                                                     gTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTa                                                                                  <  2:467646/150‑1 (MQ=255)
                                                                          gTCTTGCGCCAATTTGATCATATTGTTCCATTTGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTa                                                                             <  1:40868/150‑1 (MQ=255)
                                                                          gTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTa                                                                             <  1:1152294/150‑1 (MQ=255)
                                                                            cTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTAcc                                                                           <  1:1424916/150‑1 (MQ=255)
                                                                              tGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGATTAAATTACcaa                                                                         >  2:455884/1‑149 (MQ=255)
                                                                                 gcCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACAcc                                                                      <  2:608994/150‑1 (MQ=255)
                                                                                   cAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                    <  2:113032/150‑1 (MQ=255)
                                                                                       ttGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTgg                                                                <  2:468101/150‑1 (MQ=255)
                                                                                       ttGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTgg                                                                >  2:946185/1‑150 (MQ=255)
                                                                                        tGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGa                                                               >  2:996984/1‑150 (MQ=255)
                                                                                          aTCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGata                                                             <  2:471470/150‑1 (MQ=255)
                                                                                           tCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGatat                                                            >  2:379972/1‑150 (MQ=255)
                                                                                             atatTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGAtatta                                                          <  2:417357/150‑1 (MQ=255)
                                                                                             atatTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGAtatta                                                          >  2:601877/1‑150 (MQ=255)
                                                                                              tatTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGAtattat                                                         >  2:764883/1‑150 (MQ=255)
                                                                                               atTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTAACTTGCAGCATCTTTATCATCATTTTATTTCAAGGCTCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCCAACCAGCTGCATaatatt                                                        >  2:444366/1‑150 (MQ=255)
                                                                                                  gTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTc                                                     <  1:1551386/150‑1 (MQ=255)
                                                                                                   ttCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCt                                                    >  1:1560027/1‑150 (MQ=255)
                                                                                                   ttCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCt                                                    <  2:599640/150‑1 (MQ=255)
                                                                                                     ccATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCttt                                                  >  2:565013/1‑150 (MQ=255)
                                                                                                       aTTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGc                                                <  1:127567/150‑1 (MQ=255)
                                                                                                       aTTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGc                                                <  2:620363/150‑1 (MQ=255)
                                                                                                           ggTGCCTCGTATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAAGGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCaaaa                                            <  2:53598/150‑1 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGTTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGc                                         <  2:1101237/150‑1 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGc                                                <  1:585386/143‑1 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAAtac                                         >  1:331982/1‑148 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGc                                         >  2:563918/1‑150 (MQ=255)
                                                                                                              gCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGc                                         >  1:1317977/1‑150 (MQ=255)
                                                                                                                  cGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATc                                     <  2:897988/150‑1 (MQ=255)
                                                                                                                   ggATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATcc                                    <  1:225203/150‑1 (MQ=255)
                                                                                                                     aTTTGCTAAAAAGGTAGCAATTGCGGAAAATTCACTTGCCGCATCTTTATCATAATTTTATTTAAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCAc                                  >  2:1648913/1‑150 (MQ=255)
                                                                                                                       ttGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCActc                                <  1:601877/150‑1 (MQ=255)
                                                                                                                          cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                             >  1:554425/1‑150 (MQ=255)
                                                                                                                          cTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGc                             <  1:1589380/150‑1 (MQ=255)
                                                                                                                           tAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCa                            >  2:554251/1‑150 (MQ=255)
                                                                                                                             aaaaGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCAcc                          <  2:79680/150‑1 (MQ=255)
                                                                                                                                 ggTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGt                      <  1:1060856/150‑1 (MQ=255)
                                                                                                                                    taaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   >  2:1025480/7‑150 (MQ=255)
                                                                                                                                    aGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   <  1:342548/150‑1 (MQ=255)
                                                                                                                                     aaagatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAtt                  >  1:840968/6‑150 (MQ=255)
                                                                                                                                      cAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAttt                 <  2:1396953/150‑1 (MQ=255)
                                                                                                                                      cAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAttt                 <  1:764883/150‑1 (MQ=255)
                                                                                                                                        gatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   <  2:1082191/144‑1 (MQ=255)
                                                                                                                                        gatGCGGCAAGTGAATTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAt                   >  1:1082191/3‑146 (MQ=255)
                                                                                                                                        aTTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTAt               <  2:1639351/150‑1 (MQ=255)
                                                                                                                                          tGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCa             <  2:1303360/150‑1 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                    <  1:364693/94‑1 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAg                                                                    >  2:364693/1‑94 (MQ=255)
                                                                                                                                           gCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAAAAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAg            <  2:905651/150‑1 (MQ=255)
                                                                                                                                            cGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGt           >  1:1479484/1‑150 (MQ=255)
                                                                                                                                             ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt          >  1:191041/1‑150 (MQ=255)
                                                                                                                                             ggCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtt          >  1:534776/1‑150 (MQ=255)
                                                                                                                                              gCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtta         >  1:187945/1‑150 (MQ=255)
                                                                                                                                              gCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGtta         >  2:430135/1‑150 (MQ=255)
                                                                                                                                               cAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGAttt                 <  2:1219461/141‑1 (MQ=255)
                                                                                                                                                 aaTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatta      <  1:498266/150‑1 (MQ=255)
                                                                                                                                                  aTTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGttatta      <  1:1411997/149‑1 (MQ=255)
                                                                                                                                                   ttCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATg    >  2:1143591/1‑150 (MQ=255)
                                                                                                                                                     cACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGta  <  2:1392807/150‑1 (MQ=255)
                                                                                                                                                     cACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGta  >  1:38684/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCTGTTAATAACAGACTAAAGCTTATAATACCAAATATACGATAATGAAACCGGGCCAAAGGCAAGCGGTCCAGTCTTGCGCCAATTTGATCATATTGTTCCATTAGGTGCCTCGGATTTGCTAAAAAGGTAGCAATTGCGGCAAATTCACTTGCCGCATCTTTATCATCATTTTATTTCAAGGCGCAGGAACAAATAATATTGACAACGTGTGCTTAAATTACCACACCAGCTGGATATTATTGTCTTTGCAAAATGCTATCCACTCTGCACCCGGTGATTTATCAGTTATTATGTA  >  NC_000913/1853875‑1854173

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