Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,993 T→A 6.1% K167* (AAA→TAA)  nfeF ← NADPH‑dependent ferric chelate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,215,9930TA 6.1% 405.9 / 10.5 129K167* (AAA→TAA) nfeFNADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base T (68/53);  new base A (2/6);  total (70/59)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 131 total aligned reads displayed.

ATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCG  >  NC_000913/3215844‑3216139
                                                                                                                                                     |                                                                                                                                                  
atcatcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGttt                                                                                                                                                    <  1:716687/150‑1 (MQ=255)
 tcatcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTa                                                                                                                                                   <  1:1154059/150‑1 (MQ=255)
   atcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAc                                                                                                                                                 <  2:1007666/150‑1 (MQ=255)
    tcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACg                                                                                                                                                <  2:1254934/150‑1 (MQ=255)
    tcGGCAGGGATTTGCATCTCCGCCAGACGCGCATCTACCGCCTGCGCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACg                                                                                                                                                <  2:1246363/150‑1 (MQ=255)
     cGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACgg                                                                                                                                               <  1:953268/150‑1 (MQ=255)
     cGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACgg                                                                                                                                               <  2:203926/150‑1 (MQ=255)
     cGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACgg                                                                                                                                               <  1:593305/150‑1 (MQ=255)
            gATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGttt                                                                                                                                        <  1:141696/150‑1 (MQ=255)
             aTTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttag                                                                                                                                       <  1:47855/150‑3 (MQ=255)
              tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGc                                                                                                                                      >  2:1388169/1‑150 (MQ=255)
               ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCt                                                                                                                                     <  2:621231/150‑1 (MQ=255)
               ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                                     <  2:1022376/150‑5 (MQ=255)
               ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                                     <  1:1561489/150‑5 (MQ=255)
                tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtgc                                                                                                                                    <  1:1708207/150‑6 (MQ=255)
                 gCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCa                                                                                                                                   >  1:651850/1‑150 (MQ=255)
                  cATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCaa                                                                                                                                  >  1:1051304/1‑150 (MQ=255)
                       gcgcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGAGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGttt                                                                                                                             <  2:688486/150‑1 (MQ=255)
                       gcgcCAGACGCGCATCTACCGCCTGCTCATCATGTGCAAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGttt                                                                                                                             >  2:2093732/1‑150 (MQ=255)
                        cgcCAGCCGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTGTCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTc                                                                                                                            >  1:1604122/1‑150 (MQ=255)
                                gcgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggc                                                                                                                    <  2:2423692/150‑1 (MQ=255)
                                  gcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggcgg                                                                                                                  <  2:254828/150‑1 (MQ=255)
                                     tCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGgcgc                                                                                                               <  1:2455538/150‑1 (MQ=255)
                                       tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                             >  1:453191/1‑150 (MQ=255)
                                       tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                             >  2:1569546/1‑150 (MQ=255)
                                        aCCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAATCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAAt                                                                                                            >  2:685593/1‑150 (MQ=255)
                                         ccGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATg                                                                                                           <  2:650048/150‑1 (MQ=255)
                                         ccGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATg                                                                                                           <  1:1286805/150‑1 (MQ=255)
                                         ccGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATg                                                                                                           >  1:1444158/1‑150 (MQ=255)
                                           gCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCa                                                                                                         >  1:1078211/1‑150 (MQ=255)
                                             cTGCTCATCATGTGCCAGCCATTTGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAgg                                                                                                       <  2:1552410/150‑1 (MQ=255)
                                                cTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAt                                                                                                    >  2:1939454/1‑150 (MQ=255)
                                                cTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAt                                                                                                    >  1:299275/1‑150 (MQ=255)
                                                  catcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTc                                                                                                  >  1:1748907/1‑150 (MQ=255)
                                                   atcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTcc                                                                                                 <  1:572971/150‑1 (MQ=255)
                                                    tcatGCGCCAGCCACTCGATATTAAAACCATCAACGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCg                                                                                                <  1:1980086/150‑1 (MQ=255)
                                                      atGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGa                                                                                              <  1:2482551/150‑1 (MQ=255)
                                                       tgtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAt                                                                                             <  1:1673695/150‑1 (MQ=255)
                                                         tgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTc                                                                                           >  1:702518/1‑150 (MQ=255)
                                                         tgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTc                                                                                           >  1:1521314/1‑150 (MQ=255)
                                                          gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                          >  1:63316/1‑150 (MQ=255)
                                                          gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                          >  2:425518/1‑150 (MQ=255)
                                                          gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                          >  2:891502/1‑150 (MQ=255)
                                                            cagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATc                                                                                        <  2:518481/150‑1 (MQ=255)
                                                            cagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATc                                                                                        <  1:2270993/150‑1 (MQ=255)
                                                            cagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATc                                                                                        >  1:2332188/1‑150 (MQ=255)
                                                              gccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGc                                                                                      <  2:1927581/150‑1 (MQ=255)
                                                              gccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTAGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGc                                                                                      >  2:1057420/1‑150 (MQ=255)
                                                               ccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCa                                                                                     <  2:1051304/150‑1 (MQ=255)
                                                                caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                                    <  1:2376139/150‑1 (MQ=255)
                                                                caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                                    <  1:2367380/150‑1 (MQ=255)
                                                                  ttCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                                  <  1:99813/150‑1 (MQ=255)
                                                                   tCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTACCGGCAAGTTTGCTCAACGTTTCCCGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACa                                                                                 >  1:2448542/1‑150 (MQ=255)
                                                                     gATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAAGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATa                                                                               <  2:662741/150‑1 (MQ=255)
                                                                        atTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                            <  2:1521314/150‑1 (MQ=255)
                                                                        atTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                            >  1:1647163/1‑150 (MQ=255)
                                                                           aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                         >  1:1322885/1‑150 (MQ=255)
                                                                               ccATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAc                                                                     >  1:29601/1‑150 (MQ=255)
                                                                               ccATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAc                                                                     >  2:1652769/1‑150 (MQ=255)
                                                                                cATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcg                                                                    >  2:1334095/1‑150 (MQ=255)
                                                                                   cAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGt                                                                 <  1:1117356/150‑1 (MQ=255)
                                                                                   cAAGGTGCGCGAGATAATCCTGACAGGCGTTATCACGCACGCTAACCAGCGCACTAACTTGCGGTTTAACTGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCTGATTCATCGCAGACATACAGCTGATACGAGt                                                                 >  2:1380164/1‑150 (MQ=255)
                                                                                    aaGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTa                                                                >  2:2226861/1‑150 (MQ=255)
                                                                                     aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                               >  2:1283619/1‑150 (MQ=255)
                                                                                       gTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATc                                                             <  2:299275/150‑1 (MQ=255)
                                                                                        tGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCt                                                            <  1:766696/150‑1 (MQ=255)
                                                                                         gcgcGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCtt                                                           >  1:19665/1‑150 (MQ=255)
                                                                                          cgcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGAATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTc                                                          >  2:786456/1‑150 (MQ=255)
                                                                                           gcgAGATAATCCTGACGGGGGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                         >  2:616749/1‑150 (MQ=255)
                                                                                              agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGc                                                      <  1:2257796/150‑1 (MQ=255)
                                                                                              agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTACTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGc                                                      <  2:1976414/150‑1 (MQ=255)
                                                                                                 tAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacc                                                   <  2:1908155/150‑1 (MQ=255)
                                                                                                 tAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacc                                                   <  1:156331/150‑1 (MQ=255)
                                                                                                  aaTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacca                                                  <  1:335514/150‑1 (MQ=255)
                                                                                                  aaTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTCTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacca                                                  >  2:1484083/1‑150 (MQ=255)
                                                                                                    tcatGACAGGCGTTATCACGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTATCCAGGCGGCGGCGCAATCAAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCcccacc                                                >  2:1578188/4‑150 (MQ=255)
                                                                                                    tCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCaccacc                                                >  1:1628935/1‑150 (MQ=255)
                                                                                                     ccTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCAccacca                                               >  2:706792/1‑150 (MQ=255)
                                                                                                      cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAg                                              <  2:1384821/150‑1 (MQ=255)
                                                                                                            ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                        >  1:444609/1‑150 (MQ=255)
                                                                                                            ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                        >  1:1893248/1‑150 (MQ=255)
                                                                                                              cGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgc                                      >  2:1681373/1‑150 (MQ=255)
                                                                                                                ttATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcgc                                    <  2:1444158/150‑1 (MQ=255)
                                                                                                                     ccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAcc                               <  1:685593/150‑1 (MQ=255)
                                                                                                                      cGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATAc                                                                              >  2:2141084/1‑102 (MQ=255)
                                                                                                                      cGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATAc                                                                              <  1:2141084/102‑1 (MQ=255)
                                                                                                                          cGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCa                          >  2:2002597/1‑150 (MQ=255)
                                                                                                                          cGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCa                          >  2:1993978/1‑150 (MQ=255)
                                                                                                                           gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACGTGCCAc                         <  2:1322885/150‑1 (MQ=255)
                                                                                                                            cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                        >  1:581377/1‑150 (MQ=255)
                                                                                                                            cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                        >  2:113432/1‑150 (MQ=255)
                                                                                                                              aaCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGt                      >  2:2263088/1‑150 (MQ=255)
                                                                                                                                ccAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTaa                    >  1:873111/1‑150 (MQ=255)
                                                                                                                                     gcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTa               <  1:2093732/150‑1 (MQ=255)
                                                                                                                                     gcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTa               <  1:1681373/150‑1 (MQ=255)
                                                                                                                                       gcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATc             >  1:2481016/4‑150 (MQ=255)
                                                                                                                                        caaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCg            >  1:1593378/3‑150 (MQ=255)
                                                                                                                                            cTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCg        >  1:2287666/1‑150 (MQ=255)
                                                                                                                                                  gTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTgcg  <  1:1559522/150‑1 (MQ=255)
                                                                                                                                                  gTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTgcg  >  1:1331058/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGCG  >  NC_000913/3215844‑3216139

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: