Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 611,255 C→T A414T (GCT→ACT)  fepA ← ferrienterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913611,2550CT99.0% 381.5 / ‑4.5 98A414T (GCT→ACT) fepAferrienterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (53/44);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCAT  >  NC_000913/611108‑611401
                                                                                                                                                   |                                                                                                                                                  
tATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAgtgc                                                                                                                                                  <  1:789451/150‑1 (MQ=255)
    gCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCAcc                                                                                                                                              <  1:535510/150‑1 (MQ=255)
     ccGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCgg                                                                                                                                            <  1:199398/151‑1 (MQ=255)
        aCAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTAt                                                                                                                                         >  1:525370/1‑151 (MQ=255)
         cAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTAtt                                                                                                                                        <  1:1074460/151‑1 (MQ=255)
          aaTACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTc                                                                                                                                       <  1:958439/151‑1 (MQ=255)
           aTACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCg                                                                                                                                      >  2:188157/1‑151 (MQ=255)
                  gatgatCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGt                                                                                                                               >  1:714963/1‑151 (MQ=255)
                   atgatCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTc                                                                                                                              <  1:90398/151‑1 (MQ=255)
                     gatCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAg                                                                                                                            <  2:1408198/151‑1 (MQ=255)
                      atCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGt                                                                                                                           <  1:258474/151‑1 (MQ=255)
                       tCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTg                                                                                                                          >  1:623085/1‑151 (MQ=255)
                         gAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGcc                                                                                                                        >  1:312934/1‑151 (MQ=255)
                           aaCGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTg                                                                                                                      <  2:835492/151‑1 (MQ=255)
                              gCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGgt                                                                                                                   <  2:532096/151‑1 (MQ=255)
                                aGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGgtgt                                                                                                                 >  2:405499/1‑151 (MQ=255)
                                 gCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGtt                                                                                                                >  1:1176997/1‑151 (MQ=255)
                                  ccgcGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTg                                                                                                               <  1:361536/148‑1 (MQ=255)
                                     cGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGaa                                                                                                            <  1:925386/151‑1 (MQ=255)
                                     cGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGa                                                                                                             >  1:1411163/1‑150 (MQ=255)
                                     cGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGa                                                                                                             >  2:918798/1‑150 (MQ=255)
                                        cGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACtt                                                                                                         >  2:864903/1‑151 (MQ=255)
                                        cGTTACGATGGTGCTGCCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTGGAATACGGGCTACGTTCGGTGTTACTCACGGCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTg                                                                                                               <  2:1174311/145‑1 (MQ=255)
                                            aCGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAg                                                                                                      <  2:1053992/150‑1 (MQ=255)
                                              gATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTcc                                                                                                   >  2:313181/1‑151 (MQ=255)
                                              gATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTcc                                                                                                   <  1:1222783/151‑1 (MQ=255)
                                               aTGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCt                                                                                                  >  1:820666/1‑151 (MQ=255)
                                                 ggTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCtt                                                                                                 <  1:1415250/150‑1 (MQ=255)
                                                      tGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATcc                                                                                            >  1:1425065/1‑150 (MQ=255)
                                                      tGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCg                                                                                           >  1:486539/1‑151 (MQ=255)
                                                       gtcagtcaGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGt                                                                                          <  2:919566/151‑1 (MQ=255)
                                                          agtcagCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGc                                                                                       >  1:101377/1‑151 (MQ=255)
                                                          agtcagCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGc                                                                                       >  1:1100351/1‑151 (MQ=255)
                                                           gtcagCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCt                                                                                      >  2:473932/1‑151 (MQ=255)
                                                            tcagCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCt                                                                                      >  2:456314/1‑150 (MQ=255)
                                                              agCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGAt                                                                                   <  1:918885/151‑1 (MQ=255)
                                                                   cATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCAc                                                                              >  2:620884/1‑151 (MQ=255)
                                                                     tgttgtTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTc                                                                            >  2:573746/1‑151 (MQ=255)
                                                                      gttgttTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTc                                                                            >  1:1158254/1‑150 (MQ=255)
                                                                        tgttTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCg                                                                          <  1:1359729/150‑1 (MQ=255)
                                                                         gttTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGt                                                                         >  1:1361030/1‑150 (MQ=255)
                                                                             tCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTAc                                                                       >  2:184273/1‑148 (MQ=255)
                                                                              cGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTAccc                                                                     >  2:1419887/1‑149 (MQ=255)
                                                                              cGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcag                                                                   <  1:188181/151‑1 (MQ=255)
                                                                              cGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCca                                                                    <  2:687118/150‑1 (MQ=255)
                                                                                 cAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgt                                                                >  2:437207/1‑151 (MQ=255)
                                                                                 cAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgt                                                                >  2:1195280/1‑151 (MQ=255)
                                                                                    aCAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgtc                                                               >  1:247852/1‑149 (MQ=255)
                                                                                     cAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgtcagc                                                            >  1:182895/1‑151 (MQ=255)
                                                                                     cAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgtcagc                                                            <  2:1057910/151‑1 (MQ=255)
                                                                                       gCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgtcagcgt                                                          >  1:1419516/1‑151 (MQ=255)
                                                                                       gCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCcagcgtcagcgt                                                          >  2:426702/1‑151 (MQ=255)
                                                                                        cGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCagcgtcagcgtc                                                         <  2:1137403/151‑1 (MQ=255)
                                                                                        cGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCagcgtcagcgtc                                                         >  1:1156352/1‑151 (MQ=255)
                                                                                         gagaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCt                                                        <  1:1055388/151‑1 (MQ=255)
                                                                                          agaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTg                                                       <  1:864991/151‑1 (MQ=255)
                                                                                           gaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTgg                                                      >  1:481326/1‑151 (MQ=255)
                                                                                           gaAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTgg                                                      >  2:1455092/1‑151 (MQ=255)
                                                                                             aaaTTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGtt                                                    >  1:846553/1‑151 (MQ=255)
                                                                                              aaTTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGtt                                                    >  1:298884/1‑150 (MQ=255)
                                                                                                tttCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAAc                                                 >  1:1326739/1‑151 (MQ=255)
                                                                                                 ttCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACg                                                >  2:1009148/1‑151 (MQ=255)
                                                                                                     gCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGa                                               <  1:537892/148‑1 (MQ=255)
                                                                                                     gCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAg                                              <  2:717528/149‑1 (MQ=255)
                                                                                                        tttGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGaaa                                          <  1:789960/150‑1 (MQ=255)
                                                                                                        tttGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGa                                            <  2:312896/148‑1 (MQ=255)
                                                                                                         ttGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaat                                         <  1:1449106/150‑1 (MQ=255)
                                                                                                          tGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatc                                        >  1:139983/1‑150 (MQ=255)
                                                                                                           gAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatca                                       >  1:1399346/1‑150 (MQ=255)
                                                                                                           gAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatca                                       >  2:246025/1‑150 (MQ=255)
                                                                                                            aaTACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatcaat                                     >  1:286821/1‑151 (MQ=255)
                                                                                                            aaTACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAaatcaat                                     >  1:380952/1‑151 (MQ=255)
                                                                                                               aCGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCgg                                  >  2:1439563/1‑151 (MQ=255)
                                                                                                               aCGGGCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCgg                                  >  2:239804/1‑151 (MQ=255)
                                                                                                                 gggCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGc                                 >  2:920815/1‑150 (MQ=255)
                                                                                                                 gggCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCa                                <  2:239377/151‑1 (MQ=255)
                                                                                                                  ggCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAg                               <  2:1043037/151‑1 (MQ=255)
                                                                                                                   gCTACGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAgg                              <  2:1344521/151‑1 (MQ=255)
                                                                                                                      aCGGTCGGTGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGt                             <  1:681072/149‑1 (MQ=255)
                                                                                                                       cggtcggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTa                           >  2:1327837/1‑150 (MQ=255)
                                                                                                                        ggtcggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTaa                          >  1:700657/1‑150 (MQ=255)
                                                                                                                        ggtcggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAAc                         <  1:958702/151‑1 (MQ=255)
                                                                                                                        ggtcggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAAc                         >  2:591607/1‑151 (MQ=255)
                                                                                                                           cggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTc                      <  2:12633/151‑1 (MQ=255)
                                                                                                                           cggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTc                      <  1:811954/151‑1 (MQ=255)
                                                                                                                           cggtGGTACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTc                      <  1:230005/151‑1 (MQ=255)
                                                                                                                              tggtACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCt                   >  1:177436/1‑151 (MQ=255)
                                                                                                                               ggtACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCt                   <  1:1378490/150‑1 (MQ=255)
                                                                                                                               ggtACTCACGCCAGCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtg                  <  2:1005146/151‑1 (MQ=255)
                                                                                                                                gtACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgt                 <  2:1020054/151‑1 (MQ=255)
                                                                                                                                gtACTCACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCtgt                 <  2:527076/151‑1 (MQ=255)
                                                                                                                                     cACGCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTgca              <  1:313218/149‑1 (MQ=255)
                                                                                                                                        gCCATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATc         >  1:1205297/1‑151 (MQ=255)
                                                                                                                                         ccATCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCa        <  1:898910/151‑1 (MQ=255)
                                                                                                                                           aTCAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACg      >  2:708764/1‑151 (MQ=255)
                                                                                                                                             cAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACg      >  2:647397/1‑149 (MQ=255)
                                                                                                                                             cAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGtc    >  2:930598/1‑151 (MQ=255)
                                                                                                                                             cAATAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGtc    >  2:954227/1‑151 (MQ=255)
                                                                                                                                               aTAGTGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGtcat  <  2:1255441/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATTGCCGACAATACTGTGATGATCGAAACGCAGCCCCGGCGTTACGATGGTGCTGTCAGTCAGCTCCATGTTGTTTTCGGCAAACAGCGAGAAAATTTCTGCTTTTGAATACGGGCTACGGTCGGTGGTACTCACGCCATCAATAGCGCCACCGGTATTCGTTCCGGTCAGTGCCTGGGTGTTGGAACTTAAGTCCTTCATCCGTTGCTGATTCCACTCCGTACCCAGCGTCAGCGTCTGGTTAACGAGGAAATCAATCGGCAGGTTAACTTCGCTGTGCAGCATCACGTCAT  >  NC_000913/611108‑611401

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: