Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 242,344 C→G 5.7% P320P (CCG→CCC fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913242,3440CG 5.7% 615.7 / 20.0 176P320P (CCG→CCCfadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (74/92);  new base G (6/4);  total (80/96)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 183 total aligned reads displayed.

GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242196‑242491
                                                                                                                                                    |                                                                                                                                                   
ggCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGccc                                                                                                                                                     <  3:884550/149‑1 (MQ=255)
 gCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCg                                                                                                                                                    >  3:86782/1‑149 (MQ=255)
 gCACTCCACCAACATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACGGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGAGGCGACGACCAATTTCCACGCCCg                                                                                                                                                    >  3:1035453/1‑149 (MQ=255)
  cACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCgg                                                                                                                                                   <  2:886432/149‑1 (MQ=255)
   aCTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGc                                                                                                                                                  >  3:1376953/1‑149 (MQ=255)
    cTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCg                                                                                                                                                 >  3:391434/1‑149 (MQ=255)
     tCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                                >  1:498074/1‑149 (MQ=255)
      ccaccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtg                                                                                                                                               <  3:636473/149‑1 (MQ=255)
      ccaccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtg                                                                                                                                               <  4:942003/149‑1 (MQ=255)
        accaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtggt                                                                                                                                             >  1:1458056/1‑149 (MQ=255)
         ccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtggtg                                                                                                                                            >  1:765676/1‑149 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  4:1376953/149‑4 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGtggtgg                                                                                                                                           >  3:33383/1‑149 (MQ=255)
           aGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTggtggt                                                                                                                                          <  1:1253980/149‑1 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         >  2:1358022/1‑149 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         <  1:414397/149‑1 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         >  1:70381/1‑149 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                         <  2:765676/149‑1 (MQ=255)
             cATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        <  2:1224943/149‑1 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaattgg                                                                                                                                       <  3:1145006/149‑8 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTgg                                                                                                                                       <  2:251367/149‑1 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTgg                                                                                                                                       <  3:1398908/149‑1 (MQ=255)
               tCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTggg                                                                                                                                      >  3:203271/1‑149 (MQ=255)
                ccGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     <  4:1227628/149‑1 (MQ=255)
                ccGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     <  4:362978/149‑1 (MQ=255)
                ccGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTAGGa                                                                                                                                     <  2:132360/149‑1 (MQ=255)
                   ccAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCa                                                                                                                                  >  1:226344/1‑149 (MQ=255)
                   ccAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCa                                                                                                                                  <  2:247828/149‑1 (MQ=255)
                   ccAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCa                                                                                                                                  <  2:1375069/149‑1 (MQ=255)
                        cTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGca                                                                                                                             >  4:1168986/1‑149 (MQ=255)
                        cTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGca                                                                                                                             >  2:1426021/1‑149 (MQ=255)
                        cTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGca                                                                                                                             >  2:1055895/1‑149 (MQ=255)
                           gCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAg                                                                                                                          <  2:434169/149‑1 (MQ=255)
                             ccGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGt                                                                                                                        >  1:1033158/1‑149 (MQ=255)
                               ggCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                                      <  4:426615/149‑1 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                                     >  1:1248719/1‑149 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                                     <  1:1075173/149‑1 (MQ=255)
                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   >  1:581868/1‑149 (MQ=255)
                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   <  1:800880/149‑1 (MQ=255)
                                   aTTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                                  >  2:649343/1‑149 (MQ=255)
                                   aTTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                                  >  4:214217/1‑149 (MQ=255)
                                   aTTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGcc                                                                                                                  >  1:725233/1‑149 (MQ=255)
                                    ttttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCt                                                                                                                 >  4:17678/1‑149 (MQ=255)
                                       tCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTaaa                                                                                                              <  2:1458056/149‑1 (MQ=255)
                                        cGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAt                                                                                                             <  4:86782/149‑1 (MQ=255)
                                           cccGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtctt                                                                                                          >  4:553517/1‑149 (MQ=255)
                                                     atgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgcc                                                                                                >  2:1557239/1‑149 (MQ=255)
                                                     atgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgcc                                                                                                <  4:1425409/149‑1 (MQ=255)
                                                      tgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     <  3:55570/111‑1 (MQ=255)
                                                      tgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     >  4:55570/1‑111 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                              <  2:584660/149‑1 (MQ=255)
                                                          atcgatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTa                                                                                           <  3:1168986/149‑1 (MQ=255)
                                                          atcgatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTa                                                                                           >  2:891285/1‑149 (MQ=255)
                                                             gatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAACTCTTCTGCACCGCCGCGTAAtt                                                                                        <  4:1304761/149‑1 (MQ=255)
                                                              atcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   >  2:1442687/1‑121 (MQ=255)
                                                              atcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAttt                                                                                       <  2:1017308/149‑1 (MQ=255)
                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  4:38419/149‑1 (MQ=255)
                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  3:1488754/149‑1 (MQ=255)
                                                                    aCGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCgg                                                                                 <  1:466906/149‑1 (MQ=255)
                                                                      gAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCTGACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGt                                                                               >  2:1252500/1‑149 (MQ=255)
                                                                      gAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGt                                                                               <  2:104896/149‑1 (MQ=255)
                                                                       aaGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTc                                                                              >  1:608399/1‑149 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         >  1:531165/1‑149 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         >  3:767988/1‑149 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGCGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         <  1:968093/149‑1 (MQ=255)
                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       <  4:1107065/149‑1 (MQ=255)
                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATAAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       <  2:509309/149‑1 (MQ=255)
                                                                                ttACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGttt                                                                     <  4:1389184/149‑1 (MQ=255)
                                                                                             gTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCTTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCa                                                        >  4:318346/1‑149 (MQ=255)
                                                                                             gTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  4:1240934/122‑1 (MQ=255)
                                                                                                 gTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGAATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCa                                                    >  4:460716/1‑149 (MQ=255)
                                                                                                  ttCTGGAACGGTACGTTCAGCGGAAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCGAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCAc                                                   <  2:231528/149‑1 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCATCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                               >  2:690580/1‑149 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                               <  3:429729/149‑1 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                               <  2:71455/149‑1 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                               <  1:649343/149‑1 (MQ=255)
                                                                                                        aaCGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         >  1:1049397/1‑121 (MQ=255)
                                                                                                                gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTg                                     >  1:775410/1‑149 (MQ=255)
                                                                                                                gTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTg                                     >  1:672736/1‑149 (MQ=255)
                                                                                                                 ttCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTgt                                    <  4:1030268/149‑2 (MQ=255)
                                                                                                                          aaGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAAt                           <  3:1529705/149‑1 (MQ=255)
                                                                                                                             tGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTa                        <  1:1426021/149‑1 (MQ=255)
                                                                                                                               gCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       <  2:549883/100‑1 (MQ=255)
                                                                                                                               gCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       >  1:549883/1‑100 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc              >  3:598387/1‑149 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc              >  4:287673/1‑149 (MQ=255)
                                                                                                                                        aaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAAt                           <  4:1251273/130‑1 (MQ=255)
                                                                                                                                        aaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAAt                           >  3:1251273/6‑135 (MQ=255)
                                                                                                                                        aaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc             >  1:1387051/6‑149 (MQ=255)
                                                                                                                                        aaTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc             <  1:1441437/149‑1 (MQ=255)
                                                                                                                                          tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCTTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  3:318346/149‑1 (MQ=255)
                                                                                                                                          tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  1:1003264/149‑1 (MQ=255)
                                                                                                                                          ccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           >  3:264481/4‑149 (MQ=255)
                                                                                                                                           caccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                >  1:249246/3‑143 (MQ=255)
                                                                                                                                            accACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGTCAGCGTAATGTAGCGTTTGTTCCAGGt         >  4:70719/2‑149 (MQ=255)
                                                                                                                                              cACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCa       >  4:613570/1‑149 (MQ=255)
                                                                                                                                               aCGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAg      >  4:1442817/1‑149 (MQ=255)
                                                                                                                                                  cccGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACg   >  1:174760/1‑149 (MQ=255)
                                                                                                                                                   ccGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGc  >  1:1190469/1‑149 (MQ=255)
                                                                                                                                                    cGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  3:958355/139‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242196‑242491

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