Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,951,074 Δ1 bp 7.6% coding (426/642 nt) PP_3485 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,951,0740T. 7.6% 622.9 / 43.6 183coding (426/642 nt)PP_3485hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (80/89);  new base . (5/9);  total (85/98)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 188 total aligned reads displayed.

TCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGA  >  NC_002947/3950925‑3951223
                                                                                                                                                      |                                                                                                                                                      
tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                                                                                                        >  2:4124503/1‑150 (MQ=255)
tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                                                                                                        >  1:2157290/1‑150 (MQ=255)
tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                                                                                                        <  2:3836708/150‑1 (MQ=255)
tCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGAGTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                                                                                                        >  2:2122878/1‑150 (MQ=255)
  ccGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgc                                                                                                                                                      >  1:645313/1‑150 (MQ=255)
   cGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgcg                                                                                                                                                     >  1:2931127/1‑150 (MQ=255)
       gCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGa                                                                                                                                                 <  1:305059/150‑1 (MQ=255)
           gATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAg                                                                                                                                            <  1:3287903/150‑1 (MQ=255)
            aTGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                                           <  2:3191408/150‑1 (MQ=255)
             tGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgaggg                                                                                                                                          <  1:4037251/150‑6 (MQ=38)
                 cGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTa                                                                                                                                      <  2:1291462/150‑1 (MQ=255)
                   ttGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTaca                                                                                                                                    <  1:2929789/150‑1 (MQ=255)
                   ttGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTaca                                                                                                                                    <  2:1570402/150‑1 (MQ=255)
                      cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGCgaagtgga                                                                                                                                 <  1:2463667/150‑8 (MQ=37)
                       ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTg                                                                                                                                >  2:2264233/1‑150 (MQ=255)
                       ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTg                                                                                                                                >  1:2200733/1‑150 (MQ=255)
                       ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGCgaagtgga                                                                                                                                 <  2:527449/150‑8 (MQ=37)
                        cGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGc                                                                                                                               <  2:4483260/150‑1 (MQ=255)
                          aGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTg                                                                                                                             <  1:721069/150‑1 (MQ=255)
                           gCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTgg                                                                                                                            <  2:3230644/150‑1 (MQ=255)
                            ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTggg                                                                                                                           <  1:1788413/150‑1 (MQ=255)
                             cAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgc                                                                                                                          <  2:1858695/150‑1 (MQ=255)
                              acggccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCCGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                                         <  1:272298/148‑1 (MQ=255)
                              aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                                         >  1:4230782/1‑150 (MQ=255)
                              aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                                         >  2:435119/1‑150 (MQ=255)
                              aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                                         >  2:2578216/1‑150 (MQ=255)
                                 gccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgccgcc                                                                                                                      <  1:1230797/150‑1 (MQ=255)
                                 gccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgccgcc                                                                                                                      <  2:3860276/150‑1 (MQ=255)
                                 gccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgccgcc                                                                                                                      <  1:3015644/150‑1 (MQ=255)
                                   caggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCAt                                                                                                                    <  2:3073157/150‑1 (MQ=255)
                                         gtgtTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATg                                                                                                              <  2:2191260/150‑1 (MQ=255)
                                             tCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                          <  2:2057113/150‑1 (MQ=255)
                                             tCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                          <  2:501129/150‑1 (MQ=255)
                                             tCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                          <  1:1374163/150‑1 (MQ=255)
                                               gcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                        <  1:2609598/150‑1 (MQ=255)
                                                     gACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTc                                                                                                  >  1:1152808/1‑150 (MQ=255)
                                                      aCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                                                 <  1:1463721/150‑1 (MQ=255)
                                                       ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                                                <  1:3377044/150‑1 (MQ=255)
                                                       ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                                                <  1:3377168/150‑1 (MQ=255)
                                                       ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                                                <  1:2287141/150‑1 (MQ=255)
                                                        cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                                               >  2:3722166/1‑150 (MQ=255)
                                                        cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                                               >  2:1446904/1‑150 (MQ=255)
                                                        cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                                               >  1:3561463/1‑150 (MQ=255)
                                                         aaTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                                              <  1:3817260/150‑1 (MQ=255)
                                                         aaTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                                              >  1:3792521/1‑150 (MQ=255)
                                                         aaTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                                              >  1:1357876/1‑150 (MQ=255)
                                                           tGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                        <  2:3570341/138‑1 (MQ=255)
                                                           tGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTa                                                                                            >  1:3621602/1‑150 (MQ=255)
                                                              aTTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTg                                                                                         <  1:2802551/150‑1 (MQ=255)
                                                                  tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGcc                                                                                     >  1:3923492/1‑150 (MQ=255)
                                                                  tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGcc                                                                                     >  1:3629614/1‑150 (MQ=255)
                                                                     cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTg                                                                                  >  1:2777010/1‑150 (MQ=255)
                                                                      tttGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                                 >  1:2696277/1‑150 (MQ=255)
                                                                      tttGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                                 >  2:4061768/1‑150 (MQ=255)
                                                                      tttGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                                 >  2:4061850/1‑150 (MQ=255)
                                                                       ttGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCTGGTGTATTCAGTGTACTGCGCCCTGTc                                                                                <  2:1717192/150‑1 (MQ=255)
                                                                          gTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGaa                                                                             <  1:2874236/150‑1 (MQ=255)
                                                                          gTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGaa                                                                             <  2:1680182/150‑1 (MQ=255)
                                                                           ttcttcGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCTCTTCGCCCTCGCTCGGCGCTTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCTAAg                                                                            <  1:375360/150‑1 (MQ=255)
                                                                           ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAg                                                                            <  1:173852/150‑1 (MQ=255)
                                                                           ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAg                                                                            >  2:2320496/1‑150 (MQ=255)
                                                                           ttCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAg                                                                            >  2:719058/1‑150 (MQ=255)
                                                                            tCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAATTGCCTCCACTTCGCCCTCGCTGGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGc                                                                           <  1:380310/150‑1 (MQ=255)
                                                                                    ccTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCctct                                                                   <  1:2255512/150‑1 (MQ=255)
                                                                                    ccTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCctct                                                                   <  1:3418801/150‑1 (MQ=255)
                                                                                     cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTg                                                                  <  2:242740/150‑1 (MQ=255)
                                                                                      ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                 >  2:2764965/1‑150 (MQ=255)
                                                                                      ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                 <  2:1274851/150‑1 (MQ=255)
                                                                                        gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCa                                                               <  1:3854794/150‑1 (MQ=255)
                                                                                                 ccTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCttt                                                      >  1:4408070/1‑150 (MQ=255)
                                                                                                 ccTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCttt                                                      >  1:1124387/1‑150 (MQ=255)
                                                                                                  cTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                     <  2:3389186/150‑1 (MQ=255)
                                                                                                  cTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                     <  2:1214048/150‑1 (MQ=255)
                                                                                                  cTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                     <  1:163806/150‑1 (MQ=255)
                                                                                                     gTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                                  <  1:1076022/150‑1 (MQ=255)
                                                                                                           aTCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                            <  2:3305508/150‑1 (MQ=255)
                                                                                                           aTCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                            >  2:777643/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:107865/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:1818310/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:2065369/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:2320236/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:4117290/1‑150 (MQ=255)
                                                                                                            tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                           >  1:855412/1‑150 (MQ=255)
                                                                                                             cAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCg                                          >  2:624738/1‑150 (MQ=255)
                                                                                                                         ccACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCGCTGTCGAAGCCCGACTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATg                              >  2:2722746/1‑150 (MQ=255)
                                                                                                                          cACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGc                             <  2:1896862/150‑1 (MQ=255)
                                                                                                                           aCTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGcc                            >  2:2703968/1‑150 (MQ=255)
                                                                                                                           aCTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGcc                            <  2:2846599/150‑1 (MQ=255)
                                                                                                                            cTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCg                           >  1:923696/1‑150 (MQ=255)
                                                                                                                            cTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCg                           <  2:645313/150‑1 (MQ=255)
                                                                                                                                 cccTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCg                      >  1:1519607/1‑150 (MQ=255)
                                                                                                                                  ccTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGa                     >  1:1610156/1‑150 (MQ=255)
                                                                                                                                   cctgctcaGCGCAGTCGACGCGCCGAGCGAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGa                     >  2:3532348/4‑150 (MQ=255)
                                                                                                                                   cctgctcaGCGCAGTCGACGCGCCGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAt                    >  1:1500109/4‑150 (MQ=255)
                                                                                                                                      gctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtctc                 <  2:3932152/150‑1 (MQ=255)
                                                                                                                                       ctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCg                <  2:2363743/150‑1 (MQ=255)
                                                                                                                                        tcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGc               >  2:1298208/1‑150 (MQ=255)
                                                                                                                                           gcgcGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCaa            >  2:1074714/1‑150 (MQ=255)
                                                                                                                                                 tCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTgc       >  1:4391018/1‑150 (MQ=255)
                                                                                                                                                      tGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGa  <  2:1198652/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGA  >  NC_002947/3950925‑3951223

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