Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,282,331 A→G 10.7% L26L (CTT→CTC PP_2883 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,282,3310AG10.7% 228.4 / 14.4 84L26L (CTT→CTCPP_2883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base G (4/5);  total (45/39)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGT  >  NC_002947/3282184‑3282478
                                                                                                                                                   |                                                                                                                                                    
cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                    <  1:1980801/150‑1 (MQ=255)
cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                    <  1:1953186/150‑1 (MQ=255)
   ccAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                                 <  1:1322252/150‑1 (MQ=255)
       tgtTCGAGGCCGGCCTGTTCCCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGc                                                                                                                                            >  2:37007/1‑150 (MQ=255)
        gtTCGAGGCCGGCCTGTTCGCCGTCGGCGGCTTGGCAACGAACGTAATGCGCCAGGCCTTCGGTATCCAGGGCCAACTTCCGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCACCACCTTGAggtggctgca                                                                                                                                           <  2:1188138/150‑6 (MQ=255)
            gAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGcc                                                                                                                                       >  2:1146081/1‑150 (MQ=255)
             aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCa                                                                                                                                      >  1:220319/1‑150 (MQ=255)
             aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCa                                                                                                                                      >  2:817611/1‑150 (MQ=255)
              ggccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCaa                                                                                                                                     <  1:981060/150‑1 (MQ=255)
              ggccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCAGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCaa                                                                                                                                     >  1:1250379/1‑150 (MQ=255)
                ccggccTGTTCGCCTTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCAAGGGCTAACGGCAGTTGCCAACCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAgg                                                                                                                                   <  1:781818/150‑1 (MQ=255)
                ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGTGGCGTGCACCACCTTGAGGTGGCTGCACGccgcga                                                                                                                                   <  2:140069/150‑5 (MQ=38)
                        ttCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGCTGCAc                                                                                                                                          >  1:2259733/1‑136 (MQ=255)
                        ttCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGCTGCAc                                                                                                                                          <  2:2259733/136‑1 (MQ=255)
                          cGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTgc                                                                                                                         >  2:92344/1‑150 (MQ=255)
                           gCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTgcg                                                                                                                        <  1:37007/150‑1 (MQ=255)
                            ccGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACGTCGATCAGGCCGCCATCGCGGCGTGCAGCCAAATCAAACTGCTATACTCCATGGGCCTCAtttttt                                                                                                                       >  1:1133504/1‑146 (MQ=255)
                              gTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAg                                                                                                                     <  1:1633686/150‑1 (MQ=255)
                              gTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAg                                                                                                                     >  1:1241205/1‑150 (MQ=255)
                                  ggggCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCAc                                                                                                                 <  2:843852/150‑1 (MQ=255)
                                    ggCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAgg                                                                                                                                                   <  2:1505138/115‑1 (MQ=255)
                                    ggCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAgg                                                                                                                                                   >  1:1505138/1‑115 (MQ=255)
                                         ggCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGa                                                                                                          <  1:934622/150‑1 (MQ=255)
                                            aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGcc                                                                                                                                       <  2:2226718/118‑1 (MQ=255)
                                            aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGcc                                                                                                                                       >  1:2226718/1‑118 (MQ=255)
                                            aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                       >  1:194242/1‑150 (MQ=255)
                                                aacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCg                                                                                                   >  1:951920/1‑150 (MQ=255)
                                                 acgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                  <  1:1558463/150‑1 (MQ=255)
                                                   gCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGc                                                                                                <  1:1309904/150‑1 (MQ=255)
                                                     aaTGCGCCAGGCCTTCGGTATCCAGGGCCACCGAACGTTGACCCGATTGGTGGACGACATCGATCAGGCCGCCGTCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGt                                                                                              >  1:2224542/1‑150 (MQ=255)
                                                     aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                    >  1:1977163/1‑97 (MQ=255)
                                                     aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                    <  2:1977163/97‑1 (MQ=255)
                                                       tGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCATGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTc                                                                                            >  1:413100/1‑150 (MQ=255)
                                                              ggCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGCCATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                                     <  1:92344/150‑1 (MQ=255)
                                                               gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGACGTGCAGCCACCTTAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCCTTCGATGAACg                                                                                    >  2:1116247/1‑150 (MQ=255)
                                                                  ttCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCt                                                                                 >  2:1858995/1‑150 (MQ=255)
                                                                  ttCGGTATCCAGGGCCAACGTCAGATGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCt                                                                                 >  1:666942/1‑150 (MQ=255)
                                                                      gTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                             <  2:1241205/150‑1 (MQ=255)
                                                                       tATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTggg                                                                            >  1:355276/1‑150 (MQ=255)
                                                                          ccAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCa                                                                         >  2:1211203/1‑150 (MQ=255)
                                                                            aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                                       >  1:754994/1‑150 (MQ=255)
                                                                                ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                   >  2:1329909/1‑150 (MQ=255)
                                                                                ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                   <  2:951920/150‑1 (MQ=255)
                                                                                 cAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAg                                                                  <  1:1858995/150‑1 (MQ=255)
                                                                                     gTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAAt                                                              >  2:422629/1‑150 (MQ=255)
                                                                                             ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtg                                                                                >  2:1788026/1‑124 (MQ=255)
                                                                                             ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtg                                                                                <  1:1788026/124‑1 (MQ=255)
                                                                                             ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGc                                                      >  2:1095786/1‑150 (MQ=255)
                                                                                               aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                                    <  1:189055/150‑1 (MQ=255)
                                                                                                gCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgtg                                                   >  1:387679/1‑150 (MQ=255)
                                                                                                 cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgtgt                                                  >  1:1726871/1‑150 (MQ=255)
                                                                                                       aTGCCGACATCGATCAGGCCGCCATCGAGGGGTGCTGCCACCGCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGgcgc                                            <  2:194242/150‑1 (MQ=255)
                                                                                                          ccGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                         <  2:387679/150‑1 (MQ=255)
                                                                                                          ccGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                         <  1:817611/150‑1 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                  <  1:837526/92‑1 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                  >  2:837526/1‑92 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                        >  1:140129/1‑150 (MQ=255)
                                                                                                            gACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgatg                                       >  1:144551/1‑150 (MQ=255)
                                                                                                                   aTCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCgg                                >  1:1411197/1‑150 (MQ=255)
                                                                                                                        gccgccATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                    >  2:160826/1‑61 (MQ=255)
                                                                                                                        gccgccATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAgg                                                                                                                    <  1:160826/61‑1 (MQ=255)
                                                                                                                         ccgccATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                             >  1:126159/1‑99 (MQ=255)
                                                                                                                         ccgccATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                             <  2:126159/99‑1 (MQ=255)
                                                                                                                             cATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCgg                      <  2:1250379/150‑1 (MQ=255)
                                                                                                                               tCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGc                    <  2:355276/150‑1 (MQ=255)
                                                                                                                               ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGa                                                                 >  1:1981264/5‑105 (MQ=25)
                                                                                                                               ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGa                                                                 >  1:1981422/5‑105 (MQ=25)
                                                                                                                               ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGa                                                                 <  2:1981422/101‑1 (MQ=25)
                                                                                                                               ccttggcgTGCACCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGa                                                                 <  2:1981264/101‑1 (MQ=25)
                                                                                                                                 ttggcgTGCACCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg                  >  1:2329023/3‑150 (MQ=255)
                                                                                                                                     cgTGCACCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                  >  2:1759807/1‑66 (MQ=18)
                                                                                                                                     cgTGCACCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCgg                                                                                                  <  1:1759807/66‑1 (MQ=18)
                                                                                                                                         cAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt          <  2:1411197/150‑1 (MQ=255)
                                                                                                                                          aGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTg         >  1:913970/1‑150 (MQ=255)
                                                                                                                                           gCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGa        >  1:2094809/1‑150 (MQ=255)
                                                                                                                                              aCCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCTGGCCGATGCTCTTGgaaa     >  2:870280/1‑149 (MQ=255)
                                                                                                                                              aCCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt                >  1:583889/1‑139 (MQ=255)
                                                                                                                                              aCCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt                <  2:583889/139‑1 (MQ=255)
                                                                                                                                                cTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt          <  2:1710984/143‑1 (MQ=255)
                                                                                                                                                cTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt          >  1:1710984/1‑143 (MQ=255)
                                                                                                                                                cTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAAACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtg   >  1:1758088/1‑150 (MQ=255)
                                                                                                                                                 tCAAGTTGCTGCACGCAAAGGGCCACATTGCGCAGTCACTGGTAGACTTCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGCCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtgt  <  1:1116247/150‑1 (MQ=255)
                                                                                                                                                 tCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtgt  >  1:589847/1‑150 (MQ=255)
                                                                                                                                                 tCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtgt  <  1:1095786/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGGTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTGT  >  NC_002947/3282184‑3282478

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