Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,473 C→G 5.7% intergenic (+108/+36) galU → / ← hns glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,292,4730CG 5.7% 616.0 / 23.0 175intergenic (+108/+36)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base C (78/87);  new base G (6/4);  total (84/91)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 179 total aligned reads displayed.

TTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGT  >  NC_000913/1292325‑1292621
                                                                                                                                                     |                                                                                                                                                    
   aaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTgg                                                                                                                                                   <  3:376507/149‑1 (MQ=255)
        tGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGa                                                                                                                                              <  2:1358453/149‑1 (MQ=255)
        tGGCTTGAACAAGAGATGGGCATTAAGAGGTAACATCCGTATCGGTGGTCTCCACGAAACGGCGTTGAGCAATCGTCGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAGTAAAAAATCCCGCCGCTGGCGGGa                                                                                                                                              <  3:1253487/149‑1 (MQ=255)
         ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                             <  1:1009954/149‑1 (MQ=255)
         ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                             <  1:1130583/149‑1 (MQ=255)
         ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                             >  1:607917/1‑149 (MQ=255)
              gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                        <  2:1330431/149‑1 (MQ=255)
              gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                        <  3:499173/149‑1 (MQ=255)
                agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                                                                                                      >  3:903596/1‑149 (MQ=255)
                 gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAttttttatt                                                                                                                                     <  2:607917/149‑6 (MQ=255)
                  aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCaa                                                                                                                                    <  1:616703/149‑1 (MQ=255)
                  aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCaa                                                                                                                                    <  3:1237083/149‑1 (MQ=255)
                    gagaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                  <  1:230772/149‑1 (MQ=255)
                    gagaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                  >  3:573461/1‑149 (MQ=255)
                    gagaGGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                  <  2:195742/149‑1 (MQ=255)
                     agaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTg                                                                                                                                 <  2:1295836/149‑1 (MQ=255)
                        tGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                  <  2:652955/145‑1 (MQ=255)
                        tGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                  >  1:652955/1‑145 (MQ=255)
                                 agaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                     <  3:1064177/149‑1 (MQ=255)
                                   aagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                                   <  2:527279/149‑1 (MQ=255)
                                    attAACATCCGTATCGGTGTTATCCACAAAACGGCGTTGAGCAATCGACGCCGTTTTTTTTTAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGat                                                                                                                  <  4:1230257/147‑1 (MQ=255)
                                       aaCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat                                                                                                               <  4:699513/149‑1 (MQ=255)
                                        aCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCTGTTTTTTTATAGCTTATTCTTATTAAATTGTATTAAACCGGACCATAACAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAAGCTACAAAAGattatt                                                                                                              <  2:1281639/149‑1 (MQ=255)
                                          aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattatt                                                                                                              <  2:1390028/147‑1 (MQ=255)
                                          aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattatt                                                                                                              >  1:1390028/1‑147 (MQ=255)
                                          aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                                            <  3:928144/149‑1 (MQ=255)
                                           tCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                           <  4:1215451/149‑1 (MQ=255)
                                            ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                                          <  4:635913/149‑1 (MQ=255)
                                            ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                                          <  1:1497459/149‑1 (MQ=255)
                                              gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGa                                                                                                        >  3:324612/1‑149 (MQ=255)
                                                aTCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                      <  4:840936/149‑1 (MQ=255)
                                                aTCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                      <  2:757612/149‑1 (MQ=255)
                                                  cGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAg                                                                                                    <  3:38492/149‑1 (MQ=255)
                                                       gctataaaaAAACGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCaa                                                                                                                              >  3:952852‑M1/10‑118 (MQ=255)
                                                        tATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtc                                                                                              <  2:468182/149‑1 (MQ=255)
                                                         aTCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcg                                                                                             >  2:152494/1‑149 (MQ=255)
                                                         aTCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcg                                                                                             >  3:776659/1‑149 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                                          >  3:682158/1‑149 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                                          >  2:1226240/1‑149 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCTCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                         >  4:1369022/1‑149 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                         <  2:451020/149‑1 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                         <  4:1207647/149‑1 (MQ=255)
                                                              cGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                        <  2:67552/149‑1 (MQ=255)
                                                                 aacggcgtcgattgctcaacgccgtttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                     >  4:208585‑M2/28‑149 (MQ=255)
                                                                 AACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                     <  3:117246‑M1/149‑1 (MQ=255)
                                                                    GGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                                  <  4:1293209‑M1/149‑1 (MQ=255)
                                                                    GGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGATTTa                                                                                 <  4:880102‑M1/149‑1 (MQ=255)
                                                                       GTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAcc                                                                               >  1:741686‑M1/1‑149 (MQ=255)
                                                                              aaTCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCa                                                                        >  1:485089/1‑149 (MQ=255)
                                                                                tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                      <  2:1397825/149‑1 (MQ=255)
                                                                                   acgccgtttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCAt                                                                   <  3:1022458‑M2/140‑1 (MQ=255)
                                                                                     gccgtttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                                                 >  3:618178‑M2/8‑149 (MQ=255)
                                                                                     gccgtttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                                                 <  2:959066‑M2/142‑1 (MQ=255)
                                                                                     gccgtttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCACCCATTg                                                                 >  3:1372756‑M2/8‑149 (MQ=255)
                                                                                         tttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCttt                                                             >  2:334487‑M2/4‑149 (MQ=255)
                                                                                         tttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCttt                                                             >  3:192273‑M2/4‑149 (MQ=255)
                                                                                          ttTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                                            <  3:120408‑M2/147‑1 (MQ=255)
                                                                                            ttttATCGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                                          <  1:400191/149‑1 (MQ=255)
                                                                                            ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                                          <  3:1360260/149‑1 (MQ=255)
                                                                                                   gCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACa                                                   >  2:819154/1‑149 (MQ=255)
                                                                                                   gCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACa                                                   >  4:858185/1‑149 (MQ=255)
                                                                                                       attcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTg                                               <  4:1365024/149‑1 (MQ=255)
                                                                                                        ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat                                                                                                               <  3:378235/84‑1 (MQ=255)
                                                                                                        ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat                                                                                                               >  4:378235/1‑84 (MQ=255)
                                                                                                         tcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGa                                             <  4:573461/149‑1 (MQ=255)
                                                                                                           ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                                                                                       <  1:1307055/57‑1 (MQ=255)
                                                                                                            tattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTa                                          >  2:287743/1‑149 (MQ=255)
                                                                                                             attaAATTGTCTTAAACCGGACAATTAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCCATCTACAAAAGATTATTGTTTAATCTGAAAATCGACTAGGGATTTACCTTGCTCATCCATTGCTCTTTTGATTACAGTTGGAGTag                                         >  3:544784/1‑148 (MQ=255)
                                                                                                             attaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAg                                                                                                    <  2:1004140/90‑1 (MQ=255)
                                                                                                             attaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAg                                                                                                    >  1:1004140/1‑90 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                        >  3:1000187/1‑149 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                        >  3:892977/1‑149 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                        >  1:1428747/1‑149 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                        >  4:1343009/1‑149 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                        >  1:699087/1‑149 (MQ=255)
                                                                                                                  aTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCt                                    >  3:473776/1‑149 (MQ=255)
                                                                                                                  aTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATATACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCt                                    >  3:856329/1‑149 (MQ=255)
                                                                                                                   ttGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCACCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt                                   >  1:86145/1‑149 (MQ=255)
                                                                                                                    tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTg                                  <  2:800178/149‑1 (MQ=255)
                                                                                                                     gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                                 <  3:1293209/149‑1 (MQ=255)
                                                                                                                        ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCa                              <  3:70289/149‑1 (MQ=255)
                                                                                                                        ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCa                              >  4:1131567/1‑149 (MQ=255)
                                                                                                                          aaaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGt                            >  3:1193911/1‑149 (MQ=255)
                                                                                                                              cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                        >  2:775260/1‑149 (MQ=255)
                                                                                                                               ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAt                                                                                    >  4:1079182/1‑88 (MQ=255)
                                                                                                                               ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAt                                                                                    <  3:1079182/88‑1 (MQ=255)
                                                                                                                                gACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGt                      <  4:1044264/149‑1 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt                   <  4:470045/149‑1 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt                   >  2:1409276/1‑149 (MQ=255)
                                                                                                                                    gcttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                  >  4:1161791/5‑149 (MQ=255)
                                                                                                                                    aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                  <  4:384252/149‑1 (MQ=255)
                                                                                                                                    aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                  <  4:776659/149‑1 (MQ=255)
                                                                                                                                     tAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                        >  1:790597/1‑142 (MQ=255)
                                                                                                                                     tAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAg                 >  3:1358885/1‑149 (MQ=255)
                                                                                                                                      ttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGt                >  1:934507/3‑149 (MQ=255)
                                                                                                                                      aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGa                                                        >  1:26201/1‑109 (MQ=255)
                                                                                                                                      aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGt                >  1:1400367/1‑149 (MQ=255)
                                                                                                                                        aaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAg                 <  2:553824/146‑1 (MQ=255)
                                                                                                                                        aaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGttt              >  1:1378779/1‑149 (MQ=255)
                                                                                                                                          aaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg            >  2:257063/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGT  >  NC_000913/1292325‑1292621

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