Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,858,331 T→A 6.7% K25I (AAA→ATA)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,858,3310TA 6.7% 280.7 / 6.8 112K25I (AAA→ATA) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (48/53);  new base A (2/6);  total (51/61)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 113 total aligned reads displayed.

TTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGT  >  NC_000913/3858184‑3858477
                                                                                                                                                   |                                                                                                                                                  
ttgttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGttt                                                                                                                                                   >  2:531285/1‑149 (MQ=255)
 tgttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGGACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTa                                                                                                                                                  >  1:834225/1‑149 (MQ=255)
   ttgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGATTTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTAcc                                                                                                                                                >  1:472857/1‑149 (MQ=255)
    tgGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAATTAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATATTGCGTAGTCAGCGGGTTCCCg                                                                                                                                               >  2:927440/1‑149 (MQ=255)
     ggTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTTAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGc                                                                                                                                              >  1:310408/1‑149 (MQ=255)
     ggTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCTATATTTTGCGTATTCAGCGGTTTACCGc                                                                                                                                              >  1:270817/1‑149 (MQ=255)
      gTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCCCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCt                                                                                                                                             <  1:380882/149‑1 (MQ=255)
      gTCCACGCGCCGGACTGGTTAGCGGAGATACCGGTAAATAGTCAGGCCCGTGCAGGCGTACAAACCGGTGACCACGGCTCGGCGGAAATTAAGCGACACCCATCCGCCGCCATACTATCCATCTTTTGCGTATTCAGCGGGTTCCCGCt                                                                                                                                             >  2:757025/1‑149 (MQ=255)
       tCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTTCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGATTGAAAGCGAATACCTTCCGCCGCCAGACTATCATTATTTTGCGTATTCAGCGGTTTAACGcgt                                                                                                                                            >  2:522007/1‑147 (MQ=255)
        ccACGGGCCGGACTGGTTCTCGTAGATAACGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTATGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTAAACCGCtga                                                                                                                                           <  2:911254/149‑3 (MQ=255)
        ccACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCtat                                                                                                                                           <  2:310408/149‑1 (MQ=255)
         cACGGGCCGGACTGGTTCGCGTAGATACCGGTAAATAGACCGTCCAGTGCAGGCGTACAACCATTTGAACAGGTGTAGGCGTAATTAAAGCGAATACCTTCCTCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTAAACCGCtgaa                                                                                                                                          <  1:595079/149‑4 (MQ=255)
         cACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTAAACCGCtgaa                                                                                                                                          <  1:1018728/149‑4 (MQ=255)
         cACGGGCCGGACTGGGTAGCGCCGAGACCGGTAAATAGCCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATCTTTTGCGTATTCAGCGGTAAACCGcggaa                                                                                                                                          <  1:591876/149‑5 (MQ=255)
          aCGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAg                                                                                                                                         >  2:229196/1‑149 (MQ=255)
            gggCCGGACTGGTTAGCGTAGCTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCa                                                                                                                                       <  2:95586/149‑1 (MQ=255)
               ccGGTCTGGCTAGCGTACATTCCGGTAAATAGTCCGGCCCGTGCAGGCGTACCAACCCTGTAACAGGTTTAGGCGTAATTAAAGCGAATACATTCCGTCGCCAGACTATCAATATTTTGAGTATTCAGCGGTTTACCGCTATAGCAAcc                                                                                                                                    <  2:675522/149‑1 (MQ=255)
                 ggACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTAAACCGCtgaatac                                                                                                                                       >  2:384146/1‑138 (MQ=255)
                   aCTGGTTAGCGTCGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGAGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCACACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAcc                                                                                                                                <  1:771975/149‑1 (MQ=255)
                   aCTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGAcc                                                                                                                                >  2:520506/1‑149 (MQ=255)
                    cTGGTTAGCTTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGAGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCa                                                                                                                               <  2:469655/149‑1 (MQ=255)
                       ggtAGCGTAGAGACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGCGAACAGGTGTAGGCGGTATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCatat                                                                                                                            <  1:189869/147‑1 (MQ=255)
                       gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCCGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAAAACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCatat                                                                                                                            >  2:49766/1‑149 (MQ=255)
                       gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCatat                                                                                                                            >  2:39394/1‑149 (MQ=255)
                       gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCatat                                                                                                                            >  2:866511/1‑149 (MQ=255)
                       gTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCatat                                                                                                                            >  1:290802/1‑149 (MQ=255)
                           gCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggt                                                                                                                        <  2:129234/149‑1 (MQ=255)
                           gCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggt                                                                                                                        <  1:443632/149‑1 (MQ=255)
                           gCGTAGATACCGGTAAATAGACCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTACAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggt                                                                                                                        <  2:75488/149‑1 (MQ=255)
                            ggTAGATAGCGGTATATAGTGCGGCGCGTGCAGGCGTACTAACCGGTGAACAGGTGTAGGCTGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCATCCGACCATATtggtg                                                                                                                       <  1:215572/148‑1 (MQ=255)
                             gTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGAGTATTCAGCGGGTTACCGCTATAGCAACCGACCATATtggtgg                                                                                                                      >  1:869982/1‑149 (MQ=255)
                             gTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGCGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATtggtgg                                                                                                                      <  2:32504/149‑1 (MQ=255)
                                gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGCCTAGAAATATTTTGCGTATTAAGCGATTTACAGCTATTGCAACCGACCTTATTGGTGGcgg                                                                                                                   >  1:122248/1‑147 (MQ=255)
                                gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAGCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt                                                                                                                   <  1:770111/149‑1 (MQ=255)
                                 aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAAAACACCATATTGGTGGCCTg                                                                                                                  >  2:207940/1‑149 (MQ=255)
                                  ttccGGTAAATAGTCCGGCCCGTGCAGGCGTACATACCGGTGAACAGGTGTAGGCGGAATTAAAGCGACTACCTTCCGCCGCCAGACTATTAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTgg                                                                                                                 <  2:611376/147‑1 (MQ=255)
                                  tACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTgg                                                                                                                 <  1:447411/149‑1 (MQ=255)
                                   aCCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATGTTTTGCGTATTCAGCGGTTTTCCGCTATAGCAACCGACCATTTTGGTTGCCTggg                                                                                                                >  2:777646/1‑149 (MQ=255)
                                        tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCg                                                                                                           >  2:1573/1‑149 (MQ=255)
                                            tAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCa                                                                                                       >  1:630935/1‑149 (MQ=255)
                                              gTCCGGCCCGTGCAGGCGTACAAACCGCTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCTCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATg                                                                                                     <  2:449983/149‑1 (MQ=255)
                                               tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATTCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGa                                                                                                    >  2:136834/1‑149 (MQ=255)
                                               tCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGa                                                                                                    >  2:865014/1‑149 (MQ=255)
                                               gccGGCCCGTGCACGCGTACAAACCGGTGAACAGGTGGAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTATACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGa                                                                                                    <  1:250229/148‑1 (MQ=255)
                                                   gCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                <  1:219590/149‑1 (MQ=255)
                                                   gCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                <  1:866511/149‑1 (MQ=255)
                                                    cccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAAc                                                                                               <  2:293562/149‑1 (MQ=255)
                                                    cccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAAc                                                                                               <  1:229196/149‑1 (MQ=255)
                                                     ccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGCATCCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGcaca                                                                                              >  1:952731/1‑149 (MQ=255)
                                                      cGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCAAATTGGTGGCCGGGGTATCGGTCATGACGAACAg                                                                                             >  1:150196/1‑149 (MQ=255)
                                                           aGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAt                                                                                        >  2:406559/1‑149 (MQ=255)
                                                             gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTg                                                                                      <  2:550670/149‑1 (MQ=255)
                                                             gCGAACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTg                                                                                      <  1:176424/149‑1 (MQ=255)
                                                              cGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTgg                                                                                     <  1:541286/149‑1 (MQ=255)
                                                              cGTACAAACCGGTGAACAGGTGTAGGAGGAATTAAAGCGAAAACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTgg                                                                                     <  2:834225/149‑1 (MQ=255)
                                                                tACAAACCTGTGAACAGGTGTAGTCGGAATTAAAGAGAATTCCTTCCGCCGCCAGACTCTCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGGCTGGGTATCGGTCATGACGAACAGAAAATTgggg                                                                                   <  2:854720/149‑1 (MQ=255)
                                                                   aaaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCCGTTTACCGCTATAGCAACCGACCATATTCGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGt                                                                                >  1:298353/1‑149 (MQ=255)
                                                                    aaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCTTATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGtt                                                                               <  1:972798/149‑1 (MQ=255)
                                                                    aaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTAGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGtt                                                                               <  1:504886/149‑1 (MQ=255)
                                                                       cGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGGATTCCGCGGTTTACCGCTATAGCAACCGACGATATTGGGGGCCTGGGTAGCGGACATGACGAACAGAAAATTGGGGCGTTTCa                                                                            <  1:749838/149‑1 (MQ=255)
                                                                       cGGTGAACAGGTGTAGGCGGAATTAAAGAGAATACCTACCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCTCTATAGCAACCGACCATATTGGTGTCCTGGGTATCGGTCATGAAGAAAAGAAAATTGGGGCGTTTCa                                                                            >  1:213728/1‑149 (MQ=255)
                                                                         gTGAACAGGTGTAGTCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCTACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATa                                                                          <  1:64596/149‑1 (MQ=255)
                                                                           gAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGGCATGACGAACAGAAAATTGGGGCGTTTCATAgg                                                                        <  2:69196/149‑1 (MQ=255)
                                                                             aCAGGTGTACGCGGAATTAAAGCGAATACCATCCCCCGCCAGACTATCAAAATTTTGCTTATTAATCGGATTACCCCTAAAGCAACCGACCATATTGGTGCCCTGGGTATCGGTCATCACGAACAGAAAATTGGGGCGTTTCATAGGCt                                                                      >  2:16336/1‑149 (MQ=255)
                                                                              cAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCATTATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTCGCCTGGGTATCGGTCATGACGAACAGAAAATTGGCTCGTTTGATAGGCTa                                                                     >  2:321200/1‑149 (MQ=255)
                                                                                 gtgtAGGCGGAATTAAAGCGAATACCTTACGCCGCCAGACTCTCAATATTTTGCGTATTCAGCGGATTACCGCTATAGCAACCGACCATATTGGGGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATcc                                                                  <  1:151392/149‑1 (MQ=255)
                                                                                    tAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                               <  1:816824/149‑1 (MQ=255)
                                                                                    tAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTGTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                               <  1:933293/149‑1 (MQ=255)
                                                                                      ggCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCATGGGTATCGGTCATGACGAACAGATAATTGGGGCGGTTCATAGGCTATCCTTccc                                                             >  2:427591/1‑149 (MQ=255)
                                                                                           aaTTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCgg                                                        <  2:290802/149‑1 (MQ=255)
                                                                                             ttAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGa                                                      <  1:379161/149‑1 (MQ=255)
                                                                                             ttAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGa                                                      >  1:349374/1‑149 (MQ=255)
                                                                                                          cTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                         >  2:178156/1‑149 (MQ=255)
                                                                                                          cTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                         >  1:268741/1‑149 (MQ=255)
                                                                                                          cTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACgg                                         <  2:216676/149‑1 (MQ=255)
                                                                                                           ttCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGc                                        >  1:863539/1‑149 (MQ=255)
                                                                                                            tCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCt                                       <  2:7908/149‑1 (MQ=255)
                                                                                                              cgccgcCAGAGTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATAATGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTgg                                     <  1:6466/149‑1 (MQ=255)
                                                                                                                   ccAGACTATCAATATTTTGCGTATTCAGCGGTGAACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTATCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGt                                <  2:49368/149‑1 (MQ=255)
                                                                                                                      gACTATCAATAATTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCAGGACGAACAGAAAATTGGGGCGGGTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGAc                             <  1:966283/149‑1 (MQ=255)
                                                                                                                         tATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGAt                          >  1:160934/1‑149 (MQ=255)
                                                                                                                           tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgc                        >  2:6801/1‑149 (MQ=255)
                                                                                                                             cttatTTTCCGTATTCCGCGGTTTACCGCTATACCAACCGACCATATTTGTGGCCTGTGTATCGGTCCTGTCGAACAGAAAATTGGGGCGTTTCATAGGCTTTCGTTCCGGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccg                      <  2:708980/147‑1 (MQ=255)
                                                                                                                             aaTATTTTGCGTATTCAGCGGTTTACCGCTATAGCCACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTCCGGCTGGGTCGTGACGATgccg                      >  1:544104/1‑149 (MQ=255)
                                                                                                                             aaTATTTTGCGTATTCAGCGGTTTACCGCTATAGCACCCGACCATATTGGTGGCCTGGGTATCGGGCATGCCGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccg                      >  2:412124/1‑149 (MQ=255)
                                                                                                                             aaTATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGATACGGCAGGCTCGTGACGATgccg                      >  2:172742/1‑149 (MQ=255)
                                                                                                                              atatTTTGCGTATTCCGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCGTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccgc                     >  1:779540/1‑149 (MQ=255)
                                                                                                                                   ttGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCt                <  2:243116/149‑1 (MQ=255)
                                                                                                                                   ttGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGAGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCt                <  1:49766/149‑1 (MQ=255)
                                                                                                                                    tGCGTATTCTGCTGTTTACCGCTATAGCAACCGACCAAAATGGTGGCCGGGGTATCGGTCCTGACGAACAGAAAATTGGGGCGTTTGATAGGCTATCCTTCCCGTCGGGTGTAAACCTGAGGGCTGGGTCGTGACGATGCCGCCAGCTg               <  1:576197/149‑1 (MQ=255)
                                                                                                                                      cGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTc             <  1:319779/149‑1 (MQ=255)
                                                                                                                                       gTAGTCAGCGGTTTACCGCTATAGCCACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGGCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCa            <  1:139700/149‑1 (MQ=255)
                                                                                                                                          ttCAGGGGTGTAGCGGTATAGCAACCGACCATATTGGTGGCCTGGGTATGGGGGAAGACCAACAGAAAATTGGGGAGTTGGAGAGGCTGTCATTCCCGTCGGGAGGATAGGGTACGGCTGGGGCGTCTCGATGCGACCAGCTGGGaaaa         >  2:41399/1‑149 (MQ=255)
                                                                                                                                          ttCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCaaaa         <  1:1573/149‑1 (MQ=255)
                                                                                                                                            cAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAAt       >  1:384304/1‑149 (MQ=255)
                                                                                                                                             aGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATc      >  1:402511/1‑149 (MQ=255)
                                                                                                                                             aGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATc      >  1:833604/1‑149 (MQ=255)
                                                                                                                                              ccGGTTTACCGCTATAGCAACCGTCCATATTGGTGGCCTGGGTATCGGAGATGACGTACAGAAAATTGGGGCGTTTCATAGGCTTTCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCa     <  1:321200/148‑1 (MQ=255)
                                                                                                                                               cGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAg    <  2:349374/149‑1 (MQ=255)
                                                                                                                                                 gTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGt  <  1:804604/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGT  >  NC_000913/3858184‑3858477

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: