Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,309,870 C→A 5.4% G371C (GGC→TGC)  alsA ← D‑allose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,309,8700CA 5.4% 290.6 / 6.5 109G371C (GGC→TGC) alsAD‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (61/41);  new base A (3/4);  total (64/45)
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.88e-01

Only 100 of 114 total aligned reads displayed.

GATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAC  >  NC_000913/4309721‑4310019
                                                                                                                                                     |                                                                                                                                                     
gaTCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                                                                       >  1:1476688/1‑150 (MQ=255)
 aTCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGccgccg                                                                                                                                                      <  1:1244317/150‑1 (MQ=255)
   cAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTc                                                                                                                                                    <  1:48891/150‑1 (MQ=255)
   cAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTc                                                                                                                                                    >  1:1223751/1‑150 (MQ=255)
   cAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTc                                                                                                                                                    <  2:1337625/150‑1 (MQ=255)
    aGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATAATCTGGTTTTCCGAATGACATTTCAGCGCCAGCAGATCGCGTTGATTTTCAGCGGTACGTTGGTCGTCTACTACATGAATAAAGCCCATCGCGCCTTTAAAGCCGCCGTCt                                                                                                                                                   >  2:88930/1‑150 (MQ=255)
       aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                <  2:1516592/150‑1 (MQ=255)
       aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTc                                                                                                                                                <  2:52872/150‑1 (MQ=255)
        cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                               >  1:1330608/1‑150 (MQ=255)
        cTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCa                                                                                                                                               <  2:941856/150‑1 (MQ=255)
          ttCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGa                                                                                                                                             <  1:564734/150‑1 (MQ=255)
           tCTGCGGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTGCGCGTTGATTTTCAGCGGTACGTTGCGCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGAc                                                                                                                                            <  2:122558/150‑1 (MQ=255)
             tgctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAACCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTg                                                                                                                                          >  1:825385/1‑150 (MQ=255)
               ctgATTTCCCCCGGAGAGATCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCTCCGTCTTTCAGACTGCg                                                                                                                                        <  1:1116883/150‑1 (MQ=255)
                 gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTAGCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGc                                                                                                                                      >  2:1290663/1‑150 (MQ=255)
                 gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGc                                                                                                                                      >  2:247741/1‑150 (MQ=255)
                 gATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACggcggc                                                                                                                                      <  1:1431461/150‑1 (MQ=255)
                  aTTTCCCCCGGAGAGTGCGGTGATCTTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGGTGATTTTCAGCGGTACGTTGCGCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCt                                                                                                                                     <  1:888006/150‑1 (MQ=255)
                   tttCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTg                                                                                                                                    >  1:1454053/1‑150 (MQ=255)
                    ttCCCCCGGAGAGTTCGGTGATATTCTGGGTTACCGAATGACATTTCATCGCCAGCAGTTCGCGTTGATTTTCAGAGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGGGCCTTTATAGACGCCGTGTTGCAGACTGCGGCTGa                                                                                                                                   >  2:1304589/1‑150 (MQ=255)
                       ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                                                                                <  1:1290663/150‑1 (MQ=255)
                       ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                                                                                >  2:742941/1‑150 (MQ=255)
                       ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCAATCGCGCCTTTATAGCCGCCGTCTTTCAGCCTGCGGCTGATCg                                                                                                                                >  1:1360689/1‑150 (MQ=255)
                        cccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCGAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGc                                                                                                                               <  2:994889/150‑1 (MQ=255)
                         ccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTGAAAGACGGCGGCtataaagg                                                                                                                              <  2:613336/150‑8 (MQ=255)
                           ggAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCAt                                                                                                                            >  2:733043/1‑150 (MQ=255)
                               agTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTc                                                                                                                        <  1:1509436/150‑1 (MQ=255)
                                  tCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGAGTTGATTTTCAGAGGTACGTTGCTCGTCATCTTAATGAAACAAGCACATCGAGCCTTTATAGCCACAGTCTTTCATTCTGCTGCTGATCACCATATTCtta                                                                                                                     >  1:911213/1‑148 (MQ=255)
                                   cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                    <  1:1029912/150‑1 (MQ=255)
                                   cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                                                                    <  1:475883/150‑1 (MQ=255)
                                    ggTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGc                                                                                                                   <  2:498012/150‑1 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              >  2:1187804/1‑150 (MQ=255)
                                         tatTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGGCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATgg                                                                                                              <  1:1507700/150‑1 (MQ=255)
                                           ttCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaa                                                                                                            >  1:1875844/1‑150 (MQ=255)
                                             cTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGataa                                                                                                          >  2:844880/1‑150 (MQ=255)
                                               ggTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGt                                                                                                        >  2:1648232/1‑150 (MQ=255)
                                                    aCCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTgggg                                                                                                   <  1:632686/150‑1 (MQ=255)
                                                      cGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTCCTCGTCAACTTCATGAAACACGCCCATCGCGCCTTTATAGCCGCCGTCATTCAGACTGAGGCTGCTCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaa                                                                                                 >  2:1136648/1‑150 (MQ=255)
                                                             cATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                                                                                          >  1:1505855/1‑150 (MQ=255)
                                                                ttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAATCAAGCCCATCGCGCCTTTATAGACGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGttt                                                                                       >  2:151049/1‑149 (MQ=255)
                                                                ttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCACATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTATGAGCGAGGGAAAAGTTGGGGAAAAAACCGTTa                                                                                       >  2:156667/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGCCTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATcc                                                                                    >  2:663634/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAATAAACCGTTATcc                                                                                    >  2:1655834/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATcc                                                                                    >  1:645538/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATcc                                                                                    >  1:1742564/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATcc                                                                                    >  1:141359/1‑150 (MQ=255)
                                                                   aGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATcc                                                                                    <  1:298425/150‑1 (MQ=255)
                                                                    gcgcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCAATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATccc                                                                                   >  1:1409647/1‑150 (MQ=255)
                                                                       cccgcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGACACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAATTTGGGGAAAAAACCGGTATCCcggc                                                                                >  2:731773/4‑150 (MQ=255)
                                                                         agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcgg                                                                              <  1:1687310/150‑1 (MQ=255)
                                                                          gcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCACCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCggcggc                                                                             <  1:197269/150‑1 (MQ=255)
                                                                           cagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCt                                                                            >  1:868432/1‑150 (MQ=255)
                                                                             gTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                                                                          >  2:1019623/1‑150 (MQ=255)
                                                                             gTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                                                                          >  1:397490/1‑150 (MQ=255)
                                                                              ttCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTc                                                                         <  2:1379004/150‑1 (MQ=255)
                                                                                cgcgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAg                                                                       >  1:1588135/1‑150 (MQ=255)
                                                                                 gcgTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGt                                                                      >  2:740133/1‑150 (MQ=255)
                                                                                    ttGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGAt                                                                   >  1:1498387/1‑150 (MQ=255)
                                                                                       agtttCAGCGGTACGTTGCTCGGCAACGTCAGGAAACAAGCCCAGCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGCTGGCAAAGTTGGTGAAAAACCCGTTATCCCGGCTGGTTTCAGTGATGTa                                                                >  2:1577290/3‑150 (MQ=255)
                                                                                        ttttCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCATACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTaa                                                               >  2:1277710/1‑150 (MQ=255)
                                                                                           tCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCCCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGATATCCCGGCGGCTTTCAGTGATGTAAGca                                                            >  2:1664993/1‑149 (MQ=255)
                                                                                              gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGACATc                                                         >  1:984169/1‑150 (MQ=255)
                                                                                              gCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGACATc                                                         >  1:214786/1‑150 (MQ=255)
                                                                                                ggTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATccc                                                       <  2:892914/150‑1 (MQ=255)
                                                                                                 gTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCAGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCt                                                      >  1:1493155/1‑150 (MQ=255)
                                                                                                  tACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCtt                                                     >  1:1795224/1‑150 (MQ=255)
                                                                                                    cGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCtttt                                                   >  2:152060/1‑150 (MQ=255)
                                                                                                        gCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCa                                               >  2:68290/1‑150 (MQ=255)
                                                                                                        gCTCGTCAACTTCATGAAACAAGCCAATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTcc                                               >  1:893250/1‑149 (MQ=255)
                                                                                                         cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                              >  1:1237573/1‑150 (MQ=255)
                                                                                                         cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                              <  1:1648232/150‑1 (MQ=255)
                                                                                                         cTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCAc                                              <  1:502931/150‑1 (MQ=255)
                                                                                                           cGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACgg                                            >  1:1565340/1‑150 (MQ=255)
                                                                                                                aCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATcc                                       <  2:1444218/150‑1 (MQ=255)
                                                                                                                 cTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAACCCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCTTCCCTTTTTTCACGGCATCCa                                      >  2:1174719/1‑150 (MQ=255)
                                                                                                                 cTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCa                                      >  1:879274/1‑150 (MQ=255)
                                                                                                                  ttCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAg                                     >  2:481496/1‑150 (MQ=255)
                                                                                                                    cATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAggg                                   >  1:1382279/1‑150 (MQ=255)
                                                                                                                     aTGAAACAAGCCCCTCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCAGGCGGCTTTCAGTGATGTAAGACATACATTTTTTCACGGCATCCAgggg                                  >  1:1406807/1‑150 (MQ=255)
                                                                                                                        aaaCAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGa                               >  2:447626/1‑150 (MQ=255)
                                                                                                                         aaCAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGaa                              <  1:880281/150‑1 (MQ=255)
                                                                                                                          aCAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAAc                             <  1:247741/150‑1 (MQ=255)
                                                                                                                             aGCCCATCGCGCCTTTATAGCCTCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCAGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTg                          >  2:450162/1‑150 (MQ=255)
                                                                                                                              gCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCGTTTTTTCACGGCATCCAGGGGTGAACGTgg                         >  2:1708348/1‑150 (MQ=255)
                                                                                                                               cccATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGCGAAAAAACCGATATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGga                        <  2:519300/150‑1 (MQ=255)
                                                                                                                                ccATCGCGCCTTTATAGCCGCCGGCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGATAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       <  2:1625230/150‑1 (MQ=255)
                                                                                                                                ccATCGCCCCTTTATAGCCGCCGGCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGGTGTGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgag                       <  2:868758/150‑1 (MQ=255)
                                                                                                                                  aTCGCGCCTTTATAGCCTCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGGTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagag                     <  2:1026052/150‑1 (MQ=255)
                                                                                                                                   tCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGgagaga                    <  2:877442/150‑1 (MQ=255)
                                                                                                                                    cgcgCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGat                   >  1:1043710/1‑150 (MQ=255)
                                                                                                                                     gcgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGata                  >  1:1603647/1‑150 (MQ=255)
                                                                                                                                      cgcCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAATAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGatat                 >  1:459727/1‑150 (MQ=255)
                                                                                                                                           ttATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGCTGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTg            <  1:75339/150‑1 (MQ=255)
                                                                                                                                           ttATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTg            <  2:1409647/150‑1 (MQ=255)
                                                                                                                                              taGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCa         <  2:1688318/150‑1 (MQ=255)
                                                                                                                                               aGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAt        <  1:699699/150‑1 (MQ=255)
                                                                                                                                               aGCCGCAGTCTGAAAGCCGGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCAt        >  2:1288031/1‑150 (MQ=255)
                                                                                                                                                 ccgccgTCTTTCAGCCTGCGGCTGATCGCCATGGTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTa      >  2:1519572/1‑150 (MQ=255)
                                                                                                                                                  cgccgTCTTTCAGACTGCGGCTGCTCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAACCCGTTATACCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGCACGTGGAGAGATATCTTTGCCATTaa     >  1:146123/1‑150 (MQ=255)
                                                                                                                                                     cgTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGACGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTATCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGGGGAGAGATATATTTGGCATTAAGAc  >  1:351039/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTTTTTTCACGGCATCCAGGGGTGAACGTGGAGAGATATCTTTGCCATTAAGAC  >  NC_000913/4309721‑4310019

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: