Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 300,692 A→C 5.0% I81M (ATT→ATG paoC ← PaoABC aldehyde oxidoreductase, Moco‑containing subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,6920AC 5.0% 985.4 / 28.8 271I81M (ATT→ATGpaoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (112/144);  new base C (12/2);  total (124/147)
Fisher's exact test for biased strand distribution p-value = 4.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

TGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGT  >  NC_000913/300543‑300841
                                                                                                                                                     |                                                                                                                                                     
tGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTa                                                                                                                                                       >  1:2821061/1‑150 (MQ=255)
 gCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTaa                                                                                                                                                      <  2:2556691/150‑1 (MQ=255)
   aCCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATg                                                                                                                                                    >  2:4354319/1‑150 (MQ=255)
   aCCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATg                                                                                                                                                    >  1:83975/1‑150 (MQ=255)
    ccAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGa                                                                                                                                                   <  1:93762/150‑1 (MQ=255)
    ccAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGa                                                                                                                                                   >  1:3169161/1‑150 (MQ=255)
     cAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAc                                                                                                                                                  >  1:4531094/1‑150 (MQ=255)
           ggCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAg                                                                                                                                            <  2:1829193/150‑1 (MQ=255)
            gccgccGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGt                                                                                                                                           >  1:1178996/1‑150 (MQ=255)
            gccgccGCTCGCGCCGGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGt                                                                                                                                           >  1:1532178/1‑150 (MQ=255)
             ccgccgCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTa                                                                                                                                          <  1:4192124/150‑1 (MQ=255)
             ccgccgCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTa                                                                                                                                          <  1:2517387/150‑1 (MQ=255)
             ccgccgCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTa                                                                                                                                          >  2:33234/1‑150 (MQ=255)
                 cgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGcc                                                                                                                                      <  2:1856360/150‑1 (MQ=255)
                 cgcTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGcc                                                                                                                                      <  1:3639533/150‑1 (MQ=255)
                   cTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCg                                                                                                                                    <  2:2732053/150‑1 (MQ=255)
                    tCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCgg                                                                                                                                   <  1:4754676/150‑1 (MQ=255)
                           tGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                                                            >  2:1956011/1‑150 (MQ=255)
                           tGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCtttt                                                                                                                            >  2:1072638/1‑150 (MQ=255)
                              tCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgc                                                                                                                         >  1:2261227/1‑150 (MQ=255)
                               cGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATTGTGGGGCCGCCTAACAGCCTGGCGGTGTGTGTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcg                                                                                                                        >  1:784385/1‑150 (MQ=255)
                                gAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcgc                                                                                                                       >  1:4421576/1‑150 (MQ=255)
                                gAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcgc                                                                                                                       >  1:1756464/1‑150 (MQ=255)
                                 aaGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGcgcg                                                                                                                      <  1:749488/150‑1 (MQ=255)
                                   ggTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggc                                                                                                                    >  2:218858/1‑150 (MQ=255)
                                    gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCgcggcg                                                                                                                   <  1:1416222/150‑1 (MQ=255)
                                      ctcGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTc                                                                                                                 >  2:2454043/1‑150 (MQ=255)
                                        cGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCg                                                                                                               >  1:1292043/1‑150 (MQ=255)
                                          gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTa                                                                                                             >  2:3081374/1‑150 (MQ=255)
                                           ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTAt                                                                                                            >  1:4101857/1‑150 (MQ=255)
                                            cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATc                                                                                                           <  2:492391/150‑1 (MQ=255)
                                              cTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCaa                                                                                                         <  1:457515/150‑1 (MQ=255)
                                                   aGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCg                                                                                                    <  1:1954340/150‑1 (MQ=255)
                                                   aGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCg                                                                                                    <  1:4411016/150‑1 (MQ=255)
                                                       cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGa                                                                                                <  1:2889678/150‑1 (MQ=255)
                                                             ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTc                                                                                          <  1:3145563/150‑1 (MQ=255)
                                                             ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTc                                                                                          <  2:2975033/150‑1 (MQ=255)
                                                             ccTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTc                                                                                          <  1:3831207/150‑1 (MQ=255)
                                                             ccTGATGAGAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTc                                                                                          <  2:4403429/150‑1 (MQ=255)
                                                              cTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGAATCAAGGGCGGTGAGGCGTcc                                                                                         >  1:1088504/1‑150 (MQ=255)
                                                               tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCt                                                                                        >  2:1618413/1‑150 (MQ=255)
                                                                  tgatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTg                                                                                     <  2:2811809/150‑1 (MQ=255)
                                                                    atAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGc                                                                                   <  1:2192000/150‑1 (MQ=255)
                                                                        tGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATg                                                                               >  1:2674314/1‑150 (MQ=255)
                                                                         gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTAGTGGCGGTAATGACAGCCTGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                              <  2:3188697/150‑1 (MQ=255)
                                                                         gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATgg                                                                              <  2:1369685/150‑1 (MQ=255)
                                                                          cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc                                                                             <  1:888919/150‑1 (MQ=255)
                                                                          cTCAATAGAGGGGCCGCCTAACAGCCAGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTAGCGCCGCGTCCGTATCAACTGCGGTGAGGCGTCCaatggcaatggc                                                                             >  1:4660099/1‑150 (MQ=255)
                                                                          cTCAAGAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGc                                                                             <  2:2777707/150‑1 (MQ=255)
                                                                           tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                            <  1:4112980/150‑1 (MQ=255)
                                                                           tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                            <  1:4062949/150‑1 (MQ=255)
                                                                           tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCg                                                                            <  1:4453986/150‑1 (MQ=255)
                                                                           tCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTGTTTGCGCGGCGTCCGTATCACGGGCGGTGAGGCGTCCTTTGGCACTGGCg                                                                            >  2:3050557/1‑150 (MQ=255)
                                                                            cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGACCCGGCGTTACTGGCGTTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCgg                                                                           >  2:356313/1‑150 (MQ=255)
                                                                             aaTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCAGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGa                                                                          <  2:2846088/150‑1 (MQ=255)
                                                                              aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGaa                                                                         <  1:669663/150‑1 (MQ=255)
                                                                              aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGaa                                                                         >  1:2795160/1‑150 (MQ=255)
                                                                              aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGaa                                                                         <  1:2976656/150‑1 (MQ=255)
                                                                               tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                        <  1:1003433/150‑1 (MQ=255)
                                                                               tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAc                                                                        <  1:2328525/150‑1 (MQ=255)
                                                                                  tGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcga                                                                     >  2:3739806/1‑150 (MQ=255)
                                                                                    gggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacg                                                                   <  1:4792573/150‑1 (MQ=255)
                                                                                    gggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACcgacg                                                                   >  1:4572991/1‑150 (MQ=255)
                                                                                       ccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGata                                                                >  2:1633480/1‑150 (MQ=255)
                                                                                       ccgccTAACAGCCTGGCGGTGTTTTTGTAGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGata                                                                >  2:82348/1‑150 (MQ=255)
                                                                                        cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCGTTGGCAATGGCGGCACCGACGatat                                                               >  1:1711338/1‑150 (MQ=255)
                                                                                          ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAg                                                             <  2:2919936/150‑1 (MQ=255)
                                                                                          ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAg                                                             <  2:4580227/150‑1 (MQ=255)
                                                                                                   cTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATTGCGGAACCGACGATATAGCCATAggcg                                                    >  1:3059851/1‑150 (MQ=255)
                                                                                                    tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAggcgg                                                   <  2:2681005/150‑1 (MQ=255)
                                                                                                       cGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGt                                                >  2:4675684/1‑150 (MQ=255)
                                                                                                           gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTggg                                            <  2:150017/150‑1 (MQ=255)
                                                                                                           gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTggg                                            <  2:1808397/150‑1 (MQ=255)
                                                                                                            tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTgggg                                           >  1:3480176/1‑150 (MQ=255)
                                                                                                              tttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGc                                         >  1:4117820/1‑150 (MQ=255)
                                                                                                                      cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCAt                                 <  1:2862116/150‑1 (MQ=255)
                                                                                                                      cTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCAt                                 <  2:4445395/150‑1 (MQ=255)
                                                                                                                       tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATg                                <  2:3275306/150‑1 (MQ=255)
                                                                                                                        ttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGc                               <  1:485783/150‑1 (MQ=255)
                                                                                                                        ttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGc                               <  1:2036492/150‑1 (MQ=255)
                                                                                                                             gAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTc                          >  2:2436509/1‑150 (MQ=255)
                                                                                                                                 gCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAg                      <  1:3920464/150‑1 (MQ=255)
                                                                                                                                  ccccGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAgg                     <  1:1267104/150‑1 (MQ=255)
                                                                                                                                      ggCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTa                 <  1:3666639/150‑1 (MQ=255)
                                                                                                                                       ggcgtACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                >  2:660519/5‑150 (MQ=255)
                                                                                                                                       ggcgtACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                >  2:2220059/5‑150 (MQ=255)
                                                                                                                                       gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                <  1:2546830/150‑1 (MQ=255)
                                                                                                                                        cGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGc               <  1:2933787/150‑1 (MQ=255)
                                                                                                                                         gTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcg              <  1:4290947/150‑1 (MQ=255)
                                                                                                                                          ttACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgt             >  1:3376991/1‑150 (MQ=255)
                                                                                                                                          gtACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgt             >  2:595639/2‑150 (MQ=255)
                                                                                                                                          gtACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTAGTAGGCGTAGcgt             >  2:845276/2‑150 (MQ=255)
                                                                                                                                           tACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtg            >  2:1566796/1‑150 (MQ=255)
                                                                                                                                           tACTGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtg            >  2:344211/1‑150 (MQ=255)
                                                                                                                                             cTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgcc          <  2:3904248/150‑1 (MQ=255)
                                                                                                                                              tGGCTGTCATTACCGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccg         >  1:860207/1‑150 (MQ=255)
                                                                                                                                              tGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCGTTGGCAATGGCGGAACCGACGATATAGCCCTAGGCGGCGGTGGGGGCTTCTTCATGCCAGTCGTAGGCGTAGcgtgccg         >  1:2059114/1‑150 (MQ=255)
                                                                                                                                                 cGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccgtgc      <  1:3747424/150‑1 (MQ=255)
                                                                                                                                                  ggTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc     >  2:4193807/1‑150 (MQ=255)
                                                                                                                                                  ggTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCgtgccgtgcc     <  2:2332869/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCACCAGCGAGGCCGCCGCTCGCGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGT  >  NC_000913/300543‑300841

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