Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,432,050 T→G 5.8% intergenic (+37/‑147) PP_3924 → / → PP_3925 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,432,0500TG 5.8% 583.9 / 19.9 174intergenic (+37/‑147)PP_3924/PP_3925hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (70/93);  new base G (6/4);  total (76/98)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

Only 100 of 175 total aligned reads displayed.

GGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATG  >  NC_002947/4431903‑4432199
                                                                                                                                                   |                                                                                                                                                      
        aaGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCg                                                                                                                                              <  1:2711402/150‑1 (MQ=255)
         aGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCgg                                                                                                                                             <  2:4340343/150‑1 (MQ=255)
         aGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCgg                                                                                                                                             <  1:4115644/150‑1 (MQ=255)
           cggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGtt                                                                                                                                          <  2:1604896/150‑1 (MQ=255)
            ggcgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttt                                                                                                                                         <  1:3401590/150‑1 (MQ=255)
              cgAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttt                                                                                                                                       <  1:1290094/150‑1 (MQ=255)
                aaaCCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGttttttt                                                                                                                                     >  2:2192803/1‑150 (MQ=255)
                   ccGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGc                                                                                                                                  <  1:3730177/150‑1 (MQ=255)
                    cGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGcc                                                                                                                                 >  1:648847/1‑150 (MQ=255)
                      aCGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCt                                                                                                                               <  1:1500921/150‑1 (MQ=255)
                        ggCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGc                                                                                                                             <  2:1271249/150‑1 (MQ=255)
                           aTGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCAtt                                                                                                                          >  2:607894/1‑150 (MQ=255)
                             gCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCg                                                                                                                        <  2:1354783/150‑1 (MQ=255)
                                 gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGc                                                                                                                    >  1:3079236/1‑150 (MQ=255)
                                 gCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGc                                                                                                                    <  2:1225347/150‑1 (MQ=255)
                                  ccTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCa                                                                                                                   <  1:3242170/150‑1 (MQ=255)
                                  ccTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCa                                                                                                                   <  2:1608966/150‑1 (MQ=255)
                                  ccTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCa                                                                                                                   <  1:2070185/150‑1 (MQ=255)
                                    ttCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGt                                                                                                                 <  1:1518356/150‑1 (MQ=255)
                                       aTGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGAt                                                                                                              >  2:2281724/1‑150 (MQ=255)
                                         gggTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATcc                                                                                                            >  1:2390220/1‑150 (MQ=255)
                                            tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATccacc                                                                                                         <  1:3348723/150‑1 (MQ=255)
                                            tGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTAGGTCAAGTGATccaac                                                                                                         >  1:959437/1‑148 (MQ=255)
                                              gCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAg                                                                                                       <  1:3055679/150‑1 (MQ=255)
                                               cAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGc                                                                                                      <  1:226943/150‑1 (MQ=255)
                                                          ttGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCAc                                                                                           >  2:4027741/1‑150 (MQ=255)
                                                             aaGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTa                                                                                        <  1:107804/150‑1 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTTCCCGGCATTGGGTGCAGTGATCCACCAGCAACACCCACAACGAttt                                                                                     >  1:1613077/1‑150 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                     <  2:4509002/150‑1 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                     <  1:3582912/150‑1 (MQ=255)
                                                                gCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTAttt                                                                                     <  2:4329292/150‑1 (MQ=255)
                                                                 cGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTg                                                                                    >  1:632730/1‑150 (MQ=255)
                                                                       gatgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACAt                                                                              >  2:314894/1‑150 (MQ=255)
                                                                         tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATgg                                                                            >  1:2758425/1‑150 (MQ=255)
                                                                         tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATgg                                                                            >  2:2298811/1‑150 (MQ=255)
                                                                         tgCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATgg                                                                            <  1:1654926/150‑1 (MQ=255)
                                                                           cTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGa                                                                          <  2:3238165/150‑1 (MQ=255)
                                                                           cTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGa                                                                          <  1:444627/150‑1 (MQ=255)
                                                                                 aTGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCACACCGATTTGGCACATGGGACGATTg                                                                    >  1:120494/1‑150 (MQ=255)
                                                                                  tGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGt                                                                   >  2:4028534/1‑150 (MQ=255)
                                                                                   gAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGtt                                                                  <  1:1174634/150‑1 (MQ=255)
                                                                                    aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 >  2:139398/1‑150 (MQ=255)
                                                                                    aaaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTc                                                                 >  2:1160689/1‑150 (MQ=255)
                                                                                     aaCTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTcc                                                                <  2:850192/150‑1 (MQ=255)
                                                                                        tCAGCGGGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTATTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAAc                                                             <  1:4240588/150‑1 (MQ=255)
                                                                                        tCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAAc                                                             >  1:2829189/1‑150 (MQ=255)
                                                                                           tcGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGAACGAAAAACCGAAGCCAGACGGCTTCGGTAATTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCt                                                          <  2:4131047/149‑1 (MQ=255)
                                                                                           gCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCt                                                          <  1:1145853/150‑1 (MQ=255)
                                                                                                 tGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCGTCAGGCTTCGGtttttc                                                                                                                                      <  1:594088/68‑2 (MQ=21)
                                                                                                  gTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAc                                                   <  2:3719790/150‑1 (MQ=255)
                                                                                                   ttCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAcc                                                  <  1:1701802/150‑1 (MQ=255)
                                                                                                   ttCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGAcc                                                  >  2:4150738/1‑150 (MQ=255)
                                                                                                     cccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAt                                                >  1:3895869/1‑150 (MQ=255)
                                                                                                     cccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAt                                                <  2:4477720/150‑1 (MQ=255)
                                                                                                      ccGTCCTGAGGGCTGGGAATACCGGTTGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                               >  2:1462193/1‑150 (MQ=255)
                                                                                                      ccGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAtt                                               <  1:1723242/150‑1 (MQ=255)
                                                                                                       cGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTa                                              <  2:1038269/150‑1 (MQ=255)
                                                                                                       cGTCCTGAGGCCTGGGACTACCGGACGCTCGAACCACCGCAGCCTGCCGGCTTCGGTTTTTTTGGCCCGGCATTCGGTGCAGTGATCCACCAGCAACACCCTCCCCTAGTTGGCACAGGGGACGATTGTGCCAACCCTGGCGGACCATTa                                              <  2:4030764/150‑1 (MQ=255)
                                                                                                         tCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                            <  2:3947622/150‑1 (MQ=255)
                                                                                                         tCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                            <  2:2028600/150‑1 (MQ=255)
                                                                                                         tCCAGAGGGCTGGGAATACCGGATGATCGAAAAACCCAAGCCTGACGGCTACGGATATTTTTGCCCTGCATTCGGTGCAGTGATCCACCACCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGa                                            <  2:3404979/150‑1 (MQ=255)
                                                                                                              aGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTAcc                                       <  1:3888958/150‑1 (MQ=255)
                                                                                                               gggCGGGGATTCCGGGATGATCGAAATAGCGAAGCGGTAGGGCGTGGGTTTTTTTTGCCCTTCATTCGGTGCCGTGTTCGACCAGCTACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGTTTaggg                                      <  2:2602272/150‑4 (MQ=255)
                                                                                                                ggCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGCTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCa                                     >  2:697464/1‑150 (MQ=255)
                                                                                                                 gCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCGTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCaa                                    >  1:586732/1‑150 (MQ=255)
                                                                                                                  cTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAAt                                   <  1:2806717/150‑1 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  >  1:3766598/1‑150 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:363237/150‑1 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  >  2:173041/1‑150 (MQ=255)
                                                                                                                   tGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:1921078/150‑1 (MQ=255)
                                                                                                                   tGGCAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATa                                  <  1:428662/150‑1 (MQ=255)
                                                                                                                     ggAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTc                                                                                             >  1:1131402/1‑89 (MQ=255)
                                                                                                                     ggAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTc                                                                                             <  2:1131402/89‑1 (MQ=255)
                                                                                                                          aCCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTg                           <  1:746964/150‑1 (MQ=255)
                                                                                                                            cGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCt                         >  1:1254426/1‑150 (MQ=255)
                                                                                                                             ggATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTa                        >  1:2528898/1‑150 (MQ=255)
                                                                                                                               atgatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATc                      <  1:3628724/150‑1 (MQ=255)
                                                                                                                                 gatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                    <  2:648847/150‑1 (MQ=255)
                                                                                                                                 gatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                    <  1:1461685/150‑1 (MQ=255)
                                                                                                                                 gatCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCga                    <  2:648823/150‑1 (MQ=255)
                                                                                                                                    cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGaa                 >  1:2440778/1‑150 (MQ=255)
                                                                                                                                    cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGaa                 >  1:4014780/1‑150 (MQ=255)
                                                                                                                                    cGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCAGTAGATTACCCAATAGTACTTGCTATCgagac                 >  1:2440824/1‑149 (MQ=255)
                                                                                                                                     gAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAAt                >  1:1063016/1‑150 (MQ=255)
                                                                                                                                        aaaCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCACTAGTACTTGCTATCGAGAATAAg             >  1:2360409/1‑150 (MQ=255)
                                                                                                                                         aaCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAgg            >  2:3404668/1‑150 (MQ=255)
                                                                                                                                         aaCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAgg            >  2:545542/1‑150 (MQ=255)
                                                                                                                                         aaCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAgg            <  1:3458198/150‑1 (MQ=255)
                                                                                                                                         aaCCGAAGCCGTCAGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAgg            >  1:1269382/1‑150 (MQ=255)
                                                                                                                                         aaCCGAAGCCGTCAGGCTTCGGTTTTTTATGCCCTGCATTCGGTGCAGTGAACCACCAGCAACACCCTCACCTATTTGTCACATGGGACGAATGTTCCAACCCAGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAgg            >  1:2873116/1‑150 (MQ=255)
                                                                                                                                          aCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGa           <  2:1360972/150‑1 (MQ=255)
                                                                                                                                          aCCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGa           >  1:442975/1‑150 (MQ=255)
                                                                                                                                           ccGAAGCCTGACGGCTTCGGTTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGa           <  2:2820341/150‑1 (MQ=255)
                                                                                                                                            cGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAg         >  1:1756511/1‑150 (MQ=255)
                                                                                                                                             gAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCCAGAATAAGGAAGa        <  1:3398502/150‑1 (MQ=255)
                                                                                                                                              aaGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGAc       <  1:2418043/150‑1 (MQ=255)
                                                                                                                                                gCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGc     <  2:2790068/150‑1 (MQ=255)
                                                                                                                                                  cTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCAt   <  1:158821/150‑1 (MQ=255)
                                                                                                                                                  cTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCAt   <  2:4068218/150‑1 (MQ=255)
                                                                                                                                                   tGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATg  >  1:538557/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GGTATCGTAAGCGGCGAAACCGACGGCATGCAGGCCTTCATGGGTGGCAAGCTGCGGGTTGAAGGCGACATGATGCTGTCGATGAAACTCAGCGAGCTGTTCCCGTCCTGAGGGCTGGGAATACCGGATGATCGAAAAACCGAAGCCTGACGGCTTCGGTTTTTTTTGCCCTGCATTCGGTGCAGTGATCCACCAGCAACACCCTCACCTATTTGGCACATGGGACGATTGTTCCAACCCTGGCGGACCATTAGATTACCCAATAGTACTTGCTATCGAGAATAAGGAAGACGCATG  >  NC_002947/4431903‑4432199

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