Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 17,528 G→C 7.2% G14R (GGA→CGA)  nhaA → Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,5280GC 7.2% 417.6 / 15.9 124G14R (GGA→CGA) nhaANa(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base G (59/55);  new base C (6/3);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 130 total aligned reads displayed.

AATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGAT  >  NC_000913/17380‑17676
                                                                                                                                                    |                                                                                                                                                    
aaTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                     <  1:617189/150‑1 (MQ=255)
aaTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                     <  2:1382073/150‑1 (MQ=255)
aaTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                     <  2:1511708/150‑1 (MQ=255)
aaTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                     <  2:193205/150‑1 (MQ=255)
aaTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                     <  1:1579911/150‑1 (MQ=255)
 aTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGa                                                                                                                                                    >  2:769084/1‑150 (MQ=255)
 aTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGa                                                                                                                                                    >  2:763361/1‑150 (MQ=255)
 aTAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGa                                                                                                                                                    >  2:1100844/1‑150 (MQ=255)
         tAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                            <  2:100402/150‑1 (MQ=255)
         tAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                            >  1:544277/1‑150 (MQ=255)
         tAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCattat                                                                                                                                            >  2:1116065/1‑150 (MQ=255)
            aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTc                                                                                                                                          >  2:749964/1‑149 (MQ=255)
            aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTc                                                                                                                                          <  1:749964/149‑1 (MQ=255)
            aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCt                                                                                                                                         <  1:1483816/150‑1 (MQ=255)
            aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCCCGAGGCatcactgc                                                                                                                                         <  1:841491/150‑7 (MQ=255)
            aaaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTCAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCt                                                                                                                                         <  1:891260/150‑1 (MQ=255)
            aaaaCGAATTATTCATACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCt                                                                                                                                         <  1:739200/150‑1 (MQ=255)
             aaaCGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGAGGCATTATTCtt                                                                                                                                        <  2:42771/150‑1 (MQ=255)
                  aaTTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCa                                                                                                                                    >  2:298056/1‑149 (MQ=255)
                  aaTTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCa                                                                                                                                    >  2:292246/1‑149 (MQ=255)
                  aaTTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCa                                                                                                                                    <  1:298056/149‑1 (MQ=255)
                  aaTTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCa                                                                                                                                    <  1:292246/149‑1 (MQ=255)
                    tttttCCTAAACTATTATGTGATTTTTATGATTTTTGGTGCGTGGTAATATATTTAAAACGATGTATTCTGCTGAAAGAGAAATAAATATTGATACATGTGGATCGATTCTTTAGCAGTGATTCGTGGGGAGGCATTATTCTTATCATTg                                                                                                                                 <  2:1245562/147‑1 (MQ=255)
                        tCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGc                                                                                                                             >  2:1660057/1‑150 (MQ=255)
                        tCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGc                                                                                                                             <  2:424981/150‑1 (MQ=255)
                         ccTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCt                                                                                                                            >  2:1112328/1‑150 (MQ=255)
                           tACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTAt                                                                                                                          >  2:1233032/1‑150 (MQ=255)
                            acacTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATc                                                                                                                         >  2:842666/1‑150 (MQ=255)
                                 ataATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGc                                                                                                                    <  1:1660057/150‑1 (MQ=255)
                                   aaTCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAACAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTg                                                                                                                      <  1:215979/146‑1 (MQ=255)
                                   aaTCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAAAAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTg                                                                                                                      >  2:215979/1‑146 (MQ=255)
                                    aTCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgat                                                                                                                 <  2:410712/150‑1 (MQ=255)
                                     tCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgatg                                                                                                                <  2:474683/150‑1 (MQ=255)
                                         aTTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTa                                                                                                            <  1:1457615/150‑1 (MQ=255)
                                          ttttAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTAt                                                                                                           >  1:703630/1‑150 (MQ=255)
                                           tttAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATg                                                                                                          <  1:1100844/150‑1 (MQ=255)
                                           tttAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATg                                                                                                          >  2:580672/1‑150 (MQ=255)
                                            ttAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAAGATCCGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATgg                                                                                                         <  1:38810/150‑1 (MQ=255)
                                              aaCGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                                                                                                     <  1:360089/120‑1 (MQ=255)
                                              aaCGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                                                                                                     >  2:360089/1‑120 (MQ=255)
                                                 gatgatTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGTTATCCTGGCGATGATTATGGCCAAc                                                                                                    >  2:1515757/1‑150 (MQ=255)
                                                      ttCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcgg                                                                                               <  2:432382/150‑1 (MQ=255)
                                                        cGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAgcggcg                                                                                             <  2:1537585/150‑1 (MQ=255)
                                                         gTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGGGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGgcgc                                                                                            >  1:1067733/1‑150 (MQ=255)
                                                           gCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCaa                                                                                          <  1:657714/150‑1 (MQ=255)
                                                           gCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCaa                                                                                          <  1:842666/150‑1 (MQ=255)
                                                             ggggTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAAcc                                                                                        <  2:852367/150‑1 (MQ=255)
                                                                 tAAAATAGTAAAATCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                    <  2:1135999/150‑1 (MQ=255)
                                                                 tAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACCTCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                    <  2:1130288/150‑1 (MQ=255)
                                                                 tAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCTATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                    <  2:1038040/150‑1 (MQ=255)
                                                                 tAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTg                                                                                    >  2:759250/1‑150 (MQ=255)
                                                                   aaaTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGa                                                                                  >  1:119482/1‑150 (MQ=255)
                                                                       aGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGATATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGt                                                                              >  1:1574753/1‑150 (MQ=255)
                                                                        gTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTa                                                                             >  1:596767/1‑150 (MQ=255)
                                                                          ataaaCTATCTATTCACCAGAAAGTCAAATATAAAGTGAAACATCTGCATCGATGCTTTAGCTCGGATGCCTCGCGAGGTATAATTCTTTTGATTGCCGCTATCGTGGCGATGATTTTGGCCATGAGGGGGGCAACCAGTTGATGGGATc                                                                           <  2:1116411/148‑1 (MQ=255)
                                                                          aaaaaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATc                                                                           >  2:819131/1‑150 (MQ=255)
                                                                             aaCGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACg                                                                        <  1:495185/150‑1 (MQ=255)
                                                                                gATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACt                                                                     <  2:1373227/150‑1 (MQ=255)
                                                                                   cTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCAATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTc                                                                  <  2:1007842/150‑1 (MQ=255)
                                                                                   cTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATGCTTTAGCAGTGATGCCTCGGGCGGCATTATTCTTATCATTGCCGCTATCCTGGCAATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTc                                                                  <  2:1001762/150‑1 (MQ=255)
                                                                                    tATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCa                                                                          >  1:1340686/1‑141 (MQ=255)
                                                                                    tATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCa                                                                          <  2:1340686/141‑1 (MQ=255)
                                                                                    tATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCt                                                                 >  2:1188398/1‑150 (MQ=255)
                                                                                     aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTg                                                                >  1:226322/1‑150 (MQ=255)
                                                                                     aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTg                                                                <  1:769084/150‑1 (MQ=255)
                                                                                     aTTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTg                                                                <  1:763361/150‑1 (MQ=255)
                                                                                      ttCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                               <  1:631833/150‑1 (MQ=255)
                                                                                      ttCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTgg                                                               <  2:540162/150‑1 (MQ=255)
                                                                                        cACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGgag                                                             >  1:1154742/1‑150 (MQ=255)
                                                                                                   gaAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCa                                                  >  2:1055140/1‑150 (MQ=255)
                                                                                                       tAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTc                                              >  2:1669358/1‑150 (MQ=255)
                                                                                                                aaaCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGggttggt                                     <  2:1662947/150‑1 (MQ=255)
                                                                                                                  aCATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTc                                   <  1:1116065/150‑1 (MQ=255)
                                                                                                                      cTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTc                               >  2:20987/1‑150 (MQ=255)
                                                                                                                        gCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGTGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGa                             >  2:441369/1‑150 (MQ=255)
                                                                                                                        gCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGa                             >  1:1187510/1‑150 (MQ=255)
                                                                                                                         cATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGTTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaa                            <  1:1515757/150‑1 (MQ=255)
                                                                                                                          aTCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGaaa                           >  2:4848/1‑150 (MQ=255)
                                                                                                                                 cTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                    >  1:245655/1‑150 (MQ=255)
                                                                                                                                 cTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATcaacaa                    >  1:134193/1‑150 (MQ=255)
                                                                                                                                   ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACaaaa                  >  2:1182225/1‑150 (MQ=255)
                                                                                                                                   ttAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACaaaa                  >  1:703611/1‑150 (MQ=255)
                                                                                                                                     aagaataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc                >  1:510077/8‑150 (MQ=255)
                                                                                                                                     aagaataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc                >  1:410143/8‑150 (MQ=255)
                                                                                                                                     aagaataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc                >  1:1430316/8‑150 (MQ=255)
                                                                                                                                     aGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc                <  2:148375/150‑1 (MQ=255)
                                                                                                                                     aGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAAc                <  1:170293/150‑1 (MQ=255)
                                                                                                                                      gCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGTGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACa               >  2:834153/1‑150 (MQ=255)
                                                                                                                                      gCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGAGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACa               >  2:1603066/1‑150 (MQ=255)
                                                                                                                                       gaataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACAt              >  1:231298/6‑150 (MQ=255)
                                                                                                                                       gaataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACAt              >  2:904160/6‑150 (MQ=255)
                                                                                                                                        aataaTGCCTCCCGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATg             >  1:305186/5‑150 (MQ=255)
                                                                                                                                        aGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATg             >  2:953580/1‑150 (MQ=255)
                                                                                                                                            aTGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGt         <  1:436188/150‑1 (MQ=255)
                                                                                                                                               ccTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTAt      >  1:1469502/1‑150 (MQ=255)
                                                                                                                                               ccTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTAt      >  1:1654179/1‑150 (MQ=255)
                                                                                                                                                cTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATg     <  1:1233032/150‑1 (MQ=255)
                                                                                                                                                 tCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATgg    <  2:544277/150‑1 (MQ=255)
                                                                                                                                                  cGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGa   >  2:1054751/1‑150 (MQ=255)
                                                                                                                                                   gggAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGAt  >  1:1520745/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGAT  >  NC_000913/17380‑17676

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: