Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,554,283 T→A 9.0% L400Q (CTG→CAG)  malT → mal regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,2830TA 9.0% 1055.4 / 9.1 380L400Q (CTG→CAG) malTmal regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (127/201);  new base A (46/2);  total (177/203)
Fisher's exact test for biased strand distribution p-value = 3.29e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 384 total aligned reads displayed.

CTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAG  >  NC_000913/3554134‑3554432
                                                                                                                                                     |                                                                                                                                                     
cTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAATGCCCt                                                                                                                                                       >  1:3596002/1‑150 (MQ=255)
cTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCt                                                                                                                                                       <  2:751439/150‑1 (MQ=255)
  ggATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGc                                                                                                                                                     >  1:985316/1‑150 (MQ=255)
   gATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGcc                                                                                                                                                    >  1:2207555/1‑150 (MQ=255)
      ggCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTg                                                                                                                                                 <  2:543762/150‑1 (MQ=255)
       gCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTgg                                                                                                                                                <  2:2410482/150‑1 (MQ=255)
         ccAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGa                                                                                                                                              >  1:1988008/1‑150 (MQ=255)
          cAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGAc                                                                                                                                             >  2:949473/1‑150 (MQ=255)
          cAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGAc                                                                                                                                             <  1:3433937/150‑1 (MQ=255)
          cAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGAc                                                                                                                                             >  1:2882827/1‑150 (MQ=255)
           ccagttttCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGCAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGAcc                                                                                                                                            >  2:3279852/6‑149 (MQ=255)
           aGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACa                                                                                                                                            <  1:3267875/150‑1 (MQ=255)
                 ttCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGc                                                                                                                                      <  1:1236716/150‑1 (MQ=255)
                   cccAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTg                                                                                                                                    <  2:313310/150‑1 (MQ=255)
                     cAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCGTAAGGCCCAGCCGTGGGACAGCTTGCTGGa                                                                                                                                  >  2:605662/1‑150 (MQ=255)
                      aGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGCGGGACAGCTTGCTGGaa                                                                                                                                 >  2:1161938/1‑150 (MQ=255)
                       gCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGaaa                                                                                                                                >  2:412751/1‑150 (MQ=255)
                          aaGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATc                                                                                                                             >  2:2204698/1‑150 (MQ=255)
                          aaGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCGTAAGGCCCTGCCGTGGGACAGCTTTCTGGAAAATc                                                                                                                             >  2:1750340/1‑150 (MQ=255)
                           aGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATcc                                                                                                                            <  1:3678850/150‑1 (MQ=255)
                            gCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGTTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCg                                                                                                                           <  2:1395006/150‑1 (MQ=255)
                            gCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCg                                                                                                                           <  1:1785994/150‑1 (MQ=255)
                            gCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCg                                                                                                                           <  2:3779725/150‑1 (MQ=255)
                             cAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGc                                                                                                                          >  1:1605383/1‑150 (MQ=255)
                                ttCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGt                                                                                                                       >  1:2795327/1‑150 (MQ=255)
                                 tcatcaTGCGCTGGCGGCGGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGtt                                                                                                                      <  1:2937551/150‑1 (MQ=255)
                                   atcatGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTgg                                                                                                                    >  2:50863/1‑150 (MQ=255)
                                   atcatGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTgg                                                                                                                    >  2:655706/1‑150 (MQ=255)
                                    tcatGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGgt                                                                                                                   >  1:3158935/1‑150 (MQ=255)
                                    tcatGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGgt                                                                                                                   <  2:1833066/150‑1 (MQ=255)
                                          gcTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTAtt                                                                                                             <  1:3832281/150‑1 (MQ=255)
                                            tGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGc                                                                                                           >  2:2921113/1‑150 (MQ=255)
                                            tGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGc                                                                                                           >  2:2780703/1‑150 (MQ=255)
                                            tGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGc                                                                                                           >  1:2686813/1‑150 (MQ=255)
                                                 gcAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCg                                                                                                      >  2:1664033/1‑150 (MQ=255)
                                                 gcAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCg                                                                                                      >  1:257898/1‑150 (MQ=255)
                                                   aGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCGTGCTGGGAAATCCGCAGTTGGTGTTATTGCAGGCGTg                                                                                                    >  1:1204335/1‑150 (MQ=255)
                                                     gCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGc                                                                                                  <  1:509703/150‑1 (MQ=255)
                                                        aTGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGa                                                                                               <  2:2244439/150‑1 (MQ=255)
                                                        aTGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGa                                                                                               <  2:709407/150‑1 (MQ=255)
                                                        aTGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTAATGCAGGCGTGGCAGa                                                                                               >  2:2607552/1‑150 (MQ=255)
                                                         tGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGAt                                                                                              <  2:186654/150‑1 (MQ=255)
                                                          gcgcTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCGGCCGTGGGACAGCGTGCTGGCAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATg                                                                                             >  1:3381773/1‑150 (MQ=255)
                                                             cTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTTCAGGCGTGGCTGATGCaa                                                                                          >  2:254913/1‑150 (MQ=255)
                                                                aTGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGc                                                                                       <  2:2396132/150‑1 (MQ=255)
                                                                aTGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGc                                                                                       <  1:727927/150‑1 (MQ=255)
                                                                aTGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCGTGCAGGAAAAAACGCCGTTTGTGTTATTTCAGGCGTGGCTGATGCAAAGc                                                                                       >  1:1112700/1‑150 (MQ=255)
                                                                   ctgcGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCaa                                                                                    <  1:1080661/150‑1 (MQ=255)
                                                                   ctgcGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCaa                                                                                    <  2:796100/150‑1 (MQ=255)
                                                                    tgcGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCGTGCTGGAAAATCCGCAGTTGGTGGTATTGCAGGCGTGGCTGATGCAAAGCCAAc                                                                                   >  2:1566977/1‑150 (MQ=255)
                                                                    tgcGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACCGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAAc                                                                                   >  2:3382230/1‑150 (MQ=255)
                                                                       gcgATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATc                                                                                <  1:4097702/150‑1 (MQ=255)
                                                                       gcgATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATc                                                                                <  1:244827/150‑1 (MQ=255)
                                                                        cgATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCg                                                                               <  2:2012565/150‑1 (MQ=255)
                                                                         gATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGAAAATCCCCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGc                                                                              >  1:3650075/1‑150 (MQ=255)
                                                                           tatTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTa                                                                            <  2:3083055/150‑1 (MQ=255)
                                                                            atTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTAc                                                                           <  2:3596002/150‑1 (MQ=255)
                                                                              tCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACgg                                                                         <  1:3137596/150‑1 (MQ=255)
                                                                              tCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACgg                                                                         <  1:1003988/150‑1 (MQ=255)
                                                                                  cTTAATCTGGCCTGGTGTCTGTTCAAGCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATGCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGaa                                                                     <  2:2295856/150‑1 (MQ=255)
                                                                                  cTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGaa                                                                     <  2:3430080/150‑1 (MQ=255)
                                                                                  cTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGaa                                                                     <  1:3944218/150‑1 (MQ=255)
                                                                                      aTCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTa                                                                 <  2:879605/150‑1 (MQ=255)
                                                                                       tCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTaa                                                                <  2:390306/150‑1 (MQ=255)
                                                                                        cACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGGTGCAAATCCAACATCGCTACCGCGAAGTTGac                                                               >  1:3870077/1‑150 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGGAAATACGCAGTTTGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACAAAGCTACGGCGAAGTTaaaa                                                              >  1:670404/1‑148 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAaca                                                              <  1:1381089/150‑1 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAaca                                                              <  2:2976055/150‑1 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCGGCCGTGGGACAGCTTTCTGGAAAATACGCAGGTGGTGTTATTGCAGGCGTGGCTGATTCAAAGCCCACAACGCTACGGCGAAGATAaca                                                              >  1:1864785/1‑150 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGCCAGCTTGCTGGAAAATCCCCAGGTTGTGTTATTTCAGGCGTGGCTGATGCAAAGCCAACAACGCTACGGCGAAGTTAaca                                                              >  1:4252086/1‑150 (MQ=255)
                                                                                         aCGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCAGGTAAATACCCAGTTGGTGTTATTGCAGGCGGGGCTGATGCAAATCCAAAAACGCAACGGCGAAGTTaaaa                                                              >  1:1778688/1‑148 (MQ=255)
                                                                                           gCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAAAAcc                                                            >  1:2716307/1‑150 (MQ=255)
                                                                                           gCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGCAGTTAACAcc                                                            >  2:405764/1‑150 (MQ=255)
                                                                                              tGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCGTAAGGCCCAGCCCTGGGACAGCGTGCTGGAAAAAACCCAGGTGGTTTTATTTCAGGCGTGGGTGATGCAAAGCCAACAAACCTACGGCGAACTTAACAcccag                                                         >  1:1283398/1‑148 (MQ=255)
                                                                                                   tCTGTTCATCCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGTAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGc                                                    <  2:2742688/150‑1 (MQ=255)
                                                                                                     tGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGGCAGCTTGCTGGTAAATACCCAGTTGGTGTTATTTCAGGCGTGGCTGATGCAAAGCCAACAAAGCCACGGCGAAGTTAACACCCAGCGAGccc                                                  >  1:2116287/1‑150 (MQ=255)
                                                                                                        tCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTg                                               <  1:3760581/150‑1 (MQ=255)
                                                                                                        tCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTg                                               <  1:1446063/150‑1 (MQ=255)
                                                                                                        tCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTg                                               <  1:4056570/150‑1 (MQ=255)
                                                                                                        tCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTg                                               <  1:1808131/150‑1 (MQ=255)
                                                                                                             cATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGaa                                          <  1:944757/150‑1 (MQ=255)
                                                                                                             cATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGaa                                          <  1:1703926/150‑1 (MQ=255)
                                                                                                              aTAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAAc                                         <  2:2468886/150‑1 (MQ=255)
                                                                                                               tAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACa                                        >  2:2149335/1‑150 (MQ=255)
                                                                                                               tAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACa                                        <  1:3069636/150‑1 (MQ=255)
                                                                                                                 gagAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTTACGTTGGAAACATTTCTGGAAAAAAAGCAATTGGTGTTATTTAAGGCGTGGCTGATGCAAATACAACATACCTAAGGCAAAGTTAACAACAAGCTAGAAAGTGCTGAACATg                                      >  2:3663917/3‑150 (MQ=255)
                                                                                                                    aaCTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTTCAGGCGTGGCTGATGCAAAGCCAACATCGCAACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGaaa                                   >  2:395426/1‑150 (MQ=255)
                                                                                                                         tCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAg                              <  2:457789/150‑1 (MQ=255)
                                                                                                                         tCGCTGCTGGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAg                              >  2:1196090/1‑150 (MQ=255)
                                                                                                                                tgGAAGAGTCGCTTAAGGCCCAGCCGTGGGACAGCTTGCTGGAAAATCCCCAGTTGGTGTTAATGCAGGCGGGGCAGATGCAAAGCCAACAACCCAACGGCGAAGTTAACACCCAGCAATCCCGTGCGGAACATGAAATCAAGGacaaca                       >  2:1782000/1‑150 (MQ=255)
                                                                                                                                  gAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCaga                     <  1:2380400/150‑1 (MQ=255)
                                                                                                                                  gAAGAGTCGCTTAAGGCCCTGCCCTGGGACAGCATGCTGGAAAAAACACAGTTTGTTTTATTGCAGGCGTTGCTGATGCAAAGCCAACATCGCAACGGCGAAGTTAAAAACCTGCTAGCCCCTGCGGAACATGAAAACCAGGACaaaaga                     >  2:251084/1‑145 (MQ=255)
                                                                                                                                     gagTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGagaa                  <  1:1155981/150‑1 (MQ=255)
                                                                                                                                           cTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGTCTGATGCAAAGCCAACATCTCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACAcc            <  1:2263961/150‑1 (MQ=255)
                                                                                                                                            ttAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCa           >  2:3578445/1‑150 (MQ=255)
                                                                                                                                              aaGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATg         <  2:1602044/150‑1 (MQ=255)
                                                                                                                                                 gCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCAc      <  2:2527400/150‑1 (MQ=255)
                                                                                                                                                  cccTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACg     >  2:1104350/1‑150 (MQ=255)
                                                                                                                                                  cccTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACg     >  2:3861446/1‑150 (MQ=255)
                                                                                                                                                     tGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAg  <  1:617746/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGGATGGCCCAGGGATTTCCCAGCGAAGCAATTCATCATGCGCTGGCGGCAGGCGATGCGCTGATGCTGCGCGATATTCTGCTTAATCACGCCTGGAGTCTGTTCAACCATAGCGAACTGTCGCTGCTGGAAGAGTCGCTTAAGGCCCTGCCGTGGGACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAG  >  NC_000913/3554134‑3554432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: