Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,903 A→G A45A (GCA→GCG yjjJ → toxin YjjJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,9030AG99.0% 356.4 / ‑5.4 99A45A (GCA→GCGyjjJtoxin YjjJ
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (49/49);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAA  >  NC_000913/4621767‑4622039
                                                                                                                                        |                                                                                                                                        
tAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCgg                                                                                                                                       >  1:716602/1‑140 (MQ=255)
  aTGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGc                                                                                                                                     <  1:141956/140‑1 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTCCTGCAAGGGCCGCGTTCTGCCCGGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    <  2:1113383/140‑1 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    >  2:818503/1‑140 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    >  1:770749/1‑140 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    <  1:818503/140‑1 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    >  1:939283/1‑140 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    >  1:992652/1‑140 (MQ=255)
   tGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCa                                                                                                                                    <  2:716602/140‑1 (MQ=255)
    gagcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCTAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCaa                                                                                                                                   <  2:941691/140‑1 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  2:996634/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  1:1279976/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  1:139592/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  1:480500/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  1:157065/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  2:645822/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  2:1067002/1‑140 (MQ=255)
     agcgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  1:41989/1‑140 (MQ=255)
       cgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgc                                                                                                                                >  2:327060/1‑140 (MQ=255)
       cgagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgc                                                                                                                                >  2:1248040/1‑140 (MQ=255)
        gggcTGACTGATCTTTAACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               >  2:699839/3‑140 (MQ=255)
        gagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               <  2:444525/140‑1 (MQ=255)
        gagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               <  1:801574/140‑1 (MQ=255)
        gagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               <  1:525638/140‑1 (MQ=255)
        gagcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               <  2:829684/140‑1 (MQ=255)
         agcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               >  1:1195980/1‑139 (MQ=255)
         agcTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAcgcg                                                                                                                               <  2:1195980/139‑1 (MQ=255)
            tgactgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGatat                                                                                                                           >  2:554213/1‑140 (MQ=255)
             gcctgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGAGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATg                                                                                                                          <  2:948092/138‑1 (MQ=255)
               ctgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        <  1:491806/140‑1 (MQ=255)
               ctgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        <  2:992652/140‑1 (MQ=255)
               ctgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        <  1:1266572/140‑1 (MQ=255)
               ctgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        <  1:70842/140‑1 (MQ=255)
                tgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  <  1:1284126/129‑1 (MQ=255)
                tgaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAAc                                                                                                                                  >  2:1284126/1‑129 (MQ=255)
                 gaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        <  1:98705/138‑1 (MQ=255)
                 gaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCa                                                                                                                        >  2:98705/1‑138 (MQ=255)
                 gaTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCAct                                                                                                                      >  2:1177995/1‑140 (MQ=255)
                  aTCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCActg                                                                                                                     <  1:71541/140‑1 (MQ=255)
                    cTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCActgct                                                                                                                   <  2:157065/140‑1 (MQ=255)
                      tttACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACtgctgc                                                                                                                 >  1:1260951/1‑140 (MQ=255)
                        tACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGt                                                                                                               >  2:997581/1‑140 (MQ=255)
                        tACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGt                                                                                                               >  1:902863/1‑140 (MQ=255)
                            gCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCtt                                                                                                           >  2:332147/1‑140 (MQ=255)
                            gCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCtt                                                                                                           >  2:607872/1‑140 (MQ=255)
                             cAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTa                                                                                                          >  1:568286/1‑140 (MQ=255)
                              aaGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTAt                                                                                                         >  1:978330/1‑140 (MQ=255)
                              aaGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTAt                                                                                                         >  1:436629/1‑140 (MQ=255)
                               aGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATc                                                                                                        >  1:1166680/1‑140 (MQ=255)
                               aGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATc                                                                                                        >  1:600120/1‑140 (MQ=255)
                                 ggCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATcgc                                                                                                      >  2:1169784/1‑140 (MQ=255)
                                   ccGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCgg                                                                                                    <  2:1166680/140‑1 (MQ=255)
                                    cgcgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGa                                                                                                   >  2:443829/1‑140 (MQ=255)
                                     gcgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGaa                                                                                                  >  1:646816/1‑140 (MQ=255)
                                     gcgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGaa                                                                                                  >  2:383618/1‑140 (MQ=255)
                                      cgTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAt                                                                                                 <  1:554213/140‑1 (MQ=255)
                                       gTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAtt                                                                                                <  1:227265/140‑1 (MQ=255)
                                       gTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAAtt                                                                                                <  1:657916/140‑1 (MQ=255)
                                          cTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTGGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATCGAg                                                                                             >  2:167384/1‑140 (MQ=255)
                                           tGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                                                                                            <  1:327060/140‑1 (MQ=255)
                                           tGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                                                                                            <  2:480500/140‑1 (MQ=255)
                                           tGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGc                                                                                            <  2:139592/140‑1 (MQ=255)
                                                cGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTAtt                                                                                       >  1:425932/1‑140 (MQ=255)
                                                         gcCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTgg                                                                              <  1:1112938/140‑1 (MQ=255)
                                                          ccAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGc                                                                             <  1:1354126/140‑1 (MQ=255)
                                                           cAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCg                                                                            >  1:16918/1‑140 (MQ=255)
                                                            aGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCgg                                                                           <  2:600120/140‑1 (MQ=255)
                                                            aGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCgg                                                                           <  1:997581/140‑1 (MQ=255)
                                                            aGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCgg                                                                           <  1:1177995/140‑1 (MQ=255)
                                                             gCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGc                                                                             <  2:317938/137‑1 (MQ=255)
                                                             gCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGc                                                                             >  1:317938/1‑137 (MQ=255)
                                                                 cTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGa                                                                      <  2:43151/140‑1 (MQ=255)
                                                                 cTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGa                                                                      <  1:1169784/140‑1 (MQ=255)
                                                                 cTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGa                                                                      <  2:978330/140‑1 (MQ=255)
                                                                  tGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGAc                                                                     <  2:41989/140‑1 (MQ=255)
                                                                    gCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGa                                                                   <  2:1260951/140‑1 (MQ=255)
                                                                    gCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGa                                                                   <  1:1251827/140‑1 (MQ=255)
                                                                    gCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGa                                                                   <  2:646816/140‑1 (MQ=255)
                                                                        tcagtcagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGTGTATTCCCGTCTGGCGGGTGGACGATAcc                                                               >  2:625243/1‑140 (MQ=255)
                                                                            tcagGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGaa                                                           <  1:332147/140‑1 (MQ=255)
                                                                              agGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAg                                                         <  1:383618/140‑1 (MQ=255)
                                                                               ggCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAgg                                                        <  2:425932/140‑1 (MQ=255)
                                                                                    cGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCAt                                                   >  1:1101557/1‑140 (MQ=255)
                                                                                    cGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCAt                                                   >  1:602828/1‑140 (MQ=255)
                                                                                      ttCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATaa                                                 >  2:327198/1‑140 (MQ=255)
                                                                                      ttCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATaa                                                 >  2:297040/1‑140 (MQ=255)
                                                                                           aCGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCg                                            <  1:443829/140‑1 (MQ=255)
                                                                                                 tgttgCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACa                                      <  1:607872/140‑1 (MQ=255)
                                                                                                       cAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGt                                >  1:1039912/1‑140 (MQ=255)
                                                                                                           gaagaTCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTtgtg                            >  1:1198153/1‑140 (MQ=255)
                                                                                                               aTCGGGTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTgg                        <  2:1198153/140‑1 (MQ=255)
                                                                                                                   ggTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGc                    >  2:699063/1‑140 (MQ=255)
                                                                                                                    gTGATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCa                   >  2:1005494/1‑140 (MQ=255)
                                                                                                                      gATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAgg                 >  1:1062856/1‑140 (MQ=255)
                                                                                                                      gATTCGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAgg                 >  2:638881/1‑140 (MQ=255)
                                                                                                                          cGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAg             <  2:902863/140‑1 (MQ=255)
                                                                                                                          cGCTTTGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAg             <  2:568286/140‑1 (MQ=255)
                                                                                                                            cTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGtt           <  2:626234/140‑1 (MQ=255)
                                                                                                                              ttGGTAAAGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGt         <  1:1170178/140‑1 (MQ=255)
                                                                                                                                     aGCGCGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTaa  <  2:1101557/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                        
TAATGAGCGAGCTGACTGATCTTTTACTGCAAGGGCCGCGTTCTGCCCCGGAATTGCGCCAGCGTCTGGCAATCAGTCAGGCGACGTTCTCACGCCTTGTTGCCAGAGAAGATCGGGTGATTCGCTTTGGTAAAGCACGGGCAACGCGATATGCACTGCTGCGTCCTTATCGCGGAATTGAGCGTATTCCCGTCTGGCGGGTGGACGATACCGGAAAGGCGCATAAATTCGCCGACATCCGGTTGTGCTGGCCGCAGGGAAGTTGTCTGGTAA  >  NC_000913/4621767‑4622039

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