Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,159,421 A→T 15.9% I346F (ATT→TTT)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,159,4210AT15.9% 164.1 / 2.2 64I346F (ATT→TTT) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base A (36/16);  new base T (3/8);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGATGTGCGCAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGTTGCAAGCCG  >  NC_000913/2159183‑2159655
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
gCAGGATGTGCGCAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGa                                                                                                                                                                                                                                            <  2:1226014/239‑1 (MQ=255)
   ggATGTGCGCAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGAt                                                                                                                                                                                                                                           >  2:11530/1‑237 (MQ=255)
    gATGTGCGCAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGAt                                                                                                                                                                                                                                           >  1:598734/1‑236 (MQ=255)
         gcgcAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGTTTCtggggg                                                                                                                                                                                                                                   >  2:559632/1‑239 (MQ=255)
         gcgcAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtggtgg                                                                                                                                                                                                                                   >  1:1232607/1‑239 (MQ=255)
         gcgcAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCtggtgg                                                                                                                                                                                                                                   >  1:729766/1‑239 (MQ=255)
           gcAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTggtggt                                                                                                                                                                                                                                  <  1:686487/238‑1 (MQ=255)
               cGCCGGGATGACCAACACCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGTGCAAAATTACCGGTGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGTTTCTGGTGGTGtttt                                                                                                                                                                                                                             >  1:690033/1‑239 (MQ=255)
                    ggATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTAtt                                                                                                                                                                                                                         >  1:72070/1‑238 (MQ=255)
                                  aaaCCGGCTATTTTACTGATGATCCGCAAACGGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTcg                                                                                                                                                                                                          <  2:108944/239‑1 (MQ=255)
                                    tccGGGTTTTGTACTGACGACTCTCAATCTGGCGGAAGCCAATATTACCCAGTCGGTTGGCCTCATCCGGGCAAAACTACCGGAGTCGCAGGACACCATTCCGGCGGCGATTGATTTGCAGATTGCCCAGGATCGCCCCTCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCCAATTCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTcgcg                                                                                                                                                                                                        <  2:861508/238‑1 (MQ=255)
                                     ccGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCgcgc                                                                                                                                                                                                       <  1:812276/239‑1 (MQ=255)
                                         cTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACt                                                                                                                                                                                                   >  1:252856/1‑239 (MQ=255)
                                                 cttatgatCCGCAAACTCCCGGACGCCAATTTTCCCCAGACGGTTGACAGCACCCGGGAAAAATTACGGGAGTTGCAGGAACCCATTCCGGCGGCGTTTGATCTGAAAATTGCCCAGGATCGCCCCCCCACCATTCGCGCCTGGGTGGAAGAAGTCGAGCAACCGCTGATTATCTCGGGGGCGCTGGTGATTCTGGGGGTGTTTTTATCCCTGCGCTCGGGTCGCGCCACTATTATTcc                                                                                                                                                                                           <  2:163595/236‑1 (MQ=255)
                                                    atgatCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCccg                                                                                                                                                                                         >  1:1154038/1‑238 (MQ=255)
                                                      gatCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCccgccg                                                                                                                                                                                      >  1:616019/1‑239 (MQ=255)
                                                      gatCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCccgccg                                                                                                                                                                                      >  1:163595/1‑239 (MQ=255)
                                                      gatCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCccgccg                                                                                                                                                                                      >  2:1024120/1‑239 (MQ=255)
                                                      gatCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCccgccg                                                                                                                                                                                      >  1:1048223/1‑239 (MQ=255)
                                                                   ccGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGt                                                                                                                                                                         >  2:417212/1‑239 (MQ=255)
                                                                                ttatCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCgcgc                                                                                                                                                                                                       <  1:745649/196‑1 (MQ=255)
                                                                                   tCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACg                                                                                                                                                         >  1:991308/1‑239 (MQ=255)
                                                                                        ggtgTTGACAGCACCGGGGAAAAATTACAGGGGTTGCAGAAACCCACTCCGGCGGCGATGGATCTCCAATTTGCCCGGGATCGCTCCCCCCCCATCCGCGCCTCGCGGAAAGAATGGGAGAAAACGTTTATTCTCTCGGTGGCGTTGGTGTTTCGGGTGGTGTTTTTATCCCTGCGCGGGGGCCGCGCCACTTTTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTg                                                                                                                                                     <  2:233704/235‑1 (MQ=255)
                                                                                        aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGTTTCTGGTGGCGTTTGTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGTTGCCGGTTTAGCTGGTTGGTACggt                                                                                                                                                       >  2:1196759/1‑234 (MQ=255)
                                                                                        aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgc                                                                                                                                                    >  1:578833/1‑239 (MQ=255)
                                                                                        aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgc                                                                                                                                                    >  1:518049/1‑239 (MQ=255)
                                                                                        aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgc                                                                                                                                                    >  1:411439/1‑239 (MQ=255)
                                                                                        aCGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTgc                                                                                                                                                    >  1:1127990/1‑239 (MQ=255)
                                                                                          ggTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGAGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTTCCGGTTTCGCTGCTTGGTACGTTTgcg                                                                                                                                                   >  2:1016494/1‑238 (MQ=255)
                                                                                              gACAGCAGCGGGGCAAAATCAGCGGAGGGGCGGAAACCCACGCCGGGGGCGATGGATCGCAAAATGGCCCAGGCTCGCCCCCACACCAGTCGCGCCGCGGGGGAAGAGGGGGAGCAAACGTTTATTACCGCGGGGGCGCGGGTGTTGTGGGTGTTGTTTTTATTCCTGCGCGCGGGGCGCGCCACTATTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGa                                                                                                                                               <  2:268178/238‑1 (MQ=255)
                                                                                                    ccgcGGGCAATAGTACGGGGGTGGAAGGAAACTCTGCGGGCGGCGATGGATCTCATAATGCCCAAGGATCGCTCCCCCCCATCTCGCCGTTCGCGGGAAGAAGGGGAGAAACGTCGGAGATTCTCGGGGGCGCGGGTGTTTTTGGGGGGTTTTTTATTCCCCCTCTGGGGTCGGGCACTTATTAGTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTAcc                                                                                                                                         <  2:292554/235‑1 (MQ=255)
                                                                                                      ccGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCTCCACTATTTTTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTAGGTTTGCGGCGATGTCCCttt                                                                                                                                      >  1:975776/1‑237 (MQ=255)
                                                                                                      ccGGGCAAAATTACCGGAGTTGCAGGAAACCATCCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCtgt                                                                                                                                      >  1:110265/1‑239 (MQ=255)
                                                                                                             aaaTTCCCGGATTTGCAGGAAACCATTCCGGCGGCGATGGATCTGCAAATTGCCCAGGATCCCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGAAAACGCGAATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGAtt                                                                                                                               <  1:11530/239‑1 (MQ=255)
                                                                                                                ttccggggggtGCAGGAACCCGTTCAGGCGGCGATTGATCTGCAAATCGCCCAGGACCCCCCCCCCCCCATTCGCGCCTCGCTGAAGGAGGTCGAGCAAACGCTGATTATCTGGGGGGCGGTGGTGATTCTGGTGGTTTTTTTTGTCCGGCGCTGGGGGCGCGCCACTATTTTTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAg                                                                                                                            <  2:518049/229‑1 (MQ=255)
                                                                                                                   ccGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGtc                                                                                                                          >  1:419606/1‑238 (MQ=255)
                                                                                                                           gCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGCTTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGCTTGGTACGTTTGCGGCGCTGTCCCTGTGCTGATTCAGTCTCAAtca                                                                                                                   >  2:76883/1‑235 (MQ=255)
                                                                                                                                 aaCCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGt                                                                                                           >  2:722837/1‑239 (MQ=255)
                                                                                                                                               caaTTGATCTGAAAATGGCCCAGGTTCGCTCCCCCCCAGTTCGCGCCGGGGTGGAAGAAGTGGAGAAAACTTGGTTTATCGCGGTGGCGCGGGTGTTTCTGGTGGTTTTTTTTTCCCGGGCCTGGGGCCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCAc                                                                                              <  2:578833/236‑1 (MQ=255)
                                                                                                                                                   tGATCTACAAATTGCCCAGGACCGCTCCCCCACCTGTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTTTCTCGGTGGCGCTGTTTATTCTGGTGTTGTTTTTTTTCCTGCGCTCGGGTCGCGCCACTATTATGCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCg                                                                                         <  2:156870/239‑1 (MQ=255)
                                                                                                                                                        ggCAAATGTGCCCGGGCTCCCCCCCCCACCATTCGGCCCTCGCTGGAAGAGGTCGAGAAAACGCTGTTTGTCTGGGTGGCGGTGGTGATTCTGGTGGTGTTTTTATCCCTGCTCTGGGGTCGCCCCATTTTTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTAc                                                                                     <  2:1151806/238‑1 (MQ=255)
                                                                                                                                                           aaaTTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGt                                                                                 >  1:336762/1‑239 (MQ=255)
                                                                                                                                                            aaTTCCCCAGGACCCCCCCCCCCCCTGCCGCCCCTCGTTGAAAGAGGTCGACCAACCGCTGTTTTTCTCGGGGGCGCTGTTTTTTCTGGTGTTTTTTTTATTCCTGCCCTCGGTCCGCGCCATTTTTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGtt                                                                                <  2:336762/239‑1 (MQ=255)
                                                                                                                                                                 gCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCgt                                                                            >  2:65833/1‑239 (MQ=255)
                                                                                                                                                                             cTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTAc                                                                                                                                          >  1:1005031/1‑166 (MQ=255)
                                                                                                                                                                             cTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATg                                                                    >  1:33881/1‑236 (MQ=255)
                                                                                                                                                                             cTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGat                                                                  >  2:1145629/1‑237 (MQ=255)
                                                                                                                                                                                        ttCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCgtggt                                                       >  1:128082/1‑238 (MQ=255)
                                                                                                                                                                                          ccgccctgcggtGGAAGAAGTGGAGAAAACGCTGTTTGTCTGGGGGCGCGTGGTGTTGTTGGTGGTTTTTTTTTTCCGGGCCTGGGGTCGCGCCACTTTTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCgtggtg                                                      <  2:252856/227‑1 (MQ=255)
                                                                                                                                                                                             cgcCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTgg                                                  >  1:655790/1‑239 (MQ=255)
                                                                                                                                                                                               ccTCGCTGGAAGAAGTCGAGCAAACGTTGATTATCTCGGTGGCGCTGGTGTTTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGaa                                                <  1:417212/239‑1 (MQ=255)
                                                                                                                                                                                                   gCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCGGGTGTTGTTTTTATTCCTGCGCTCGGGTCCCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACa                                            <  1:559632/239‑1 (MQ=255)
                                                                                                                                                                                                            agTCGAGTAAACGCTGATTATCTCGGTGGCGCTGGTTATTCGGGTGGTTTTTTTATTCCTGCGCTCGGGTCGCGCCCTTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACg                                     <  1:1016494/237‑1 (MQ=255)
                                                                                                                                                                                                            agTCGAGAAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTGGGGTCGCGCCACTATTATCCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACg                                     <  1:777917/237‑1 (MQ=255)
                                                                                                                                                                                                              tCGAGCAAACGCTGATTATCTCGGTGGCGCTGGCGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGCCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGCCGCCATCGTGGTGCTGGCACCCATTGCACGTCa                                  >  1:43980/1‑238 (MQ=255)
                                                                                                                                                                                                               cGAGCAAGCGCTGATTTGCTCGGTGTCGCTGGTGATTTGGGTGTTGTTTTTATTCCGCCGCTCGGGTCGCGCCCCTATTATGCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCAt                                 <  1:1055701/238‑1 (MQ=255)
                                                                                                                                                                                                                    aaaCGTTGATTATCGCGGTGGCGTTGGGGTTTGTGGGGGTGTTTTTCTTCCTGCGCGCGGGTCGCGCCATTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGa                            <  2:625105/238‑1 (MQ=255)
                                                                                                                                                                                                                              ttATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAAt                   <  1:722837/238‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACg                                                                 >  1:919008/1‑185 (MQ=255)
                                                                                                                                                                                                                                       tGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGtt          >  1:720785/1‑238 (MQ=255)
                                                                                                                                                                                                                                       tGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGtt          >  1:654206/1‑238 (MQ=255)
                                                                                                                                                                                                                                       tGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGtt          >  1:611322/1‑238 (MQ=255)
                                                                                                                                                                                                                                            cTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCAGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGTTGCAAg     >  1:253462/1‑238 (MQ=255)
                                                                                                                                                                                                                                                gTGATTCTGGTGGTGTTTTTTTTCCTGCGCGCGGGTGGCGCCATTTTTATGCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGTTGCAAGCCg  <  1:696026/238‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTa                                                                                                                                           >  1:236614/1‑101 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
GCAGGATGTGCGCAACGCCGGGATGACCAACGCCAAACCGGCTATTTTACTGATGATCCGCAAACTGCCGGAAGCCAATATTATCCAGACGGTTGACAGCATCCGGGCAAAATTACCGGAGTTGCAGGAAACCATTCCGGCGGCGATTGATCTGCAAATTGCCCAGGATCGCTCCCCCACCATTCGCGCCTCGCTGGAAGAAGTCGAGCAAACGCTGATTATCTCGGTGGCGCTGGTGATTCTGGTGGTGTTTTTATTCCTGCGCTCGGGTCGCGCCACTATTATTCCCGCCGTTTCGGTGCCGGTTTCGCTGATTGGTACGTTTGCGGCGATGTACCTGTGCGGATTCAGTCTCAATAACCTTTCGTTAATGGCGCTCACCATCGCTACTGGTTTCGTGGTGGATGACGCCATCGTGGTGCTGGAAAACATTGCACGTCATCTGGAAGCGGGAATGAAACCGTTGCAAGCCG  >  NC_000913/2159183‑2159655

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