Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,104 T→A 6.8% intergenic (+28/‑432) 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,1040TA 6.8% 267.7 / 9.2 104intergenic (+28/‑432)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (47/50);  new base A (3/4);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

Only 100 of 106 total aligned reads displayed.

TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTA  >  NC_000913/2986955‑2987253
                                                                                                                                                     |                                                                                                                                                     
ttttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCt                                                                                                                                                       <  2:733594/150‑1 (MQ=255)
  ttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTAc                                                                                                                                                     >  1:338592/1‑150 (MQ=255)
         tGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTg                                                                                                                                              <  1:407044/150‑1 (MQ=255)
          gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGa                                                                                                                                             <  1:1096989/150‑1 (MQ=255)
          gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGa                                                                                                                                             >  1:584370/1‑150 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                            <  2:1080896/150‑1 (MQ=255)
            cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                           <  2:1051334/150‑1 (MQ=255)
            cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                           <  1:1080655/150‑1 (MQ=255)
             tACGCTGTCTACGTATCTCTTCCCGGCGCTTAATCACTGGCTTTGCATCTGAAATCTGATTAGGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTTCTACAATTTTGAAcc                                                                                                                                          <  2:538731/150‑1 (MQ=255)
             tACGCTGTCAACGTATCTCATCCCGGCGCGTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAcc                                                                                                                                          <  1:1071316/150‑1 (MQ=255)
             tACGCTGTCAACGTATCTCATCCCGGCGCCTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAcc                                                                                                                                          <  1:434393/150‑1 (MQ=255)
                 cTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccatta                                                                                                                                      >  1:305549/1‑149 (MQ=255)
                  tGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcc                                                                                                                                     <  1:1189526/150‑1 (MQ=255)
                     cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATCGTTCAAAATTGCTACAATTTTGAAccattcaatt                                                                                                                                  >  2:185285/1‑150 (MQ=255)
                      aaCGTATCTCCTCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGAACGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTAGTGGATGGTTCAAAATTGCTAAAATATTGAACCCTTCCAttt                                                                                                                                 >  2:805078/1‑150 (MQ=255)
                          tATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCa                                                                                                                             <  2:1106962/150‑1 (MQ=255)
                             ctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataa                                                                                                                          >  1:567466/1‑150 (MQ=255)
                                       cgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAtttt                                                                                                                >  2:1272420/1‑149 (MQ=255)
                                         cTTAAACCCTGGCGTTCCATCGGAAATCTGGTTACGCTGATCGTGGGTATCCTGTGCGTATCGGTTATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGATACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                              >  2:83020/1‑150 (MQ=255)
                                          ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                             <  1:623521/150‑1 (MQ=255)
                                          ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                             <  2:609985/150‑1 (MQ=255)
                                           tAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                            >  1:402027/1‑150 (MQ=255)
                                            aaaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTAATTTGttatt                                                                                                           >  2:1297892/1‑150 (MQ=255)
                                             aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                          <  1:889256/150‑1 (MQ=255)
                                              aCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGTGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTtattat                                                                                                         >  1:950489/1‑150 (MQ=255)
                                               cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                        >  2:1227190/1‑150 (MQ=255)
                                               cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                        >  1:59484/1‑150 (MQ=255)
                                                 cTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTcc                                                                                                      <  2:409339/150‑1 (MQ=255)
                                                   ggCTTTCCATCGTTAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATAAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                    <  1:116228/150‑1 (MQ=255)
                                                   ggCGTTCCATCGGAAATCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                    <  1:566039/150‑1 (MQ=255)
                                                   ggCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                    <  2:399990/150‑1 (MQ=255)
                                                    gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                                   >  1:915397/1‑150 (MQ=255)
                                                        tCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTAt                                                                                               >  1:983645/1‑150 (MQ=255)
                                                         ccATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATc                                                                                              >  1:280496/1‑150 (MQ=255)
                                                          cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCa                                                                                             >  2:718602/1‑150 (MQ=255)
                                                           aTCGGAAATCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                            <  1:1247656/150‑1 (MQ=255)
                                                             cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                          >  1:277577/1‑150 (MQ=255)
                                                                aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctc                                                                                       >  2:566133/1‑150 (MQ=255)
                                                                  aTCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                                     <  1:1310183/150‑1 (MQ=255)
                                                                   tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGATACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACACTATCACTTCTCTCg                                                                                    >  2:1193907/1‑150 (MQ=255)
                                                                     tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                  <  1:880489/150‑1 (MQ=255)
                                                                      gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGGTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAAAATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                                 <  2:555535/150‑1 (MQ=255)
                                                                         tACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCatat                                                                              <  1:780910/150‑1 (MQ=255)
                                                                           cacTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATcc                                                                            <  1:532542/148‑1 (MQ=255)
                                                                           cGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATcc                                                                            >  2:802608/1‑150 (MQ=255)
                                                                              tGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAg                                                                         <  2:321965/150‑1 (MQ=255)
                                                                               gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGt                                                                        >  2:1065230/1‑150 (MQ=255)
                                                                                ttCGTGGGTATCCTGTGCGTATAGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGta                                                                       <  1:129431/149‑1 (MQ=255)
                                                                                 tCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtat                                                                      <  2:685244/150‑1 (MQ=255)
                                                                                 tCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtat                                                                      <  1:1053327/150‑1 (MQ=255)
                                                                                   gTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatta                                                                    <  1:106581/150‑1 (MQ=255)
                                                                                     gggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                  >  2:769732/1‑150 (MQ=255)
                                                                                       gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                                >  2:932335/1‑150 (MQ=255)
                                                                                       gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                                <  2:265185/150‑1 (MQ=255)
                                                                                         aTCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                              <  2:915397/150‑1 (MQ=255)
                                                                                           ccTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttat                                                            >  2:495072/1‑149 (MQ=255)
                                                                                             tgtgAGTATCGGTAATTTTCTTTAGCTAATATCAGTTGTGGTTGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagatagtt                                                          <  2:1155054/150‑1 (MQ=255)
                                                                                              gtgCGTATAGGTAATGTTCTTTAGCTTATATCAGTTGTTGATTGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAATTCTCTCGTCATATCCTAGTATTATCAgttagtta                                                         <  1:1068828/150‑1 (MQ=255)
                                                                                                gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                       <  1:1242024/150‑1 (MQ=255)
                                                                                                 cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                      <  2:584370/150‑1 (MQ=255)
                                                                                                  gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                                     <  2:1231332/150‑1 (MQ=255)
                                                                                                    aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCACGTCATATCCTAGTATTAGCAGTTAGTTAAacacc                                                   >  2:349787/1‑149 (MQ=255)
                                                                                                    aTCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCAGTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacaca                                                   >  2:344848/1‑150 (MQ=255)
                                                                                                      cGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAcacaca                                                 <  2:3838/150‑1 (MQ=255)
                                                                                                      cGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAcacaca                                                 >  1:925322/1‑150 (MQ=255)
                                                                                                            tGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                           >  1:1099744/1‑150 (MQ=255)
                                                                                                             gTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCATAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                          >  2:927051/1‑150 (MQ=255)
                                                                                                              ttCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                         <  2:983645/150‑1 (MQ=255)
                                                                                                                cTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTg                                       <  2:338592/150‑1 (MQ=255)
                                                                                                                atttATCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTg                                       <  1:426360/149‑1 (MQ=255)
                                                                                                                  ttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGaa                                     <  2:140405/150‑1 (MQ=255)
                                                                                                                   tAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAAc                                    <  2:363853/150‑1 (MQ=255)
                                                                                                                    aGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACa                                   >  2:125361/1‑150 (MQ=255)
                                                                                                                    aGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACa                                   >  2:1291461/1‑150 (MQ=255)
                                                                                                                     gCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACAt                                  >  1:176734/1‑150 (MQ=255)
                                                                                                                      cTAATATCAGTTGTGGATGGTTCAAAATTTCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                                <  1:984323/150‑1 (MQ=255)
                                                                                                                         atatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTTTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                             <  2:1306536/150‑1 (MQ=255)
                                                                                                                          tatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                             <  2:402027/150‑1 (MQ=255)
                                                                                                                            tCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGCACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTAGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAACTGAACAACTTAAAt                           >  2:801731/1‑150 (MQ=255)
                                                                                                                               gTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                        <  2:1268163/150‑1 (MQ=255)
                                                                                                                               gTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                        >  1:1174620/1‑150 (MQ=255)
                                                                                                                                 tgtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGAc                      <  2:1069942/150‑1 (MQ=255)
                                                                                                                                  gtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACACTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACt                     >  2:927775/1‑150 (MQ=255)
                                                                                                                                   tgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTa                    >  2:1307565/1‑150 (MQ=255)
                                                                                                                                        ggTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaa               >  1:831271/1‑150 (MQ=255)
                                                                                                                                         gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  1:1057531/1‑150 (MQ=255)
                                                                                                                                         gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  1:118139/1‑150 (MQ=255)
                                                                                                                                         gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaat              >  2:201821/1‑149 (MQ=255)
                                                                                                                                         gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  2:20494/1‑150 (MQ=255)
                                                                                                                                         gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  2:967533/1‑150 (MQ=255)
                                                                                                                                         gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGACATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg              >  2:1134888/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa             <  1:964649/150‑1 (MQ=255)
                                                                                                                                           tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt            <  1:315528/150‑1 (MQ=255)
                                                                                                                                           tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCAACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACTCAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt            <  2:33501/150‑1 (MQ=255)
                                                                                                                                            cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATg           <  1:688232/150‑1 (MQ=255)
                                                                                                                                            cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTAATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGACCATCTTAAATAAGCATATTTAATATg           >  2:531387/1‑150 (MQ=255)
                                                                                                                                             aaaaTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgt          <  2:692046/150‑2 (MQ=255)
                                                                                                                                               aaTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgttc        <  2:850293/150‑4 (MQ=255)
                                                                                                                                                   gCTACAATTTTGAACCATTCCATTTTTCATAATATTAATTTGTTATTATTCCACCCTATACCTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTTAACACACATAGCGATCTGAACATCTTAAATATGCCTATTTAAGATGGCTAACt    >  2:185279/1‑150 (MQ=255)
                                                                                                                                                     tACAATTTTGAACCATTACATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTa  <  1:332035/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTA  >  NC_000913/2986955‑2987253

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