Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,511,665 A→G G250G (GGT→GGC fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,6650AG99.6% 967.2 / ‑4.8 237G250G (GGT→GGCfecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (135/101);  minor base T (0/1);  total (135/102)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 239 total aligned reads displayed.

GCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACAC  >  NC_000913/4511567‑4511762
                                                                                                  |                                                                                                 
gCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGCTGCTACGCATCATATGCGGCACCACGAGGcc                                                                                                 >  1:8337114/1‑101 (MQ=255)
gCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGcc                                                                                                 >  1:8337111/1‑101 (MQ=255)
 cAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCa                                                                                                >  2:8337116/1‑101 (MQ=255)
  aaCCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCaa                                                                                               >  2:8337122/1‑101 (MQ=255)
  aaCCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCaa                                                                                               >  1:8337119/1‑101 (MQ=255)
  aaCCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCaa                                                                                               >  2:8337118/1‑101 (MQ=255)
    ccaccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATa                                                                                             <  1:8336458/101‑1 (MQ=255)
     caccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGTCCAATaa                                                                                            <  2:8336711/101‑1 (MQ=255)
      accaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATaaa                                                                                           >  1:8337126/1‑101 (MQ=255)
       ccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATaaa                                                                                           <  2:8336585/100‑1 (MQ=255)
       ccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAg                                                                                          >  2:8337128/1‑101 (MQ=255)
       ccaACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAg                                                                                          <  2:8336335/101‑1 (MQ=255)
        caACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGc                                                                                         >  1:8337130/1‑101 (MQ=255)
         aaCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCt                                                                                        <  1:8336476/101‑1 (MQ=255)
         aaCAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCt                                                                                        <  1:8336779/101‑1 (MQ=255)
          aCAGCAACACACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTa                                                                                       >  1:8337134/1‑101 (MQ=255)
           cAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTaa                                                                                      <  2:8336406/101‑1 (MQ=255)
           cAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTaa                                                                                      <  1:8336858/101‑1 (MQ=255)
            aGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAAt                                                                                     <  2:8336397/101‑1 (MQ=255)
              cAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCg                                                                                   <  2:8336652/101‑1 (MQ=255)
              cAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCg                                                                                   >  2:8337136/1‑101 (MQ=255)
               aaCGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCgg                                                                                  >  1:8337139/1‑101 (MQ=255)
                aCGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCggg                                                                                 <  2:8336481/101‑1 (MQ=255)
                 cGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGc                                                                                >  1:8337141/1‑101 (MQ=255)
                 cGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGc                                                                                <  1:8336337/101‑1 (MQ=255)
                  gCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGcc                                                                               >  2:8337142/1‑101 (MQ=255)
                    aCCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGc                                                                             <  2:8336890/101‑1 (MQ=255)
                    aCCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGc                                                                             >  2:8337145/1‑101 (MQ=255)
                      cTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAg                                                                           <  1:8336374/101‑1 (MQ=255)
                      cTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAg                                                                           >  1:8337149/1‑101 (MQ=255)
                      cTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAg                                                                           <  1:8336514/101‑1 (MQ=255)
                      cTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAg                                                                           <  2:8336951/101‑1 (MQ=255)
                       tGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAgg                                                                          >  2:8337151/1‑101 (MQ=255)
                        gTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATc                                                                                    <  2:8336575/90‑1 (MQ=255)
                          cAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAggcgg                                                                       >  2:8337153/1‑101 (MQ=255)
                           aGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAggcggc                                                                      <  1:8336764/101‑1 (MQ=255)
                            gggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGcc                                                                     <  1:8336414/101‑1 (MQ=255)
                             ggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAggcggc                                                                      <  1:8336685/99‑1 (MQ=255)
                             ggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCa                                                                    <  2:8336627/101‑1 (MQ=255)
                             ggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCa                                                                    >  1:8337155/1‑101 (MQ=255)
                             ggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCa                                                                    <  2:8336515/101‑1 (MQ=255)
                             ggCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCa                                                                    <  2:8336433/101‑1 (MQ=255)
                              gCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCAc                                                                   <  2:8336738/101‑1 (MQ=255)
                              gCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCAc                                                                   <  1:8336640/101‑1 (MQ=255)
                               ccGAAACAGGCAGCAGTCTGCGGTGACGGCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATaa                                                                                            <  1:8336669/75‑1 (MQ=255)
                                cGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCAc                                                                   >  1:8337158/1‑99 (MQ=255)
                                 gAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGcc                                                                <  1:8336546/101‑1 (MQ=255)
                                 gAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGcc                                                                >  2:8337160/1‑101 (MQ=255)
                                 gAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGcc                                                                >  1:8337161/1‑101 (MQ=255)
                                 gAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGcc                                                                <  2:8336462/101‑1 (MQ=255)
                                  aaaCAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCg                                                               >  1:8337163/1‑101 (MQ=255)
                                     cAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTa                                                            <  1:8336485/101‑1 (MQ=255)
                                       ggCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGa                                                          >  2:8337167/1‑101 (MQ=255)
                                        gcagcaGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGAt                                                         >  2:8337169/1‑101 (MQ=255)
                                         cagcagTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATg                                                        <  1:8336623/101‑1 (MQ=255)
                                              gTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATg                                                   >  1:8337172/1‑101 (MQ=255)
                                               tCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATgg                                                  >  2:8337174/1‑101 (MQ=255)
                                                 tGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCg                                                >  2:8337178/1‑101 (MQ=255)
                                                   cGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGAc                                              >  2:8337180/1‑101 (MQ=255)
                                                      tGACGTCCACCTGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGc                                           >  1:8337182/1‑101 (MQ=255)
                                                      tGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGc                                           >  2:8337185/1‑101 (MQ=255)
                                                      tGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGc                                           <  2:8336675/101‑1 (MQ=255)
                                                          gTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGt                                       >  2:8337189/1‑101 (MQ=255)
                                                          gTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGt                                       >  1:8337191/1‑101 (MQ=255)
                                                          gTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGt                                       >  1:8337188/1‑101 (MQ=255)
                                                           tCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTa                                      >  1:8337194/1‑101 (MQ=255)
                                                            ccACCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTaa                                     >  2:8337195/1‑101 (MQ=255)
                                                              aCCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACa                                   >  1:8337200/1‑101 (MQ=255)
                                                              aCCGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACa                                   >  2:8337201/1‑101 (MQ=255)
                                                               ccGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACaa                                  >  1:8337204/1‑101 (MQ=255)
                                                                cGGTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACaaa                                 >  2:8337205/1‑101 (MQ=255)
                                                                 ggTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTa                                                            >  2:8337208/1‑73 (MQ=255)
                                                                 ggTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGc                                           <  2:8336776/90‑1 (MQ=255)
                                                                 ggTGATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAg                                <  1:8336787/101‑1 (MQ=255)
                                                                    gATGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGccc                             >  2:8337210/1‑101 (MQ=255)
                                                                     aTGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCa                            <  1:8336723/101‑1 (MQ=255)
                                                                      tGCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAg                           <  2:8336703/101‑1 (MQ=255)
                                                                       gCTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGa                          <  1:8336879/101‑1 (MQ=255)
                                                                        cTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGaa                         >  2:8337214/1‑101 (MQ=255)
                                                                        cTACGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGaa                         >  1:8337213/1‑101 (MQ=255)
                                                                          aCGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATc                       <  2:8336558/101‑1 (MQ=255)
                                                                           cGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCg                      >  2:8337219/1‑101 (MQ=255)
                                                                           cGCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCg                      >  1:8337220/1‑101 (MQ=255)
                                                                            gCATCATATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCgg                     >  2:8337222/1‑101 (MQ=255)
                                                                             catcatATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCggg                    >  2:8337225/1‑101 (MQ=255)
                                                                              atcatATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGt                   >  2:8337228/1‑101 (MQ=255)
                                                                               tcatATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCggga                                                            >  1:8337230/1‑57 (MQ=255)
                                                                                catATGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTAt                 >  2:8337232/1‑101 (MQ=255)
                                                                                  tatGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgg               >  1:8337236/1‑101 (MQ=255)
                                                                                  tatGCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATgg               <  1:8336581/101‑1 (MQ=255)
                                                                                     gCGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAt                        <  2:8336614/89‑1 (MQ=255)
                                                                                      cGGCACCACGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCgg                     >  2:8337238/1‑91 (MQ=255)
                                                                                       ggCACCACGAGGCAAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCt                                          <  1:8336607/69‑1 (MQ=255)
                                                                                          accacGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGa       >  1:8337240/1‑101 (MQ=255)
                                                                                           ccacGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGa       <  2:8336525/100‑1 (MQ=255)
                                                                                            cacGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAAt                        <  1:8336309/82‑1 (MQ=255)
                                                                                            cacGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCgg                     <  2:8336971/85‑1 (MQ=255)
                                                                                            cacGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGaca     >  1:8337242/1‑101 (MQ=255)
                                                                                             acGAGGCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGacac    >  2:8337244/1‑101 (MQ=255)
                                                                                               gaggCCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGacacac  >  1:8337248/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTACGCATCATATGCGGCACCACGAGACCAATAAAGCTAATCGGGCCGCAGGCGGCCACGCCGGTAGATGTCATGGCGACAGCTAGTAACAAAGCCCAGAATCGGGTATGGGGCACCGACACAC  >  NC_000913/4511567‑4511762

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: