Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,695 A→T 5.5% intergenic (+22/+32) ileX → / ← yqjH tRNA‑Ile/putative siderophore interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,215,6950AT 5.5% 391.4 / 9.2 128intergenic (+22/+32)ileX/yqjHtRNA‑Ile/putative siderophore interacting protein
Reads supporting (aligned to +/- strand):  ref base A (59/62);  new base T (2/5);  total (61/67)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 132 total aligned reads displayed.

GCGGTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215548‑3215844
                                                                                                                                                   |                                                                                                                                                     
       tGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAg                                                                                                                                              >  1:1032307/1‑150 (MQ=255)
        ggCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGc                                                                                                                                             <  1:974377/150‑1 (MQ=255)
        ggCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGc                                                                                                                                             >  2:82670/1‑150 (MQ=255)
        ggCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGc                                                                                                                                             >  2:1135987/1‑150 (MQ=255)
         gCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGATTTCTCGTCAGcc                                                                                                                                            <  2:457863/150‑1 (MQ=255)
            ggATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCTGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCCGATATAGCAAAGGCTGACGAGATCTCGGCAGCgttt                                                                                                                                         >  1:1767813/1‑150 (MQ=25)
            ggATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCcttt                                                                                                                                         >  2:1598446/1‑150 (MQ=255)
             gATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttt                                                                                                                                        >  2:416468/1‑150 (MQ=255)
                 gCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctt                                                                                                                                    >  2:1573899/1‑150 (MQ=255)
                 gCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTAAAGTCAAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAAACGTCAGGctttttctt                                                                                                                                    >  1:1196470/1‑150 (MQ=255)
                   gACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttctttt                                                                                                                                  <  1:535252/150‑1 (MQ=255)
                    acacGGAGTTACTTTATACTCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCCAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                                                                                                                                 <  1:216405/150‑1 (MQ=255)
                    acacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCCGTGGGTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCctttttcttttt                                                                                                                                 <  1:1121100/150‑1 (MQ=255)
                    acacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGTAATCGTCAGCctttttcttttt                                                                                                                                 <  2:1638411/150‑1 (MQ=255)
                     cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGTCGACTCAAAATCGCTTGGTCGCTGGTTCAAGTCCAGAAGGGGCCACCAGATATAGCAAAGTCAGACGAGAAATCGTCAGCctttttctttttt                                                                                                                                >  2:1707824/1‑149 (MQ=37)
                     cacGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATGTAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTta                                                                                                                                <  1:1437372/150‑1 (MQ=255)
                       cGGAGTTACTTAATAATCCGCTACCATGGCCCCTTAGCTCAGGGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTtata                                                                                                                              <  1:767978/150‑1 (MQ=255)
                           gTTACTATATAATCCGCTACGATGGCCCCTTAGCTCAGTGGTTAGAGCAGGAGACTCATATTCGCTTGGTCGATGGTTCAAGTGCAGCAGGGGCCACCAGATATAGCATAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATc                                                                                                                          <  2:1367682/150‑1 (MQ=37)
                             tACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAg                                                                                                                        <  2:1474840/150‑1 (MQ=255)
                             tACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAg                                                                                                                        <  2:1334014/150‑1 (MQ=255)
                              aCTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGATTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGCTATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGt                                                                                                                       >  1:1311693/1‑150 (MQ=255)
                               cTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATTTCAGtt                                                                                                                      >  2:45845/1‑150 (MQ=255)
                               cTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGtt                                                                                                                      >  1:720308/1‑150 (MQ=255)
                                 ttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCCCCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTAc                                                                                                                    >  2:1313438/1‑150 (MQ=255)
                                 ttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTAc                                                                                                                    >  2:953741/1‑150 (MQ=255)
                                 ttATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTAc                                                                                                                    <  1:1563846/150‑1 (MQ=255)
                                   ataATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACtt                                                                                                                  >  1:1041953/1‑150 (MQ=255)
                                     aaTCCGCTACCATGGCCCCTTAGCTCAGGGGTTATAGTAGGCGACTAATAATCGCTTGGTCGCTGGTTCAATTACAGCATGTGCCACCATATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTg                                                                                                                <  1:1619422/150‑1 (MQ=37)
                                     aaTACGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGCCTGACGAGAAATCGTCACCCTTATTCATATTATATATCAGTTACTTTg                                                                                                                >  2:880581/1‑150 (MQ=25)
                                      aTCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGc                                                                                                               >  1:208647/1‑150 (MQ=255)
                                      aTCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGc                                                                                                               >  2:1099980/1‑150 (MQ=255)
                                       tCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAAACGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTATTATATATAAGTTACTTTGCg                                                                                                              >  1:1659359/1‑150 (MQ=255)
                                        ccGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGt                                                                                                             <  2:901220/150‑1 (MQ=255)
                                           cTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGcc                                                                                                          <  1:665419/150‑1 (MQ=255)
                                            tACCATGGCCCCTTAGCTCAGTGGTTAGCGCAGGCGACTCATAATCGCTTGGGCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCa                                                                                                         <  2:199874/150‑1 (MQ=255)
                                            tACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCa                                                                                                         >  2:626624/1‑150 (MQ=255)
                                            tACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCa                                                                                                         <  2:226253/150‑1 (MQ=255)
                                             aCCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAg                                                                                                        >  1:1049445/1‑150 (MQ=255)
                                                aTGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTaa                                                                                                     >  1:629887/1‑150 (MQ=255)
                                                   gCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGcc                                                                                                  <  2:1621401/150‑1 (MQ=255)
                                                   gCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGcc                                                                                                  <  2:1690278/150‑1 (MQ=255)
                                                   gCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGcc                                                                                                  >  2:544946/1‑150 (MQ=255)
                                                    ccccTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCg                                                                                                 >  1:1054608/1‑150 (MQ=255)
                                                        ttAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGc                                                                                             <  2:1096701/150‑1 (MQ=255)
                                                         tAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCa                                                                                            >  2:602873/1‑150 (MQ=255)
                                                           gCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACg                                                                                          >  2:1095693/1‑150 (MQ=255)
                                                              cAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAc                                                                                       >  1:72948/1‑150 (MQ=255)
                                                              cAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAc                                                                                       >  2:710539/1‑150 (MQ=255)
                                                                 tGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGCTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCg                                                                                    <  1:552276/150‑1 (MQ=255)
                                                                  ggTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGc                                                                                   <  1:1489810/150‑1 (MQ=255)
                                                                        cgCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                             <  1:839754/149‑1 (MQ=255)
                                                                        agCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                             <  1:532314/150‑1 (MQ=255)
                                                                        agCAGGCGACTCATAATCGCTTGGTCGCTGGTTAAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTggg                                                                             <  2:1325480/150‑1 (MQ=255)
                                                                          cAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTc                                                                           <  2:1577116/150‑1 (MQ=255)
                                                                            ggCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCAt                                                                         >  1:1743090/1‑150 (MQ=255)
                                                                                 cTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGAAGAGGCCTCCAGAAATAGCAAAAGCTGACGAGAAAACGTCAGCCTTTTTCTTGTTATATATCAGATTCTTAGCGTCCCAGTAAGCCGCTTCACGTACCCCCTGTGGTTCATACTGt                                                                    >  1:274501/1‑150 (MQ=255)
                                                                                  tCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAACTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTAACCGCTGTGGGTCATACTGtt                                                                   >  1:1675179/1‑150 (MQ=255)
                                                                                   cATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTc                                                                  <  1:248442/150‑1 (MQ=255)
                                                                                    aTAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTcc                                                                 <  2:20168/150‑1 (MQ=255)
                                                                                         cGCTTGGTCGGTGGGTCAAGTCCAGCAGGGGCCACCAGATTTTGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTc                                                            <  2:476062/150‑1 (MQ=255)
                                                                                          gCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGCCGAGAAGTCGTCAGCCTGGTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCGGTGGGTCATACTGTTCCGCttaa                                                           >  2:299127/1‑148 (MQ=255)
                                                                                            ttGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCaaa                                                         <  2:668768/150‑1 (MQ=255)
                                                                                             tGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAg                                                        <  1:1617274/150‑1 (MQ=255)
                                                                                              ggTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgc                                                       <  2:117778/150‑1 (MQ=255)
                                                                                               gTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcg                                                      <  2:1656936/150‑1 (MQ=255)
                                                                                               gTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcg                                                      <  1:1336289/150‑1 (MQ=255)
                                                                                                 cGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGAAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCAGCTTCAAAgcggc                                                    >  2:908159/1‑150 (MQ=255)
                                                                                                  gCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAgcggcg                                                   <  1:626624/150‑1 (MQ=255)
                                                                                                   cTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGcggcgg                                                  >  2:1251517/1‑150 (MQ=255)
                                                                                                     ggTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCAGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCt                                                >  2:268305/1‑150 (MQ=255)
                                                                                                       ttCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCGTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTa                                              >  2:800420/1‑150 (MQ=255)
                                                                                                        tCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTaa                                             <  2:741385/150‑1 (MQ=255)
                                                                                                         cAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTaaa                                            <  2:1715473/150‑1 (MQ=255)
                                                                                                             tCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCt                                        <  2:1025362/150‑1 (MQ=255)
                                                                                                               cagcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTa                                      <  1:1573899/150‑1 (MQ=255)
                                                                                                               cagcagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTGTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTCCACGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTcgtc                                      >  2:1165194/1‑147 (MQ=255)
                                                                                                                agcagGGGCCACCATATATAGCAAAGGCTGTCGTGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGACAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTaa                                     <  2:787772/150‑1 (MQ=255)
                                                                                                                  cagGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACg                                   <  2:649166/150‑1 (MQ=255)
                                                                                                                   agGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGa                                  <  1:167540/150‑1 (MQ=255)
                                                                                                                    ggggCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGAc                                 >  2:1520424/1‑150 (MQ=255)
                                                                                                                     gggCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACt                                >  2:637751/1‑150 (MQ=255)
                                                                                                                     gggCCACCAGATATACCAAAGTCTGACGAGAAATCTTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACt                                <  2:1512430/150‑1 (MQ=255)
                                                                                                                       gCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACttt                              <  2:208647/150‑1 (MQ=255)
                                                                                                                       gCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACttt                              <  2:787121/150‑1 (MQ=255)
                                                                                                                        ccaccaGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTa                             <  1:1625631/150‑1 (MQ=255)
                                                                                                                                 ataGCTAAGGCTGTCGAGAAATCGTCAGCCTTTTTCTTTTAATATATCAGTTTCTTTGCGTGCCAGTTAGCCGCTGGACGTAGCCGCTGTGGGTGATACTTTTGCGCTTATAAGCTGCGGCTTAAATTCTTAACGACTTTACCTTCTcct                    <  2:791342/150‑2 (MQ=255)
                                                                                                                                  taGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTAACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCgg                   <  1:82670/150‑1 (MQ=255)
                                                                                                                                   aGCAAAGGCTGACGTGAAATCGTGAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCTCGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  1:64506/150‑1 (MQ=255)
                                                                                                                                   aGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGt                  <  2:1228762/150‑1 (MQ=255)
                                                                                                                                    gCAAAGGCTGACGAGAATTCGTCAGCCTTTTTATTTTTATATTTCAGTTACTTTGCGTGCCAGTAATCCGCTGCACATACCAGATGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAAGTCTTTACTACTTTACCTTCGCAGGtt                 >  1:125998/1‑150 (MQ=255)
                                                                                                                                    gCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGtt                 >  1:141367/1‑150 (MQ=255)
                                                                                                                                     cAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGGTa                >  1:1310948/1‑150 (MQ=255)
                                                                                                                                      aaaGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTAt               >  2:1400898/1‑150 (MQ=255)
                                                                                                                                        aGGCTGACGAGAAATCGTCAGGCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCATTAAGCCGCTGCACGTTCCCGCTGTGTGTCATACTGTTCCGCTTCATATCGGCGGCTTAAATTGTTTACGACTTTACCTTCGCCGGTTATcc             <  1:880581/150‑1 (MQ=255)
                                                                                                                                         ggCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCa            <  1:672520/150‑1 (MQ=255)
                                                                                                                                         ggCTGACGATTTCTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCa            >  1:1771691/1‑150 (MQ=255)
                                                                                                                                          gCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAg           <  2:1041953/150‑1 (MQ=255)
                                                                                                                                             gACGAGAAATCGTCTGCATTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCTTCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATg        <  1:1366731/150‑1 (MQ=255)
                                                                                                                                                 agaAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTATATTCTTTACGACTTTACCTTCGCCGGTTATCCAGATGAAGt    >  2:714403/1‑150 (MQ=255)
                                                                                                                                                   aaaTCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTaa  <  1:687102/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGGTAATGGCTGGATTGCGACACGGAGTTACTTTATAATCCGCTACCATGGCCCCTTAGCTCAGTGGTTAGAGCAGGCGACTCATAATCGCTTGGTCGCTGGTTCAAGTCCAGCAGGGGCCACCAGATATAGCAAAGGCTGACGAGAAATCGTCAGCCTTTTTCTTTTTATATATCAGTTACTTTGCGTGCCAGTAAGCCGCTGCACGTACCCGCTGTGGGTCATACTGTTCCGCTTCAAAGCGGCGGCTTAAATTCTTAACGACTTTACCTTCGCCGGTTATCCAGATGAAGTAA  >  NC_000913/3215548‑3215844

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: