Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 621,501 A→T 6.6% L230Q (CTG→CAG)  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,5010AT 6.6% 210.1 / 3.3 70L230Q (CTG→CAG) fepDferric enterobactin ABC transporter membrane subunit FebD
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (3/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

AGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGC  >  NC_000913/621352‑621650
                                                                                                                                                     |                                                                                                                                                     
agagatgatCGGCCCCAACCAGACAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGAACTACCACAAAGCAAGGTAATCGCCAGCAGACCAATCAACTGAGTGCGCGCCACGCGACTGACCa                                                                                                                                                       >  2:1421012/4‑150 (MQ=255)
 gCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAg                                                                                                                                                      >  2:836921/1‑150 (MQ=255)
       atCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGt                                                                                                                                                >  2:1268080/1‑150 (MQ=255)
            cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                                                                                                                           <  1:1605561/150‑1 (MQ=255)
            cacccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCgg                                                                                                                                           >  2:1149511/1‑150 (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:737708/1‑150 (MQ=255)
                                           catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                            <  1:1268080/150‑1 (MQ=255)
                                           catcaGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAgcgc                                                                                                            <  1:41259/150‑1 (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                                                                                                                                             >  2:1289363/1‑80 (MQ=255)
                                                                                gTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTcgc                                                                                                                                             <  1:1289363/80‑1 (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)
                                                                                                                                                     aGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGc  >  2:111962/1‑150 (MQ=255)
                                                                                                                                                     aGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCAAGCGAACCGGCTTGCCAGAAACGc  >  2:117776/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGC  >  NC_000913/621352‑621650

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: ‑