Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,989 G→C 5.2% 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.2% 333.2 / 7.1 97P168R (CCG→CGG) nfeFNADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base G (47/45);  new base C (2/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGC  >  NC_000913/3215842‑3216138
                                                                                                                                                   |                                                                                                                                                     
tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                                     <  1:32739/150‑1 (MQ=255)
     atcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAc                                                                                                                                                >  1:791124/1‑150 (MQ=255)
      tcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACg                                                                                                                                               <  1:1330319/150‑1 (MQ=255)
       cGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAAccg                                                                                                                                              <  1:1066748/150‑3 (MQ=255)
        ggCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGc                                                                                                                                             >  2:723891/1‑150 (MQ=255)
         gCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCa                                                                                                                                            >  1:279847/1‑150 (MQ=255)
           aGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAg                                                                                                                                          <  1:999840/150‑1 (MQ=255)
               aTTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttag                                                                                                                                      <  1:1549993/150‑3 (MQ=255)
                tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagt                                                                                                                                     <  2:1638736/150‑4 (MQ=255)
                  tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTc                                                                                                                                   <  2:1768908/150‑1 (MQ=255)
                      tCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                               <  2:1802063/150‑1 (MQ=255)
                      tCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                               <  1:1424825/150‑1 (MQ=255)
                       cTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGt                                                                                                                              >  1:685350/1‑150 (MQ=255)
                         gcgcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGttt                                                                                                                            >  2:1738663/1‑150 (MQ=255)
                          cgcCAGACGCGCATCTACCGCCTGCTCGTCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCCACGTTTc                                                                                                                           >  1:1788939/1‑150 (MQ=255)
                           gcCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTcc                                                                                                                          <  1:723891/150‑1 (MQ=255)
                            ccAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCa                                                                                                                         >  1:1991820/1‑150 (MQ=255)
                             cAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAg                                                                                                                        >  1:365251/1‑150 (MQ=255)
                                  gcgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggc                                                                                                                   <  2:395030/150‑1 (MQ=255)
                                  gcgcATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggc                                                                                                                   <  1:647636/150‑1 (MQ=255)
                                     cATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAggcggcggc                                                                                                                >  2:1102386/1‑150 (MQ=255)
                                       tCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGgcgc                                                                                                              <  1:903831/150‑1 (MQ=255)
                                       tCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGgcgc                                                                                                              <  1:896800/150‑1 (MQ=255)
                                          aCCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAAt                                                                                                           >  2:852806/1‑150 (MQ=255)
                                              ccTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGTGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGGGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAg                                                                                                       <  1:1986066/150‑1 (MQ=255)
                                              ccTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAg                                                                                                       <  2:791124/150‑1 (MQ=255)
                                                  cTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAt                                                                                                   >  1:1549557/1‑150 (MQ=255)
                                                      tcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCg                                                                                               <  1:269623/150‑1 (MQ=255)
                                                      tcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCg                                                                                               >  2:1347420/1‑150 (MQ=255)
                                                        atGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGa                                                                                             <  2:112436/150‑1 (MQ=255)
                                                          gtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAtt                                                                                           <  2:1567373/150‑1 (MQ=255)
                                                            gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCGGGATTCa                                                                                         >  1:1337826/1‑150 (MQ=255)
                                                            gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                         >  1:1931464/1‑150 (MQ=255)
                                                            gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                         >  1:131269/1‑150 (MQ=255)
                                                               agccaTTCGATATTAAAACCATCAATGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCg                                                                                      >  2:629773/1‑150 (MQ=255)
                                                                 ccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCa                                                                                    >  1:1938790/1‑150 (MQ=255)
                                                                  caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                                   <  1:1509892/150‑1 (MQ=255)
                                                                     tCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACa                                                                                >  1:324490/1‑150 (MQ=255)
                                                                         tatTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACa                                                                            >  1:41385/1‑150 (MQ=255)
                                                                          atTAAAACTATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                           >  2:628059/1‑150 (MQ=255)
                                                                          atTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAg                                                                           >  1:1761830/1‑150 (MQ=255)
                                                                           ttAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGc                                                                          >  2:123328/1‑150 (MQ=255)
                                                                             aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATATAGCTg                                                                        <  1:1738663/150‑1 (MQ=255)
                                                                             aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                        >  1:1758838/1‑150 (MQ=255)
                                                                              aaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGa                                                                       >  2:1689249/1‑150 (MQ=255)
                                                                              aaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGa                                                                       >  2:1682294/1‑150 (MQ=255)
                                                                               aaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGAt                                                                      <  1:300589/150‑1 (MQ=255)
                                                                               aaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGAt                                                                      <  2:292722/150‑1 (MQ=255)
                                                                                  cATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcg                                                                   >  2:1102479/1‑150 (MQ=255)
                                                                                   aTCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATAcgc                                                                  >  1:368598/1‑150 (MQ=255)
                                                                                       aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                              >  1:1385573/1‑150 (MQ=255)
                                                                                       aGGTGCGCGAGATAAACCTGACAGGCGTTATCCCGAACGCAAACCAGAGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                              >  2:907664/1‑150 (MQ=255)
                                                                                         gTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATc                                                            <  2:685350/150‑1 (MQ=255)
                                                                                         gTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATc                                                            >  2:1624745/1‑150 (MQ=255)
                                                                                             gcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                        <  2:1117757/150‑1 (MQ=255)
                                                                                             gcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                        <  1:852806/150‑1 (MQ=255)
                                                                                             gcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                        <  2:279847/150‑1 (MQ=255)
                                                                                                 gaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGca                                                    <  1:1196308/150‑1 (MQ=255)
                                                                                                   tAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacc                                                  >  1:1889267/1‑150 (MQ=255)
                                                                                                    aaTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacca                                                 <  2:1508781/150‑1 (MQ=255)
                                                                                                      tCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCaccacc                                               <  2:1454125/150‑1 (MQ=255)
                                                                                                      tCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCaccacc                                               <  1:1599906/150‑1 (MQ=255)
                                                                                                       ccTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCAccacca                                              >  2:1512067/1‑150 (MQ=255)
                                                                                                        cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                                 <  1:1484396/114‑1 (MQ=255)
                                                                                                        cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                                 >  2:1484396/1‑114 (MQ=255)
                                                                                                        cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAg                                             >  1:1804542/1‑150 (MQ=255)
                                                                                                        cTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAg                                             >  1:30020/1‑150 (MQ=255)
                                                                                                         tGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCATACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGc                                            <  2:1368155/150‑1 (MQ=255)
                                                                                                          gACAGGCGTTATCCCGCACGCTAACCAGCGCACTGACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCg                                           >  1:1132779/1‑150 (MQ=255)
                                                                                                          gACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGAAACGCGTAATCTTCCGGCACCACCAGCg                                           <  1:629773/150‑1 (MQ=255)
                                                                                                           aCAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGa                                          <  1:1682294/150‑1 (MQ=255)
                                                                                                           aCAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGa                                          <  1:1689249/150‑1 (MQ=255)
                                                                                                                cGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgc                                     >  2:398372/1‑150 (MQ=255)
                                                                                                                 gTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcg                                    <  2:1991820/150‑1 (MQ=255)
                                                                                                                   tATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACcgcgcg                                  >  2:1879648/1‑150 (MQ=255)
                                                                                                                     tCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGa                                >  1:1034993/1‑150 (MQ=255)
                                                                                                                      cccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAc                               >  2:1408301/1‑150 (MQ=255)
                                                                                                                       ccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAcc                              >  1:81814/1‑150 (MQ=255)
                                                                                                                        cGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCt                             >  1:240825/1‑150 (MQ=255)
                                                                                                                         gCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                            <  2:1337826/150‑1 (MQ=255)
                                                                                                                           aCGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGcc                          <  1:1102386/150‑1 (MQ=255)
                                                                                                                            cGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCa                         >  1:723440/1‑150 (MQ=255)
                                                                                                                             gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAc                        <  2:800513/150‑1 (MQ=255)
                                                                                                                             gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAc                        <  2:1010748/150‑1 (MQ=255)
                                                                                                                             gCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAc                        <  2:1034993/150‑1 (MQ=255)
                                                                                                                              cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                       <  1:1102479/150‑1 (MQ=255)
                                                                                                                              cTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCAcc                       >  2:961111/1‑150 (MQ=255)
                                                                                                                               tAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATTGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCg                      <  1:1925842/150‑1 (MQ=255)
                                                                                                                                aaCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGt                     <  2:131269/150‑1 (MQ=255)
                                                                                                                                 aCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTa                    >  1:1376987/1‑150 (MQ=255)
                                                                                                                                  ccAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGACCCGCGCGGACCTGCCACCGTaa                   >  1:1757721/1‑150 (MQ=255)
                                                                                                                                     gcgcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGtt                <  1:123328/150‑1 (MQ=255)
                                                                                                                                      tgagcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGttt               >  1:1860159/7‑150 (MQ=255)
                                                                                                                                       gcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTa              >  2:1154642/1‑150 (MQ=255)
                                                                                                                                        agcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTAt             >  2:260971/5‑150 (MQ=255)
                                                                                                                                             aCTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGccc        <  1:907664/150‑1 (MQ=255)
                                                                                                                                                tGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGt     <  2:1376987/150‑1 (MQ=255)
                                                                                                                                                   ggTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCTTAAGTTTATCGCCCGGTTgc  >  1:56645/1‑150 (MQ=255)
                                                                                                                                                   ggTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTgc  >  1:234028/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCCGGTTGC  >  NC_000913/3215842‑3216138

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: