Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,916 G→A 5.5% intergenic (‑191/+19) 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,9160GA 5.5% 290.1 / 4.2 90intergenic (‑191/+19)uof/fldARyhB‑regulated fur leader peptide/flavodoxin 1
Reads supporting (aligned to +/- strand):  ref base G (43/42);  new base A (3/2);  total (46/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 = 8.28e-01

CAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCAT  >  NC_000913/710767‑711062
                                                                                                                                                     |                                                                                                                                                  
cAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACGACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATg                                                                                                                                                    >  1:890575/1‑150 (MQ=255)
cAAATATTGTTAATAAAAACGAGGCACAGCTTACATTTACAACGGCAAGGCCGCTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCTATATCTATTAATACGTAGTAGACTTg                                                                                                                                                    <  2:775100/150‑1 (MQ=255)
      ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTa                                                                                                                                              <  1:1407961/150‑1 (MQ=255)
      ttGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTa                                                                                                                                              >  2:1488100/1‑150 (MQ=255)
        gTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcg                                                                                                                                            >  1:1325870/1‑150 (MQ=255)
          taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                          <  2:50521/150‑1 (MQ=255)
          taataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTAcgcc                                                                                                                                          <  1:798148/150‑1 (MQ=255)
           aataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAATAGTTACTAATTTGTAGCATTTATTTTGATTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgccg                                                                                                                                         <  2:793618/150‑1 (MQ=255)
            ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGACGTAGACATGAGTCTAcgccgc                                                                                                                                        <  1:1021921/150‑1 (MQ=255)
            ataaAAACGTGGCACAGCTTACATTTACAACGGCAAGGACACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGTTTGGCATTATCTATTAATACGGCGTAGACTCATGTCTAcgccgt                                                                                                                                        <  1:1073172/150‑2 (MQ=255)
             taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                       <  1:699250/150‑1 (MQ=255)
             taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGca                                                                                                                                       <  1:1175078/150‑1 (MQ=255)
             taaAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGACGca                                                                                                                                       >  1:1133351/1‑150 (MQ=255)
               aaaaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatc                                                                                                                                     <  1:377869/150‑1 (MQ=255)
                 aaCGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcac                                                                                                                                   <  1:1056692/150‑1 (MQ=255)
                     tGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGcatcacatca                                                                                                                               <  2:915180/150‑1 (MQ=255)
                      ggCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAg                                                                                                                              >  1:471272/1‑150 (MQ=255)
                       gCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAgg                                                                                                                             >  1:1521882/1‑150 (MQ=255)
                        cacaGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGc                                                                                                                            >  1:996746/1‑150 (MQ=255)
                           aGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCAtt                                                                                                                         >  2:584543/1‑150 (MQ=255)
                            gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                                                                                                                        <  1:1002949/150‑1 (MQ=255)
                            gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTg                                                                                                                        <  2:893818/150‑1 (MQ=255)
                            gCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCCACGCCGCATCACATCAGGCATTg                                                                                                                        >  1:807830/1‑150 (MQ=255)
                            gCTTACATTTACAACGGCAAGGCCACTTCAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCATGCATTg                                                                                                                        <  2:994647/150‑1 (MQ=255)
                             cTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTga                                                                                                                       >  2:523866/1‑150 (MQ=255)
                              ttACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTgag                                                                                                                      <  2:804898/150‑1 (MQ=255)
                                 cATTTACAACGGCAAGTCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                                                                                                                   <  2:748096/150‑1 (MQ=255)
                                 cATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAAt                                                                                                                   >  2:1004532/1‑150 (MQ=255)
                                    ttACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACGACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtc                                                                                                                <  2:890575/150‑1 (MQ=255)
                                     tACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcg                                                                                                               >  2:1082216/1‑150 (MQ=255)
                                       cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtc                                                                                                             <  2:572210/150‑1 (MQ=255)
                                       cAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtc                                                                                                             >  2:1414194/1‑150 (MQ=255)
                                        aaCGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTtcgtcg                                                                                                            >  2:1366983/1‑150 (MQ=255)
                                          cGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag                                                                                                          >  2:1044324/1‑150 (MQ=255)
                                                 gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                   >  2:1109881/1‑150 (MQ=255)
                                                 gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                   <  1:1488100/150‑1 (MQ=255)
                                                 gCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCCACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAAc                                                                                                   <  2:807830/150‑1 (MQ=255)
                                                        taaCACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTATGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAg                                                                                            <  1:304586/149‑1 (MQ=255)
                                                        gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAg                                                                                            <  1:1109881/150‑1 (MQ=255)
                                                        gAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAg                                                                                            <  2:471272/150‑1 (MQ=255)
                                                          acacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGaa                                                                                          <  1:1004532/150‑1 (MQ=255)
                                                           cacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGaaa                                                                                         <  1:934526/150‑1 (MQ=255)
                                                           cacCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGaaa                                                                                         >  1:1058952/1‑150 (MQ=255)
                                                             ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag                                                                                                          >  1:782105/1‑131 (MQ=255)
                                                             ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCgag                                                                                                          <  2:782105/131‑1 (MQ=255)
                                                             ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGatcgg                                                                                       >  1:266268/1‑146 (MQ=255)
                                                             ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATc                                                                                       <  1:1156610/150‑1 (MQ=255)
                                                             ccGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATc                                                                                       >  2:218264/1‑150 (MQ=255)
                                                                 gACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGtt                                                                                   <  1:945487/150‑1 (MQ=255)
                                                                  ctgttGTACCTGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGttt                                                                                  >  2:872844/4‑150 (MQ=255)
                                                                     ttGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAAc                                                                               <  1:853171/150‑1 (MQ=255)
                                                                        tACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCa                                                                            >  2:730033/1‑150 (MQ=255)
                                                                          cAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCAtt                                                                          >  2:130903/1‑150 (MQ=255)
                                                                                 aacaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTac                                                                   >  2:1186965/1‑150 (MQ=255)
                                                                                  acaGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTaca                                                                  >  1:647837/1‑150 (MQ=255)
                                                                                    aGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACg                                                                >  2:497973/1‑150 (MQ=255)
                                                                                     gCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAAACTGTTTAACCCATTTTACTACACGt                                                               >  2:618045/1‑150 (MQ=255)
                                                                                          aGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGc                                                          <  2:1521882/150‑1 (MQ=255)
                                                                                          aGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGc                                                          <  2:1537630/150‑1 (MQ=255)
                                                                                            ttACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCgg                                                        <  1:916114/150‑1 (MQ=255)
                                                                                                 aaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATGTGTTTAACCCATTTTTGTACAGGTTCAGCGGTAAGt                                                   >  2:1044382/1‑150 (MQ=255)
                                                                                                 aaTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGt                                                   >  2:340863/1‑150 (MQ=255)
                                                                                                  aTTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGtt                                                  >  2:752623/1‑150 (MQ=255)
                                                                                                          cAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGa                                          >  1:764003/1‑150 (MQ=255)
                                                                                                           aaTTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCGTTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGAc                                         <  1:1081933/150‑1 (MQ=255)
                                                                                                           aaTTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGAc                                         >  2:71185/1‑150 (MQ=255)
                                                                                                            attattTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg                                        <  1:1082216/150‑1 (MQ=255)
                                                                                                            attattTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACg                                        <  1:1408401/150‑1 (MQ=255)
                                                                                                             ttattTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACgg                                       <  1:1075506/150‑1 (MQ=255)
                                                                                                             aaagtttGATTCGCATAATCTACTAATACGGCGCAGAGATGAGGATACGCCTCATCACATCACGCATTCAGAATTGCGGCGATATGCAACTCTTCAGAAATATGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCTGCTGACgg                                       <  2:1175798/146‑1 (MQ=255)
                                                                                                                      ttGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgtcg                              <  1:730033/150‑1 (MQ=255)
                                                                                                                      ttGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTATACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTtcgtcg                              >  1:876813/1‑150 (MQ=255)
                                                                                                                         gCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATa                           <  1:112699/150‑1 (MQ=255)
                                                                                                                          cATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAg                          <  2:372167/150‑1 (MQ=255)
                                                                                                                            ttATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGcc                        >  1:1009671/1‑150 (MQ=255)
                                                                                                                            ttATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGcc                        >  1:1014688/1‑150 (MQ=255)
                                                                                                                               tCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGa                     >  2:27182/1‑150 (MQ=255)
                                                                                                                               tCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGa                     >  2:1472214/1‑150 (MQ=255)
                                                                                                                                cTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGAc                    >  2:1453630/1‑150 (MQ=255)
                                                                                                                                 tATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGAcc                   <  1:523866/150‑1 (MQ=255)
                                                                                                                                  aTTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCg                  >  1:604650/1‑150 (MQ=255)
                                                                                                                                  aTTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCg                  >  2:122724/1‑150 (MQ=255)
                                                                                                                                   ttAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGa                 >  2:309900/1‑150 (MQ=255)
                                                                                                                                     aaTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACa               >  1:698703/1‑150 (MQ=255)
                                                                                                                                     aaTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACa               <  1:872844/150‑1 (MQ=255)
                                                                                                                                      aTACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACaa              <  1:1366983/150‑1 (MQ=255)
                                                                                                                                       tACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACaaa             >  1:1234945/1‑150 (MQ=255)
                                                                                                                                        gcGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAg            >  1:1046875/2‑150 (MQ=255)
                                                                                                                                         cGGCGTAGACTCATGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGt           >  2:277934/1‑150 (MQ=255)
                                                                                                                                          ggCGTAGACTAATGTCTAAGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAAACAATTTTTCTAAACGTTCAGCGGTCAGTTCCGGCTGACGGTATTCGTCGATAGCCAGACCGACAAAGtt          >  2:354115/1‑149 (MQ=255)
                                                                                                                                           gCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTgg         >  1:1020075/1‑150 (MQ=255)
                                                                                                                                                  cATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCAt  <  1:145875/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAATATTGTTAATAAAAACGTGGCACAGCTTACATTTACAACGGCAAGGCCACTTGAACACCGGGACATTGTACAGGTACAACAGCAAAAGTTACAAATTTGTAGCAATTATTTTGATTGGCATTATCTATTAATACGGCGTAGACATGAGTCTACGCCGCATCACATCAGGCATTGAGAATTTCGTCGAGATGCAACTCTTCAGAAATCTGTTTAACCCATTTTTCTACACGTTCAGCGGTCAGTTCCGGCTGACGGTCTTCGTCGATAGCCAGACCGACAAAGTGGTCGTCAT  >  NC_000913/710767‑711062

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: