Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,496 T→G 6.9% S232R (AGT→CGT)  fepD ← ferric enterobactin ABC transporter membrane subunit FebD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913621,4960TG 6.9% 212.7 / 7.9 73S232R (AGT→CGT) fepDferric enterobactin ABC transporter membrane subunit FebD
Reads supporting (aligned to +/- strand):  ref base T (30/38);  new base G (3/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

ACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCA  >  NC_000913/621349‑621643
                                                                                                                                                   |                                                                                                                                                   
aCCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGc                                                                                                                                                   <  2:1523911/150‑1 (MQ=255)
 ccAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGcc                                                                                                                                                  <  1:236319/150‑1 (MQ=255)
   agagatgatCGGCCCCAACCAGACAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGAACTACCACAAAGCAAGGTAATCGCCAGCAGACCAATCAACTGAGTGCGCGCCACGCGACTGACCa                                                                                                                                                >  2:1421012/4‑150 (MQ=255)
    gCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAg                                                                                                                                               >  2:836921/1‑150 (MQ=255)
          atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGt                                                                                                                                         >  2:1268080/1‑150 (MQ=255)
               cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                                                                                                                    >  2:1149511/1‑150 (MQ=255)
               cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                                                                                                                    <  1:1605561/150‑1 (MQ=255)
                acccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtg                                                                                                                                   <  1:997626/150‑3 (MQ=255)
                 cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGgtg                                                                                                                                  <  1:293987/150‑1 (MQ=255)
                   caccaGCCAACGCGCCATTTGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCTTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCACAGCGCCGTCGCGGTGTc                                                                                                                                <  2:1345324/150‑1 (MQ=255)
                   caccaGCCAACGCGCCATATGGGGCATGATCATGCCAATAAAGGCAATCTGGCCAACTATTGCCGTGGGACTACCACAAAGGACGGTAATCGCCAGCAGACGAATCATCTGTGTGCTCGCCAGGCGACTGCCCAGCGCCGTCGCGGTGTc                                                                                                                                <  1:1552353/150‑1 (MQ=255)
                   caccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACTGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTc                                                                                                                                >  2:875096/1‑150 (MQ=255)
                   caccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGCCCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTc                                                                                                                                >  1:1332229/1‑150 (MQ=255)
                   caccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTc                                                                                                                                >  1:1535184/1‑150 (MQ=255)
                    accaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCg                                                                                                                               <  1:1252266/150‑1 (MQ=255)
                     ccaGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACGGCGCTGGGCAGTcgcgtggcgcg                                                                                                                              <  2:1228865/150‑8 (MQ=255)
                      caGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCt                                                                                                                             >  1:1395820/1‑150 (MQ=255)
                       aGCCAAGGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTg                                                                                                                            <  2:1330663/150‑1 (MQ=255)
                          cAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCg                                                                                                                         >  2:481160/1‑150 (MQ=255)
                          cAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCACTGCCg                                                                                                                         >  1:361555/1‑150 (MQ=255)
                              gcgcCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGc                                                                                                                     <  2:492650/150‑1 (MQ=255)
                                    tatGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAAc                                                                                                               >  2:1237175/1‑150 (MQ=255)
                                         ggCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTc                                                                                                          <  2:1082575/150‑1 (MQ=255)
                                           catcatCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAg                                                                                                        <  1:875096/150‑1 (MQ=255)
                                              catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                     <  1:1268080/150‑1 (MQ=255)
                                              catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                     <  1:41259/150‑1 (MQ=255)
                                              catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                     >  1:737708/1‑150 (MQ=255)
                                               atcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcg                                                                                                    <  2:1494403/150‑1 (MQ=255)
                                                tcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgcgc                                                                                                   >  2:1477937/1‑150 (MQ=255)
                                                     ccAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTc                                                                                              <  2:1108265/150‑1 (MQ=255)
                                                      cAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCa                                                                                             >  1:1097387/1‑150 (MQ=255)
                                                       aaTAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAg                                                                                            >  2:77900/1‑150 (MQ=255)
                                                         tAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTa                                                                                          <  2:1535184/150‑1 (MQ=255)
                                                          aaaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTaa                                                                                         >  2:1249828/1‑150 (MQ=255)
                                                           aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                        <  1:845430/150‑1 (MQ=255)
                                                             ggCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAg                                                                                      <  1:1214557/150‑1 (MQ=255)
                                                             ggCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAg                                                                                      <  1:695287/150‑1 (MQ=255)
                                                                aaTCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAgcgc                                                                                   <  2:470237/150‑1 (MQ=255)
                                                                               tGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGAt                                                                    <  1:836921/150‑1 (MQ=255)
                                                                                gCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                   >  1:631867/1‑150 (MQ=255)
                                                                                gCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATc                                                                   <  2:631867/150‑1 (MQ=255)
                                                                                 ccGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCa                                                                  >  1:1468661/1‑150 (MQ=255)
                                                                                   gTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcgc                                                                                                                                      <  1:1289363/80‑1 (MQ=255)
                                                                                   gTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcgc                                                                                                                                      >  2:1289363/1‑80 (MQ=255)
                                                                                   gTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGc                                                                <  1:1446089/150‑1 (MQ=255)
                                                                                          tACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGa                                                         <  2:1214085/150‑1 (MQ=255)
                                                                                            ccACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATc                                                       <  2:1274830/150‑1 (MQ=255)
                                                                                            ccACAAAGCACGGTAATCGCCAACAGACCAATCAACTGTGTGCGCGCCACCCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGCCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCCGCGATCAGCACCGGGATc                                                       >  2:1490520/1‑150 (MQ=255)
                                                                                              acaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAg                                                     <  1:593265/150‑1 (MQ=255)
                                                                                               caAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGc                                                    <  1:1506939/150‑1 (MQ=255)
                                                                                                  aGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGcacc                                                 >  2:1242257/1‑150 (MQ=255)
                                                                                                         tAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTa                                          >  2:50031/1‑150 (MQ=255)
                                                                                                         tAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTa                                          <  1:481160/150‑1 (MQ=255)
                                                                                                          aaTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTaa                                         >  1:1022591/1‑150 (MQ=255)
                                                                                                          aaTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTaa                                         >  1:893782/1‑150 (MQ=255)
                                                                                                          aaTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCACCGCCGTCGCGGTGTCGCTGCCCAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCAACACttaaaa                                         >  1:908078/1‑146 (MQ=255)
                                                                                                           aTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAg                                        <  1:77900/150‑1 (MQ=255)
                                                                                                                       ccAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTg                            >  1:1380642/1‑150 (MQ=255)
                                                                                                                           tCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGaa                        <  1:1149511/150‑1 (MQ=255)
                                                                                                                            cAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGGCGCGGTGGCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAat                       <  2:391871/150‑1 (MQ=255)
                                                                                                                              aCTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAatat                     >  2:196130/1‑150 (MQ=255)
                                                                                                                                   gtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCTCGCGACTCAGTAATAGAGGAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGc                <  1:1490520/150‑1 (MQ=255)
                                                                                                                                   gtgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGc                <  2:1332229/150‑1 (MQ=255)
                                                                                                                                    tgCGCGCCTCGCGACTGCGCAGCGCCGTCGCGGTGTCGCTGCCGAGGCGCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCg               <  2:267773/150‑1 (MQ=255)
                                                                                                                                    tgCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCTCAGTTGCTCCGGCGATCAGCACCGTGATCAGCACCACTTTTAAGGTATGTAGATTGCTAATATCCAGCg               <  2:1097387/150‑1 (MQ=255)
                                                                                                                                      cgcgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGaa             >  1:445218/1‑150 (MQ=255)
                                                                                                                                        cgcCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAAcc           >  1:565360/1‑150 (MQ=255)
                                                                                                                                          caccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCgg         >  1:1324301/4‑150 (MQ=255)
                                                                                                                                           cACGCGACTGCCCAGCGCCGTCGCGGTGTCACTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGc        <  2:361555/150‑1 (MQ=255)
                                                                                                                                            ccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGACCCGGCt       >  2:858788/2‑150 (MQ=255)
                                                                                                                                            ccgcgACGGCGCTGGGCAGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCt       >  1:1169237/2‑150 (MQ=255)
                                                                                                                                                gACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGcc   <  1:1452037/150‑1 (MQ=255)
                                                                                                                                                 aCTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCCCCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCa  >  1:1134686/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCA  >  NC_000913/621349‑621643

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