Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,105 A→G 5.3% intergenic (+29/‑431) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1050AG 5.3% 436.6 / 9.5 131intergenic (+29/‑431)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (60/64);  new base G (2/5);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 136 total aligned reads displayed.

ATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTA  >  NC_000913/2986960‑2987253
                                                                                                                                                 |                                                                                                                                                     
atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGATGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑Caa                                                                                                                                                   >  1:433790/1‑149 (MQ=255)
  catGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTACATCGGAGATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAGTGTTCTTGAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAAt                                                                                                                                                  >  4:659926/1‑148 (MQ=255)
      ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAG‑CAATTTTGa                                                                                                                                             <  1:669993/149‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAG‑CAATTTTGaa                                                                                                                                            <  1:353290/149‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGaa                                                                                                                                            <  1:465558/149‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGaa                                                                                                                                            <  2:749595/149‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGa                                                                                                                                             <  4:614420/148‑1 (MQ=255)
                cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAAccattcca                                                                                                                                    >  4:1749/1‑148 (MQ=255)
                cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAAccattcca                                                                                                                                    >  3:420182/1‑148 (MQ=255)
                  aCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAAccattccatt                                                                                                                                  >  3:1045261/1‑148 (MQ=255)
                   cGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCAttt                                                                                                                                 >  3:615596/1‑148 (MQ=255)
                   cGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCAttt                                                                                                                                 >  3:32653/1‑148 (MQ=255)
                     tATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTc                                                                                                                              >  2:411926/1‑149 (MQ=255)
                       tctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCat                                                                                                                            >  1:586808/1‑149 (MQ=255)
                        ctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAAccattccatt                                                                                                                                  >  4:942586/1‑142 (MQ=255)
                            tCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAatat                                                                                                                       >  2:426146/1‑149 (MQ=255)
                             cccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCGGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAatat                                                                                                                       >  4:1001132/1‑148 (MQ=255)
                              ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTc                                                                                                                     <  1:426146/149‑1 (MQ=255)
                                ggCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                    <  4:1037267/148‑1 (MQ=255)
                                    cTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTg                                                                                                                >  4:386813/1‑148 (MQ=255)
                                     ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATTGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                              >  1:118529/1‑149 (MQ=255)
                                       aaaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                            <  1:103419/149‑1 (MQ=255)
                                        aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                            <  4:450988/148‑1 (MQ=255)
                                        aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                            >  4:839783/1‑148 (MQ=255)
                                         aCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatt                                                                                                           <  3:1057713/148‑1 (MQ=255)
                                                  ttCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAG‑CAATTTTGa                                                                                                                                             <  1:862058/105‑1 (MQ=255)
                                                  ttCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCt                                                                                                 >  2:1076174/1‑149 (MQ=255)
                                                  ttCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTG‑TAGCAATTTTGa                                                                                                                                             >  2:862058/1‑105 (MQ=255)
                                                     cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTAt                                                                                               <  4:1045261/148‑1 (MQ=255)
                                                     cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCAGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATc                                                                                              <  2:527714/149‑1 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAAtttt                                                                                                                                               <  3:1014485/99‑1 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATc                                                                                              <  4:615596/148‑1 (MQ=255)
                                                        cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                            >  3:653032/1‑148 (MQ=255)
                                                        cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACt                                                                                           <  1:985976/149‑1 (MQ=255)
                                                           aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAAGCATTCCATTTTTCataat                                                                                                                         <  3:835954/116‑1 (MQ=255)
                                                           aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCataat                                                                                                                         >  4:835954/1‑116 (MQ=255)
                                                              tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                      <  4:858487/148‑1 (MQ=255)
                                                              tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                      <  4:44425/148‑1 (MQ=255)
                                                              tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                      <  4:612571/148‑1 (MQ=255)
                                                              tCTGATTACGATGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                                     <  2:433790/149‑1 (MQ=255)
                                                                tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTa                                                                                                <  1:263006/136‑1 (MQ=255)
                                                                tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGt                                                                                   <  2:229766/149‑1 (MQ=255)
                                                                      cGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATc                                                                             <  1:585314/149‑1 (MQ=255)
                                                                       gCTGATCGTGGGTATCCTTTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGCACCATTCCATTTTTCATAATATTCAt                                                                                                                   >  1:312156/1‑110 (MQ=255)
                                                                       gCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGAATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCAt                                                                                                                   <  2:312156/110‑1 (MQ=255)
                                                                          gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                          <  3:739750/148‑1 (MQ=255)
                                                                          gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                          >  3:1020775/1‑148 (MQ=255)
                                                                          gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAg                                                                         >  1:1011720/1‑149 (MQ=255)
                                                                          gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACACTATCACTTCTCTCGTCATATCCTa                                                                          >  4:26303/1‑148 (MQ=255)
                                                                            tCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                          >  4:517496/1‑114 (MQ=255)
                                                                               tGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatta                                                                    >  2:562271/1‑149 (MQ=255)
                                                                               tGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatt                                                                     >  3:883645/1‑148 (MQ=255)
                                                                               tGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatt                                                                     >  4:705722/1‑148 (MQ=255)
                                                                               tGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatt                                                                     <  3:598934/148‑1 (MQ=255)
                                                                                 ggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                  <  2:958410/149‑1 (MQ=255)
                                                                                      ccTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtta                                                             <  1:445316/149‑1 (MQ=255)
                                                                                      ccTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCGA‑CAATTTTGGACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCCAGTATTATCagtt                                                              <  3:1001132/148‑1 (MQ=255)
                                                                                         gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagtt                                                          <  2:688680/149‑1 (MQ=255)
                                                                                         gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                              >  2:587129/1‑145 (MQ=255)
                                                                                             gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                      >  1:64783/1‑149 (MQ=255)
                                                                                              tATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                                     <  1:1032210/149‑1 (MQ=255)
                                                                                               aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                  >  1:680629/1‑135 (MQ=255)
                                                                                               aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                  <  2:680629/135‑1 (MQ=255)
                                                                                               aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacac                                                    <  2:586808/149‑1 (MQ=255)
                                                                                                 cGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacaca                                                   >  3:889749/1‑148 (MQ=255)
                                                                                                     aaTGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAag                                              <  2:843852/149‑1 (MQ=255)
                                                                                                          tCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaa                                                        <  3:864647/134‑1 (MQ=255)
                                                                                                          tCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                         <  1:411926/149‑1 (MQ=255)
                                                                                                            tttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAag                                              <  1:560571/142‑1 (MQ=255)
                                                                                                             ttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaaa                                               <  3:21171/140‑1 (MQ=255)
                                                                                                               aGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGaa                                     <  4:1020775/148‑1 (MQ=255)
                                                                                                                  tAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGGACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaaa                                               <  2:1058718/135‑1 (MQ=255)
                                                                                                                  tAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                           >  4:569825/1‑139 (MQ=255)
                                                                                                                  tAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                           <  3:569825/139‑1 (MQ=255)
                                                                                                                    atatCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                                >  4:120123/1‑148 (MQ=255)
                                                                                                                    atatCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                                <  3:391110/148‑1 (MQ=255)
                                                                                                                       tCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAATCACACAAAGAGATCTGAACATCTTaaa                            <  1:128606/149‑1 (MQ=255)
                                                                                                                       tCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaaa                            <  2:354139/149‑1 (MQ=255)
                                                                                                                       tCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                             >  3:211582/1‑148 (MQ=255)
                                                                                                                         aGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAAt                           <  3:223346/148‑1 (MQ=255)
                                                                                                                          gTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATaa                         >  1:968005/1‑149 (MQ=255)
                                                                                                                          gTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATa                          >  3:983447/1‑148 (MQ=255)
                                                                                                                          gTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATa                          <  3:1048159/148‑1 (MQ=255)
                                                                                                                           ttGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATaa                         >  3:614232/1‑148 (MQ=255)
                                                                                                                            tgtgGATGGTTCAAAATTGCTA‑CAATTATGCACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATGCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                        >  4:648950/1‑148 (MQ=255)
                                                                                                                               ggATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTa                    >  2:264593/1‑149 (MQ=255)
                                                                                                                                  tGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATACCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAttt                 <  1:827051/149‑1 (MQ=255)
                                                                                                                                    gTTCAAAATTGTAG‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGa                       >  2:470901/1‑141 (MQ=255)
                                                                                                                                     ttCAAAATTG‑TAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  2:52298/1‑149 (MQ=255)
                                                                                                                                      tCAAAATTGTAG‑CAATTTTGAACCATCCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtattat                                                                   >  4:1058383/1‑95 (MQ=38)
                                                                                                                                      tCAAAATTGTAG‑CAATTTTGAACCATCCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtattat                                                                   <  3:1058383/95‑1 (MQ=38)
                                                                                                                                      tCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  3:101662/1‑148 (MQ=255)
                                                                                                                                      tCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa             <  1:783705/149‑1 (MQ=255)
                                                                                                                                      tCAAAATTG‑TAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  4:728471/1‑148 (MQ=255)
                                                                                                                                         aaaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATg           <  3:705722/148‑1 (MQ=255)
                                                                                                                                         aaaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATg           <  4:883645/148‑1 (MQ=255)
                                                                                                                                         aaaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGc          >  2:447823/1‑149 (MQ=255)
                                                                                                                                           aTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCAt        <  2:486829/149‑1 (MQ=255)
                                                                                                                                                 a‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACtt   >  4:216611/1‑148 (MQ=255)
                                                                                                                                                 a‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTa  <  1:1074458/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTA  >  NC_000913/2986960‑2987253

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