Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,049 A→T 7.1% intergenic (+17/+79) trpT → / ← hdfR tRNA‑Trp/flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0490AT 7.1% 338.5 / 12.7 98intergenic (+17/+79)trpT/hdfRtRNA‑Trp/flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base T (2/5);  total (45/53)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

Only 100 of 101 total aligned reads displayed.

CCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGC  >  NC_000913/3946903‑3947197
                                                                                                                                                  |                                                                                                                                                    
cctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaaa                                                                                                                                                    <  3:397435/149‑1 (MQ=255)
cctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaaa                                                                                                                                                    <  1:1781478/149‑1 (MQ=255)
   gTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCt                                                                                                                                                 <  4:763936/149‑1 (MQ=255)
     cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  2:976192/149‑1 (MQ=255)
     cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  3:566530/149‑1 (MQ=255)
     cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  1:143007/149‑1 (MQ=255)
     cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  4:1574284/149‑1 (MQ=255)
     cACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCAACGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  4:1488105/149‑1 (MQ=255)
      aCGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                              >  3:1214165/1‑149 (MQ=255)
      aCGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                              >  2:1176289/1‑149 (MQ=255)
      aCGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                              <  2:107963/149‑1 (MQ=255)
      aCGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAg                                                                                                                                              <  3:1641269/149‑1 (MQ=255)
        gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            <  3:1514649/149‑1 (MQ=255)
        gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            <  3:1149696/149‑1 (MQ=255)
        gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            <  1:182360/149‑1 (MQ=255)
        gCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGa                                                                                                                                            >  1:1442578/1‑149 (MQ=255)
         cAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAt                                                                                                                                           <  1:1176289/149‑1 (MQ=255)
             gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttttt                                                                                                                                       <  3:95402/149‑1 (MQ=255)
               gTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttttttt                                                                                                                                     <  2:1442578/149‑1 (MQ=255)
                 cgcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                   <  4:1214165/149‑1 (MQ=255)
                          gAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            >  4:1505248/1‑131 (MQ=255)
                          gAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            <  3:1505248/131‑1 (MQ=255)
                           aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttt                                                                                                                                        <  1:1040905/134‑1 (MQ=255)
                           aGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttt                                                                                                                                        >  2:1040905/1‑134 (MQ=255)
                             tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            <  3:1502364/128‑1 (MQ=255)
                             tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                            >  4:1502364/1‑128 (MQ=255)
                             tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaa                                                                                                                             <  1:1655159/143‑1 (MQ=255)
                             tCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaa                                                                                                                             >  2:1655159/1‑143 (MQ=255)
                                     gtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               >  1:470501/1‑117 (MQ=255)
                                     gtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  2:470501/117‑1 (MQ=255)
                                         cGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                           >  1:457089/1‑149 (MQ=255)
                                               ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                   <  3:1129339/119‑1 (MQ=255)
                                               ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                                                                                                   >  4:1129339/1‑119 (MQ=255)
                                               ccTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCt                                                                                                     >  4:1115336/1‑149 (MQ=255)
                                                 tAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                                                                                                   <  3:1115336/149‑1 (MQ=255)
                                                   aTTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTa                                                                                                 >  2:1560379/1‑149 (MQ=255)
                                                    ttAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAc                                                                                                <  4:1322326/149‑1 (MQ=255)
                                                         ggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTGGCGAAAGCTAAGGATTTTTTTTATCTGAAATAAc                                                                                                                        <  2:330983/120‑1 (MQ=255)
                                                         ggCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTCTCTGAAATAAc                                                                                                                        >  1:330983/1‑120 (MQ=255)
                                                               gTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTAt                                                                                                                                  <  1:68389/104‑1 (MQ=255)
                                                               gTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTAt                                                                                                                                  >  2:68389/1‑104 (MQ=255)
                                                                ttCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCt                                                                                                                                >  4:127808/1‑105 (MQ=255)
                                                                ttCAATTGGTAGAGCACCGGTCTCCAAAACCGCGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCt                                                                                                                                <  3:127808/105‑1 (MQ=255)
                                                                  cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcca                                                                                  <  1:1560379/149‑1 (MQ=255)
                                                                       ggTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAt                                                                             <  2:999848/149‑1 (MQ=255)
                                                                              cACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               >  2:917914/1‑76 (MQ=255)
                                                                              cACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  1:917914/76‑1 (MQ=255)
                                                                                   gTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATac                                                                 >  2:496023/1‑149 (MQ=255)
                                                                                      tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttttt                                                                                                                                       <  2:713552/76‑1 (MQ=255)
                                                                                      tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttttt                                                                                                                                       >  1:713552/1‑76 (MQ=255)
                                                                                      tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       <  1:1651634/140‑1 (MQ=255)
                                                                                      tcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       >  2:1651634/1‑140 (MQ=255)
                                                                                        cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                            >  3:205164/1‑149 (MQ=255)
                                                                                         aaaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               >  3:201261/1‑65 (MQ=255)
                                                                                         aaaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTaa                                                                                                                                               <  4:201261/65‑1 (MQ=255)
                                                                                          aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCCTTCATATTAATACACTTCAt                                                          >  4:745448/1‑149 (MQ=255)
                                                                                           aaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATc                                                         <  1:496023/149‑1 (MQ=255)
                                                                                             ccGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCata                                                                        <  4:275121/132‑1 (MQ=255)
                                                                                             ccGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCata                                                                        >  3:275121/1‑132 (MQ=255)
                                                                                                    ttGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCg                                                                                              >  3:113991/1‑103 (MQ=255)
                                                                                                    ttGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCg                                                                                              <  4:113991/103‑1 (MQ=255)
                                                                                                     tGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCAt                                                          >  2:134917/1‑138 (MQ=255)
                                                                                                     tGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCAt                                                          <  1:134917/138‑1 (MQ=255)
                                                                                                      gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       >  4:282935/1‑124 (MQ=255)
                                                                                                      gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       <  3:282935/124‑1 (MQ=255)
                                                                                                      gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACt                                                              >  1:1422839/1‑133 (MQ=255)
                                                                                                      gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACt                                                              <  2:1422839/133‑1 (MQ=255)
                                                                                                           ttCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                   <  2:772606/123‑1 (MQ=255)
                                                                                                           ttCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                   >  1:772606/1‑123 (MQ=255)
                                                                                                              gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                                                                                                   >  3:912564/1‑88 (MQ=255)
                                                                                                              gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                                                                                                   <  4:912564/88‑1 (MQ=255)
                                                                                                               aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     <  4:205164/149‑1 (MQ=255)
                                                                                                                gTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATAtt                                                                      <  2:3830/115‑1 (MQ=255)
                                                                                                                gTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATAtt                                                                      >  1:3830/1‑115 (MQ=255)
                                                                                                                gTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAg                                    >  4:757363/1‑149 (MQ=255)
                                                                                                                  ctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag                                                                                 <  3:1449705/102‑1 (MQ=255)
                                                                                                                  ctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag                                                                                 >  4:1449705/1‑102 (MQ=255)
                                                                                                                          cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       >  3:336104/1‑104 (MQ=255)
                                                                                                                          cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       <  4:336104/104‑1 (MQ=255)
                                                                                                                            cTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCaa                        >  1:1651228/1‑149 (MQ=255)
                                                                                                                                 aGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAt                   <  3:745448/149‑1 (MQ=255)
                                                                                                                                  gAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     >  1:892303/1‑130 (MQ=255)
                                                                                                                                  gAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     <  2:892303/130‑1 (MQ=255)
                                                                                                                                       catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTAt             <  2:457089/149‑1 (MQ=255)
                                                                                                                                        atcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAcc                                                                                               >  4:1713657/1‑66 (MQ=255)
                                                                                                                                        atcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAcc                                                                                               <  3:1713657/66‑1 (MQ=255)
                                                                                                                                         tcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                  <  2:1409676/126‑1 (MQ=255)
                                                                                                                                         tcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATAAACTTCATCCAGCACGTTAATTTTCAAAAGAt                                  >  1:1409676/1‑126 (MQ=255)
                                                                                                                                            ttAGCGAGAGCTCAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTGATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTa        >  4:1602249/1‑149 (MQ=255)
                                                                                                                                            ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCat                                                                         <  2:972827/84‑1 (MQ=255)
                                                                                                                                            ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCat                                                                         >  1:972827/1‑84 (MQ=255)
                                                                                                                                            ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTGATACACTTCATCCAGCACGTTAATTTTCaaaa                                     >  4:388910/1‑120 (MQ=255)
                                                                                                                                            ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     <  3:388910/120‑1 (MQ=255)
                                                                                                                                            ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAATGCATTTTTATCGCTa        >  2:546664/1‑149 (MQ=255)
                                                                                                                                              aGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                          <  3:653808/81‑1 (MQ=255)
                                                                                                                                              aGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                          >  4:653808/1‑81 (MQ=255)
                                                                                                                                              aGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTAtt      <  2:1651228/149‑1 (MQ=255)
                                                                                                                                                 gAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg   >  1:898531/1‑149 (MQ=255)
                                                                                                                                                 gAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg   >  4:1580317/1‑149 (MQ=255)
                                                                                                                                                  aaaGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc  >  1:1835741/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGC  >  NC_000913/3946903‑3947197

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