Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,989 G→C 5.3% P168R (CCG→CGG)  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,9890GC 5.3% 448.6 / 11.4 132P168R (CCG→CGG) nfeFNADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base G (68/57);  new base C (2/5);  total (70/62)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 134 total aligned reads displayed.

AGTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTT  >  NC_000913/3215840‑3216136
                                                                                                                                                     |                                                                                                                                                   
aGTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAACACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCg                                                                                                                                                     >  2:513491/1‑150 (MQ=255)
 gTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCgg                                                                                                                                                    >  1:1205898/1‑150 (MQ=255)
  tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                                   <  1:43570/150‑1 (MQ=255)
  tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                                   <  1:1041197/150‑1 (MQ=255)
  tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                                   <  1:1550313/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                                                                                                                                            <  2:203926/150‑1 (MQ=255)
                gATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGttt                                                                                                                                     <  1:141696/150‑1 (MQ=255)
                  tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGc                                                                                                                                   >  2:1388169/1‑150 (MQ=255)
                   ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCt                                                                                                                                  <  2:621231/150‑1 (MQ=255)
                   ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                                  <  1:1536062/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)
                              ccAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCa                                                                                                                       <  2:1013234/150‑1 (MQ=255)
                                      gcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggcgg                                                                                                               <  2:254828/150‑1 (MQ=255)
                                         tCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGgcgc                                                                                                            <  1:2455538/150‑1 (MQ=255)
                                           tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                          >  2:1569546/1‑150 (MQ=255)
                                           tACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCaa                                                                                                          >  1:453191/1‑150 (MQ=255)
                                            aCCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAATCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAAt                                                                                                         >  2:685593/1‑150 (MQ=255)
                                             ccGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATg                                                                                                        <  2:650048/150‑1 (MQ=255)
                                               gCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCa                                                                                                      >  1:1078211/1‑150 (MQ=255)
                                                 cTGCTCATCATGTGCCAGCCATTTGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAgg                                                                                                    <  2:1552410/150‑1 (MQ=255)
                                                    cTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAt                                                                                                 >  1:299275/1‑150 (MQ=255)
                                                      catcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTc                                                                                               >  1:1740184/1‑150 (MQ=255)
                                                       atcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTcc                                                                                              <  1:572971/150‑1 (MQ=255)
                                                         catGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCgg                                                                                            >  1:713734/1‑150 (MQ=255)
                                                          atGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGa                                                                                           <  1:2482551/150‑1 (MQ=255)
                                                           tgtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTACTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAt                                                                                          >  1:1976414/1‑150 (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:891502/1‑150 (MQ=255)
                                                              gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                       >  2:2179077/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)
                                                                      ttCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                               <  1:909461/150‑1 (MQ=255)
                                                                       tCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTACCGGCAAGTTTGCTCAACGTTTCCCGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACa                                                                              >  1:2448542/1‑150 (MQ=255)
                                                                        cGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACAt                                                                             >  2:2032111/1‑150 (MQ=255)
                                                                        cGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACAt                                                                             >  2:260191/1‑150 (MQ=255)
                                                                           tatTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACa                                                                          >  2:1955754/1‑150 (MQ=255)
                                                                            atTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                         >  1:1647163/1‑150 (MQ=255)
                                                                            atTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                         <  2:1521314/150‑1 (MQ=255)
                                                                              tAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCt                                                                       >  1:1586365/1‑150 (MQ=255)
                                                                               aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                      >  1:1322885/1‑150 (MQ=255)
                                                                                   ccATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAc                                                                  >  2:1652769/1‑150 (MQ=255)
                                                                                   ccATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAc                                                                  >  1:29601/1‑150 (MQ=255)
                                                                                    cATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcg                                                                 >  2:1334095/1‑150 (MQ=255)
                                                                                       cAAGGTGCGCGAGATAATCCTGACAGGCGTTATCACGCACGCTAACCAGCGCACTAACTTGCGGTTTAACTGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCTGATTCATCGCAGACATACAGCTGATACGAGt                                                              >  2:1380164/1‑150 (MQ=255)
                                                                                        aaGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTa                                                             >  2:2226861/1‑150 (MQ=255)
                                                                                         aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                            >  2:1283619/1‑150 (MQ=255)
                                                                                          ggTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTACCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCCACGTTTCCAGGCGGCGGCGCAATGAAGGCATTCCGGATTCATCGCAGCCATACAGCTGATACGCGTAAt                                                           >  2:553056/1‑150 (MQ=255)
                                                                                           gTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATc                                                          <  1:636137/150‑1 (MQ=255)
                                                                                           gTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATc                                                          <  2:299275/150‑1 (MQ=255)
                                                                                            tGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCt                                                         <  1:766696/150‑1 (MQ=255)
                                                                                             gcgcGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCtt                                                        >  1:19665/1‑150 (MQ=255)
                                                                                             gcgcGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCtt                                                        >  1:2190772/1‑150 (MQ=255)
                                                                                             gcgcGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCtt                                                        >  1:1342323/1‑150 (MQ=255)
                                                                                              cgcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGAATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTc                                                       >  2:786456/1‑150 (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)
                                                                                                          cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAg                                           <  2:1384821/150‑1 (MQ=255)
                                                                                                           tGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGc                                          >  1:1651422/1‑150 (MQ=255)
                                                                                                                ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                     >  1:444609/1‑150 (MQ=255)
                                                                                                                ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                     >  1:1893248/1‑150 (MQ=255)
                                                                                                                   gTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcg                                  <  2:702518/150‑1 (MQ=255)
                                                                                                                   gTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcg                                  >  2:1690591/1‑150 (MQ=255)
                                                                                                                    ttATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcgc                                 <  2:1444158/150‑1 (MQ=255)
                                                                                                                          cGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATAc                                                                           <  1:2141084/102‑1 (MQ=255)
                                                                                                                              cGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCa                       >  2:1993978/1‑150 (MQ=255)
                                                                                                                               gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACGTGCCAc                      <  2:1322885/150‑1 (MQ=255)
                                                                                                                                cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                     >  2:113432/1‑150 (MQ=255)
                                                                                                                                cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                     >  1:581377/1‑150 (MQ=255)
                                                                                                                                  aaCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGt                   >  2:2263088/1‑150 (MQ=255)
                                                                                                                                    ccAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTaa                 >  1:873111/1‑150 (MQ=255)
                                                                                                                                       gcgcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGtt              <  2:651850/150‑1 (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)
                                                                                                                                            cTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCg         >  1:560469/1‑150 (MQ=255)
                                                                                                                                               aCTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGccc      <  1:1166346/150‑1 (MQ=255)
                                                                                                                                                cTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCg     >  1:2287666/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTT  >  NC_000913/3215840‑3216136

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: ‑