Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,917 A→T 5.5% intergenic (‑192/+18) uof ← / ← fldA RyhB‑regulated fur leader peptide/flavodoxin 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,9170AT 5.5% 296.3 / 5.6 91intergenic (‑192/+18)uof/fldARyhB‑regulated fur leader peptide/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base T (3/2);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTC  >  NC_000913/710769‑711060
                                                                                                                                                    |                                                                                                                                               
aaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                                                                                                                <  2:1286409/150‑1 (MQ=255)
aaTATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAg                                                                                                                                                <  1:740442/150‑1 (MQ=255)
    ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTa                                                                                                                                            >  1:1129279/1‑150 (MQ=255)
       ttAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgc                                                                                                                                         >  1:1221715/1‑150 (MQ=255)
       ttAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgc                                                                                                                                         >  1:783330/1‑150 (MQ=255)
        taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                        >  1:1613231/1‑150 (MQ=255)
        taataaAAACGTCGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgcc                                                                                                                                        <  1:1490919/150‑1 (MQ=255)
         aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgccg                                                                                                                                       <  1:604027/150‑1 (MQ=255)
         aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccg                                                                                                                                       >  2:1033086/1‑150 (MQ=255)
          ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccgc                                                                                                                                      >  2:788418/1‑150 (MQ=255)
          ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgccgc                                                                                                                                      <  1:871909/150‑1 (MQ=255)
           taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                     <  1:659820/150‑1 (MQ=255)
           taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGACATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                     <  1:1276086/150‑1 (MQ=255)
                aCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGAATGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcaca                                                                                                                                <  1:1136551/150‑1 (MQ=255)
                      cacaGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGc                                                                                                                          >  1:1570582/1‑150 (MQ=255)
                       acaGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCa                                                                                                                         <  1:991530/150‑1 (MQ=255)
                       acaGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCa                                                                                                                         <  2:1647996/150‑1 (MQ=255)
                           cTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                                                                                                                     <  2:1178756/150‑1 (MQ=255)
                            ttACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgag                                                                                                                    <  2:1602528/150‑1 (MQ=255)
                             tACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgaga                                                                                                                   <  2:1485376/150‑1 (MQ=255)
                              aCATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGaa                                                                                                                  <  2:90387/150‑1 (MQ=255)
                               cATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                                                                                                                 >  1:1442273/1‑150 (MQ=255)
                                aTTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAtt                                                                                                                >  2:1162891/1‑150 (MQ=255)
                                  ttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtc                                                                                                              <  1:1632854/150‑1 (MQ=255)
                                   tACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                                                                                                             >  2:1252655/1‑150 (MQ=255)
                                     cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtc                                                                                                           <  2:1422853/150‑1 (MQ=255)
                                      aaCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtcg                                                                                                          <  2:1109168/150‑1 (MQ=255)
                                       aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCga                                                                                                         <  2:567907/150‑1 (MQ=255)
                                       aCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCga                                                                                                         >  1:743337/1‑150 (MQ=255)
                                        cGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCACTGAGAATTTCGTCgag                                                                                                        <  1:1246425/150‑1 (MQ=255)
                                            aaGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGc                                                                                                    <  2:694090/150‑1 (MQ=255)
                                             aGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCa                                                                                                   <  1:271124/150‑1 (MQ=255)
                                               gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                 <  2:1523877/150‑1 (MQ=255)
                                                      gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAg                                                                                          <  2:1434768/150‑1 (MQ=255)
                                                           ccGGGACATTGTTCAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATc                                                                                     <  2:446775/150‑1 (MQ=255)
                                                           ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATc                                                                                     >  1:1533641/1‑150 (MQ=255)
                                                           ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATc                                                                                     >  1:1291607/1‑150 (MQ=255)
                                                              ggACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGt                                                                                  >  1:1225119/1‑150 (MQ=255)
                                                                aCATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGttt                                                                                <  1:788418/150‑1 (MQ=255)
                                                                  aTTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTaa                                                                              >  2:1284142/1‑150 (MQ=255)
                                                                      tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGAATTGAGAATTTCGTCGAGATGCAACTATTCATAAATCTGATTaaaaaa                                                                          >  1:723114/1‑146 (MQ=255)
                                                                        cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAtt                                                                        >  1:1554676/1‑150 (MQ=255)
                                                                        cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAtt                                                                        >  2:1576491/1‑150 (MQ=255)
                                                                        cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAtt                                                                        >  2:325687/1‑150 (MQ=255)
                                                                         aGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAttt                                                                       <  2:466684/150‑1 (MQ=255)
                                                                            tACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTc                                                                    <  1:256306/150‑1 (MQ=255)
                                                                              caacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTa                                                                  >  1:1616049/1‑150 (MQ=255)
                                                                               aacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTac                                                                 >  2:420853/1‑150 (MQ=255)
                                                                                acaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTaca                                                                <  1:1170807/150‑1 (MQ=255)
                                                                                 caGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTacac                                                               <  2:1593465/150‑1 (MQ=255)
                                                                                 caGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTacac                                                               >  2:365886/1‑150 (MQ=255)
                                                                                  aGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACg                                                              >  2:1084450/1‑150 (MQ=255)
                                                                                     aaaaGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTc                                                           >  2:495763/1‑150 (MQ=255)
                                                                                        aGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGc                                                        <  2:1637628/150‑1 (MQ=255)
                                                                                         cttACAAATTTGTAGCAATTAATTTTATAGCAATTATCTATTAATACGGCGTAGACAATATTCTACGACGAATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCg                                                       >  2:592057/2‑150 (MQ=255)
                                                                                          ttACAAATTTGTATCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCgg                                                      <  1:1010927/150‑1 (MQ=255)
                                                                                          ttACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCgg                                                      <  1:833723/150‑1 (MQ=255)
                                                                                          ttACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCgg                                                      <  1:35511/150‑1 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCgggc                                                    <  2:573910/150‑3 (MQ=255)
                                                                                            aCAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTc                                                    <  1:1539078/150‑1 (MQ=255)
                                                                                              aaaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAg                                                  <  2:1613231/150‑1 (MQ=255)
                                                                                                   tGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCg                                             <  2:743337/150‑1 (MQ=255)
                                                                                                    gTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATGACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGCTCAgg                                            >  1:1537934/1‑150 (MQ=255)
                                                                                                     tAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGAATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGc                                           <  1:494456/150‑1 (MQ=255)
                                                                                                      aGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCCACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCt                                          >  1:533089/1‑150 (MQ=255)
                                                                                                       gCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTg                                         >  2:1396919/1‑150 (MQ=255)
                                                                                                         acttattTTGATTGGCATTATCTATTATTACTGCGTAGACATGAGTCTACGCCGCATCACATCAGGCCTTGAGAATTTCGTCGAGATGCAACTCATCAGAAATCTGTTTAACCCATATTTCTACACGTTCAGCGGTCAGTTCCGGCTGAc                                       <  2:456913/148‑1 (MQ=255)
                                                                                                             attTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTc                                   >  2:1039921/1‑150 (MQ=255)
                                                                                                             attTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTc                                   >  1:79495/1‑150 (MQ=255)
                                                                                                              ttttGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCt                                  >  2:674715/1‑150 (MQ=255)
                                                                                                                ttGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtc                                >  1:1388945/1‑150 (MQ=255)
                                                                                                                  taatcgCATAATCTATTAAAACGGCGTAGACATGAGTCTAACACGCATCACATCAGGCAATTAGAAATTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgt                              >  2:54249/6‑150 (MQ=255)
                                                                                                                  gATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgt                              <  1:1162891/150‑1 (MQ=255)
                                                                                                                   aTTGGCATTATCTATTAATACGGCGTAGACATGAGTCCACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgtc                             <  2:533089/150‑1 (MQ=255)
                                                                                                                    ttGGCATTATATATTAAAACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGCGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgtcg                            <  2:1074082/150‑1 (MQ=255)
                                                                                                                        cATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAg                        <  1:273143/150‑1 (MQ=255)
                                                                                                                            aTCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAg                    <  1:530527/150‑1 (MQ=255)
                                                                                                                            aTCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAg                    >  2:3881/1‑150 (MQ=255)
                                                                                                                               tATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGAcc                 >  1:1488523/1‑150 (MQ=255)
                                                                                                                               tATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGAcc                 >  1:1493065/1‑150 (MQ=255)
                                                                                                                                aTTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTGTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCg                >  2:549970/1‑150 (MQ=255)
                                                                                                                                aTTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCg                >  2:1031586/1‑150 (MQ=255)
                                                                                                                                 ttAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGa               <  2:1442273/150‑1 (MQ=255)
                                                                                                                                   aaTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACa             <  2:1225119/150‑1 (MQ=255)
                                                                                                                                       cGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGt         >  2:1245254/1‑150 (MQ=255)
                                                                                                                                       cGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGt         <  1:1576491/150‑1 (MQ=255)
                                                                                                                                        ggCGTAGACTCATGTCTACGCCGCAACACATCAGGCATTGAGAATATCGACGAGATGCAACTCTTCAGAAATCACTTTAACCCATTTTTATACACCTTCAGCGGTCAGTTACGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTg        >  2:163623/1‑150 (MQ=255)
                                                                                                                                        ggCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTg        >  1:1067503/1‑150 (MQ=255)
                                                                                                                                          cGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgt      <  1:487205/150‑1 (MQ=255)
                                                                                                                                              gACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgtcgtc  >  2:476609/1‑150 (MQ=255)
                                                                                                                                              gACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGgtcgtc  >  2:1047500/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTC  >  NC_000913/710769‑711060

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