Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,826 T→C 9.4% intergenic (‑259/+27) bglG ← / ← phoU transcriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,8260TC 9.4% 388.4 / 24.3 141intergenic (‑259/+27)bglG/phoUtranscriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (63/64);  new base C (7/6);  total (70/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 144 total aligned reads displayed.

ATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAAT  >  NC_000913/3906679‑3906975
                                                                                                                                                   |                                                                                                                                                     
ataACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTa                                                                                                                                                     >  2:275155/1‑150 (MQ=255)
  aaCGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAg                                                                                                                                                   <  2:1376367/150‑1 (MQ=255)
  aaCGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAg                                                                                                                                                   <  2:579364/150‑1 (MQ=255)
   aCGAACATCCAGGTTCGGAAATTAATTTAATTGCGTTCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGt                                                                                                                                                  <  1:524569/150‑1 (MQ=255)
    cGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGTAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                 >  2:362201/1‑150 (MQ=255)
      aaCATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTAt                                                                                                                                               >  1:75588/1‑150 (MQ=255)
        cATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTg                                                                                                                                             >  2:1284003/1‑150 (MQ=255)
        cATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTg                                                                                                                                             >  2:1254468/1‑150 (MQ=255)
           ccAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGa                                                                                                                                          <  2:865035/150‑1 (MQ=255)
            cAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAt                                                                                                                                         >  2:1402611/1‑150 (MQ=255)
             aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAtt                                                                                                                                        <  1:505442/150‑1 (MQ=255)
             aGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAtt                                                                                                                                        <  1:1262540/150‑1 (MQ=255)
              ggTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAttt                                                                                                                                       <  1:200639/150‑1 (MQ=255)
               gTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGTAtttt                                                                                                                                      >  1:47350/1‑149 (MQ=255)
               gTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTg                                                                                                                                      >  2:622816/1‑150 (MQ=255)
                ttCGGATATTAATTTAATTGCGTGCATAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATTCTCGCAGTTATTTACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGATTTGt                                                                                                                                     <  1:347820/150‑1 (MQ=255)
                  cGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCt                                                                                                                                   >  2:1587374/1‑150 (MQ=255)
                   ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTg                                                                                                                                  >  1:478432/1‑150 (MQ=255)
                   ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTg                                                                                                                                  >  2:1024999/1‑150 (MQ=255)
                   ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTg                                                                                                                                  <  2:330531/150‑1 (MQ=255)
                   ggAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGATTTGttat                                                                                                                                  <  1:102822/150‑4 (MQ=255)
                    gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTgg                                                                                                                                 <  1:160132/150‑1 (MQ=255)
                    gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTgg                                                                                                                                 <  2:242056/150‑1 (MQ=255)
                    gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGATTTGttata                                                                                                                                 <  2:903476/150‑5 (MQ=255)
                     aaaTTAATTTAATTGAGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGt                                                                                                                                <  2:287844/150‑1 (MQ=255)
                        ttaatttaatTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGaa                                                                                                                             <  1:722771/150‑1 (MQ=255)
                        ttaatttaatTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGaa                                                                                                                             >  2:173774/1‑150 (MQ=255)
                         taatttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATCTGTCTGGTGAat                                                                                                                            <  2:1199612/150‑1 (MQ=255)
                            tttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattat                                                                                                                         <  2:1815096/150‑1 (MQ=255)
                                attGCGTGCAGAGAAAATAGTAATGAGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGt                                                                                                                     <  2:267269/150‑1 (MQ=255)
                                 ttGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTc                                                                                                                    <  2:1132737/150‑1 (MQ=255)
                                   gCGTGCAGAGAAAATAGCAATGCTCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                  <  1:248715/150‑1 (MQ=255)
                                   gCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                  <  2:826297/150‑1 (MQ=255)
                                   gCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                  >  2:527048/1‑150 (MQ=255)
                                     gTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                <  1:1583571/150‑1 (MQ=255)
                                     gTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                >  1:858348/1‑150 (MQ=255)
                                         agagAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCtt                                                                                                            <  1:1664859/150‑1 (MQ=255)
                                          gagaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCttt                                                                                                           >  2:1825377/1‑150 (MQ=255)
                                            taaaaTATGAATTCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAAATTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTcc                                                                                                         <  1:145615/149‑1 (MQ=255)
                                               aaTAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCg                                                                                                      >  2:1330248/1‑150 (MQ=255)
                                                  aGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCca                                                                                                   >  2:407908/1‑150 (MQ=255)
                                                        gcgcTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGtt                                                                                             <  2:1168284/150‑1 (MQ=255)
                                                           cTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTAt                                                                                          <  1:499902/150‑1 (MQ=255)
                                                                gATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGc                                                                                     <  2:885866/150‑1 (MQ=255)
                                                                gATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCACGCCAGCAGTTTATCCAGc                                                                                     >  2:358601/1‑150 (MQ=255)
                                                                 aTAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCt                                                                                    >  1:1060055/1‑150 (MQ=255)
                                                                 aTAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCt                                                                                    <  1:1144413/150‑1 (MQ=255)
                                                                 aTAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCt                                                                                    <  2:679224/150‑1 (MQ=255)
                                                                    aaaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCAt                                                                                 >  1:393225/1‑150 (MQ=255)
                                                                     aaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATc                                                                                >  2:949769/1‑150 (MQ=255)
                                                                      aaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCAGCCAGCAGTTTATCCAGCTCATCg                                                                               >  2:658479/1‑150 (MQ=255)
                                                                      aaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTAAGGAATTTATAAATTTATATATAACAAATCCCAATAATTAATTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGATCAACg                                                                               >  2:84147/1‑150 (MQ=255)
                                                                        tatGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGcc                                                                             <  2:478432/150‑1 (MQ=255)
                                                                           gACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCAcc                                                                          <  1:300297/150‑1 (MQ=255)
                                                                               aTGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTGTATCCAGCTCATCGCCACCGACg                                                                      <  2:1326799/150‑1 (MQ=255)
                                                                               aTGCTAGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAATTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACg                                                                      <  1:1756182/150‑1 (MQ=255)
                                                                                   tCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCACAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTACCCGCCAGCAGTGTATACAGCTCATCGCCACCGACGTGAc                                                                  >  2:477517/1‑150 (MQ=255)
                                                                                     gCAGTTATTTACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCACTCCAGCAGTTTATCCAGCTCATCGCCACCGCCCAGACgg                                                                >  2:995504/1‑150 (MQ=255)
                                                                                     gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCAAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCACGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACgg                                                                >  1:903374/1‑150 (MQ=255)
                                                                                      cAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGa                                                               <  2:1721751/150‑1 (MQ=255)
                                                                                      cAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAAATAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTATCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGa                                                               <  2:988948/150‑1 (MQ=255)
                                                                                        gTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGaaa                                                             <  1:1304559/150‑1 (MQ=255)
                                                                                          tattaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATc                                                           >  2:565468/1‑150 (MQ=255)
                                                                                          tattaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATc                                                           >  2:929456/1‑150 (MQ=255)
                                                                                            ttaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCt                                                         >  1:668311/1‑150 (MQ=255)
                                                                                            ttaACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCt                                                         <  1:703079/150‑1 (MQ=255)
                                                                                              aaCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGc                                                       <  1:1699618/150‑1 (MQ=255)
                                                                                                  ttGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCt                                                   <  1:975486/150‑1 (MQ=255)
                                                                                                     tgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCa                                                <  1:628382/150‑1 (MQ=255)
                                                                                                            ttAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTACCCGCCAGCAGTTTATCCAGCTCATCGCAACCGACGTGACGGAAATCCTGCCCCTTCACgtagta                                         >  1:992949/1‑150 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCTACAGCAGTTTATCCAGCTCATAGCCAACGACGTGACGGAAATCATGCCCCTTCACGTAGTataa                                      >  1:304491/1‑150 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                      <  1:929456/150‑1 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                      >  2:1286391/1‑150 (MQ=255)
                                                                                                                gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCTCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaag                                     >  1:80727/1‑150 (MQ=255)
                                                                                                                gAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaag                                     <  2:1735461/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaaga                                    <  1:1052817/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaaga                                    <  1:622816/150‑1 (MQ=255)
                                                                                                                             tatatataACAAATCCCAATTATTAAGTTATTGGGATTTGTCTGGTGATTTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACaaa                        >  1:186787/1‑150 (MQ=255)
                                                                                                                               tatataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAata                      >  1:844150/1‑150 (MQ=255)
                                                                                                                                atataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAatat                     >  1:854190/1‑150 (MQ=255)
                                                                                                                                 tataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                    >  1:1692811/1‑150 (MQ=255)
                                                                                                                                 tataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                    >  1:1040741/1‑150 (MQ=255)
                                                                                                                                 tataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                    <  1:407908/150‑1 (MQ=255)
                                                                                                                                 tataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCACCGCCAGCAGTTTATCCAGCTCTTCGCCAACGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATAtt                    >  1:1752821/1‑150 (MQ=255)
                                                                                                                                  ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                   <  2:383927/150‑1 (MQ=255)
                                                                                                                                   taACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                  <  1:1216329/150‑1 (MQ=255)
                                                                                                                                   agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                  >  2:1289721/3‑150 (MQ=255)
                                                                                                                                   agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCAAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                  >  1:1222885/3‑150 (MQ=255)
                                                                                                                                    aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                 <  1:1402611/150‑1 (MQ=255)
                                                                                                                                    aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                 <  1:273368/150‑1 (MQ=255)
                                                                                                                                     aCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTgg                >  2:1669760/1‑150 (MQ=255)
                                                                                                                                        aaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg             >  1:863672/1‑150 (MQ=255)
                                                                                                                                        aaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGCGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg             >  1:1423141/1‑150 (MQ=255)
                                                                                                                                        aaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTAGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg             <  2:75588/150‑1 (MQ=255)
                                                                                                                                        aaTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg             >  2:1179044/1‑150 (MQ=255)
                                                                                                                                           cccAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCgg          >  2:688942/1‑150 (MQ=255)
                                                                                                                                           cccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTATTAGAAGATAAACTCACAAATATTATTGCAGccg          >  2:754706/1‑148 (MQ=255)
                                                                                                                                           cccAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCgg          >  2:1251848/1‑150 (MQ=255)
                                                                                                                                             cAATAATTAAGTTATTGGGCTTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc        >  1:1417690/1‑150 (MQ=255)
                                                                                                                                              aataatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAACCTCACAAATATTCTGGCAGCGGTCg       <  1:1284003/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAAT  >  NC_000913/3906679‑3906975

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: