Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,215,989 G→C 10.5% 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,9890GC10.5% 217.7 / 16.3 76P168R (CCG→CGG) nfeFNADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base G (30/38);  new base C (2/6);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

GTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCC  >  NC_000913/3215841‑3216132
                                                                                                                                                    |                                                                                                                                               
gTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCTATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCATGAAATGCGTTATCCCGCACGCTAACCAGCGAACTAACTTTCgg                                                                                                                                                <  2:639872/150‑1 (MQ=255)
 tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                               <  2:509645/150‑1 (MQ=255)
 tAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGt                                                                                                                                               <  2:20102/150‑1 (MQ=255)
      atcGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATTCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAAc                                                                                                                                          <  1:1155200/150‑1 (MQ=255)
        cGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAGCgg                                                                                                                                        <  1:270458/150‑1 (MQ=255)
             gggATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGt                                                                                                                                   >  1:69191/1‑150 (MQ=255)
              ggATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGtt                                                                                                                                  <  1:1249686/150‑1 (MQ=255)
               gATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGttt                                                                                                                                 <  1:947419/150‑1 (MQ=255)
               gATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGttt                                                                                                                                 <  2:1313090/150‑1 (MQ=255)
                aTTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTg                                                                                                                                <  2:12746/150‑1 (MQ=255)
                 tttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGc                                                                                                                               <  1:1264618/150‑1 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCt                                                                                                                              <  2:689944/150‑1 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCt                                                                                                                              <  1:1558657/150‑1 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                              <  1:871819/150‑5 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                              <  2:1374977/150‑5 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                              <  1:818428/150‑5 (MQ=255)
                  ttGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtg                                                                                                                              <  2:375709/150‑5 (MQ=255)
                   tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTc                                                                                                                             <  1:1079045/150‑1 (MQ=255)
                   tGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCCGTTAAACCGCAAGttagtgc                                                                                                                             <  1:1200983/150‑6 (MQ=255)
                    gCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCa                                                                                                                            >  1:658287/1‑150 (MQ=255)
                       tCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                         <  1:71516/150‑1 (MQ=255)
                       tCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACg                                                                                                                         >  1:1215744/1‑150 (MQ=255)
                        cTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGt                                                                                                                        <  2:569480/150‑1 (MQ=255)
                                             cGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGc                                                                                                   <  2:824079/150‑1 (MQ=255)
                                               ccTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAg                                                                                                 <  1:1376514/150‑1 (MQ=255)
                                                cTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCGAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACCAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAgg                                                                                                <  2:1285605/150‑1 (MQ=255)
                                                 tGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGc                                                                                               <  2:152705/150‑1 (MQ=255)
                                                  gCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCa                                                                                              >  2:275579/1‑150 (MQ=255)
                                                    tcatcaTGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCAtt                                                                                            >  2:1238193/1‑150 (MQ=255)
                                                      atcatGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTcc                                                                                          >  1:311735/1‑150 (MQ=255)
                                                          tgtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAt                                                                                      <  2:1384470/150‑1 (MQ=255)
                                                          tgtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAt                                                                                      <  1:191767/150‑1 (MQ=255)
                                                           ttgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGCCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTGTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAtt                                                                                     <  2:1592888/149‑1 (MQ=255)
                                                           gtgCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACTTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAtt                                                                                     >  1:640629/1‑150 (MQ=255)
                                                           gtgCCAGCAATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGTCGATATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGAtt                                                                                     >  2:1276301/1‑150 (MQ=255)
                                                             gccagccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCa                                                                                   >  2:445894/1‑150 (MQ=255)
                                                                  ccaTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCa                                                                              <  1:1487515/150‑1 (MQ=255)
                                                                   caTTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAg                                                                             <  1:1238193/150‑1 (MQ=255)
                                                                     ttCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAc                                                                           <  2:40900/150‑1 (MQ=255)
                                                                      tCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACa                                                                          >  1:325392/1‑150 (MQ=255)
                                                                       cGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACAt                                                                         <  1:445894/150‑1 (MQ=255)
                                                                        gATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATa                                                                        >  1:508629/1‑150 (MQ=255)
                                                                          tatTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACa                                                                      >  1:1552928/1‑150 (MQ=255)
                                                                            ttAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGc                                                                    >  2:1602791/1‑150 (MQ=255)
                                                                             tAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCt                                                                   <  2:1136881/150‑1 (MQ=255)
                                                                              aaaaCCATCAAGGTGCGCTAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                  <  1:462752/150‑1 (MQ=255)
                                                                              aaaaCCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTg                                                                  >  2:1230555/1‑150 (MQ=255)
                                                                                      cAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGt                                                          <  2:508629/150‑1 (MQ=255)
                                                                                        aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGTGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                        >  2:373753/1‑150 (MQ=255)
                                                                                        aGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCGGGATTCATCGCAGACATACAGCTGATACGCGTaa                                                        >  1:1323877/1‑150 (MQ=255)
                                                                                           tGCGCGAGAAAATCCAGACAGGCGTTATCCCGAACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGAGTAGCAATTCAGCCATTCCGGATTCATCGCAGACATACAGATGATACGTGTAATCt                                                     >  2:623991/1‑150 (MQ=255)
                                                                                              gcgAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTcc                                                  <  2:74234/150‑1 (MQ=255)
                                                                                                 agaTAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACTGCTGATAGGCGTAATCTTCCGGc                                               <  2:1277924/150‑1 (MQ=255)
                                                                                                   aTAATCCTGACAGGCGTTATCCGGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCGAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcac                                             <  1:1439621/150‑1 (MQ=255)
                                                                                                    tAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacc                                            <  1:823898/150‑1 (MQ=255)
                                                                                                     aaTCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGcacca                                           <  1:822290/150‑1 (MQ=255)
                                                                                                       tCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCaccacc                                         >  1:1241298/1‑150 (MQ=255)
                                                                                                          tGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGc                                      >  2:848568/1‑150 (MQ=255)
                                                                                                               ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                 >  2:523409/1‑150 (MQ=255)
                                                                                                               ggCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGCTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGAGATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAAcc                                 <  1:194059/150‑1 (MQ=255)
                                                                                                                     aTCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCgg                           >  1:1454476/1‑150 (MQ=255)
                                                                                                                       cccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATCCAGCCGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCggcc                         >  1:710185/1‑148 (MQ=255)
                                                                                                                       cccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAc                         >  2:881925/1‑150 (MQ=255)
                                                                                                                       cccGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGAc                         <  2:1241298/150‑1 (MQ=255)
                                                                                                                          gCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACTTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTg                      <  2:640629/150‑1 (MQ=255)
                                                                                                                            aCGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGcc                    <  2:1091424/150‑1 (MQ=255)
                                                                                                                                 aaCCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGt               >  1:403865/1‑150 (MQ=255)
                                                                                                                                   ccAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGCCCTGCCACCGTaa             >  2:1324669/1‑150 (MQ=255)
                                                                                                                                   ccAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTaa             >  1:1505164/1‑150 (MQ=255)
                                                                                                                                      gcgcACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTGCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGtt          <  2:1403319/150‑1 (MQ=255)
                                                                                                                                         agcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTAt       >  1:165503/5‑150 (MQ=255)
                                                                                                                                          gcaaaCTTGCCGTTAAACCGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATc      >  1:812138/4‑150 (MQ=255)
                                                                                                                                           cTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCg     >  1:1425103/1‑150 (MQ=255)
                                                                                                                                           cTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCg     >  1:1403033/1‑150 (MQ=255)
                                                                                                                                             aaCTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGcc   >  1:1630057/1‑150 (MQ=255)
                                                                                                                                                  gCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGccc  >  1:1413234/1‑146 (MQ=255)
                                                                                                                                                  gCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGccc  <  2:1413234/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTAATCATCGGCAGGGATTTGCATCTGCGCCAGACGCGCATCTACCGCCTGCTCATCATGTGCCAGCCATTCGATATTAAAACCATCAAGGTGCGCGAGATAATCCTGACAGGCGTTATCCCGCACGCTAACCAGCGCACTAACTTGCGGTTTAACGGCAAGTTTGCTCAACGTTTCCAGGCGGCGGCGCAATGCAGGCATTCCGGATTCATCGCAGACATACAGCTGATACGCGTAATCTTCCGGCACCACCAGCGAACCGCGCGGACCTGCCACCGTAAGTTTATCGCCC  >  NC_000913/3215841‑3216132

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: