Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,949,101:1 +A 9.9% coding (433/642 nt) PP_3485 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,949,1011.A 9.9% 268.3 / 23.2 71coding (433/642 nt)PP_3485hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (31/33);  new base A (4/3);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTC  >  NC_002947/3948956‑3949242
                                                                                                                                                  |                                                                                                                                              
gcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTc                                                                                                                                             >  1:8691/1‑149 (MQ=255)
     ccGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACt                                                                                                                                         <  4:217221/148‑1 (MQ=255)
       gCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgcg                                                                                                                                      <  7:133765/149‑1 (MQ=255)
       gCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgcg                                                                                                                                      >  3:107650/1‑149 (MQ=255)
        cTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTgcgc                                                                                                                                     <  6:263344/149‑1 (MQ=255)
           cccGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCt                                                                                                                                    <  8:360499/147‑1 (MQ=255)
               aTGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAg                                                                                                                             <  6:309150/149‑1 (MQ=255)
                 gAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTTTCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGTCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                            <  8:192608/148‑1 (MQ=255)
                         cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGa                                                                                                                                           >  4:40710/1‑126 (MQ=255)
                         cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACt                                                                                                                   >  4:47830/1‑149 (MQ=255)
                         cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                          <  1:156245/142‑1 (MQ=38)
                         cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                          >  2:156245/1‑142 (MQ=38)
                          ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                            >  5:283367/1‑139 (MQ=255)
                          ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                            <  6:283367/139‑1 (MQ=255)
                             aGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                            <  4:209519/136‑1 (MQ=255)
                             aGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAgg                                                                                                                            >  3:209519/1‑136 (MQ=255)
                               ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                          <  4:197592/136‑1 (MQ=37)
                               ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                          >  3:197592/1‑136 (MQ=37)
                                 aaGGCCAGGCCGTGTTCGCGTATGACCATTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgag                                                                                                                             >  1:233624/1‑128 (MQ=37)
                                 aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgc                                                                                                           >  8:12153/1‑149 (MQ=255)
                                 aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgag                                                                                                                             <  2:233624/131‑4 (MQ=38)
                                               ttCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGcc                                                                                             >  2:129290/1‑149 (MQ=255)
                                               ttCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGcc                                                                                             <  5:20423/149‑1 (MQ=255)
                                                  gcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                           <  3:47830/148‑1 (MQ=255)
                                                   cgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                         <  1:129290/149‑1 (MQ=255)
                                                   cgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGGCCTGGTGGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                         <  7:340546/149‑1 (MQ=255)
                                                            aaTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                                <  7:12153/149‑1 (MQ=255)
                                                             aTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                               <  3:289983/149‑1 (MQ=255)
                                                               ggATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCGTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                           <  6:297442/135‑1 (MQ=255)
                                                               ggATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                           >  5:297442/1‑135 (MQ=255)
                                                                   tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTgc                                                                         >  5:357513/1‑149 (MQ=255)
                                                                     tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGc                                                                                                                          <  6:204147/98‑1 (MQ=255)
                                                                     tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGc                                                                                                                          >  5:204147/1‑98 (MQ=255)
                                                                       aCTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGcc                                                                      <  4:107650/148‑1 (MQ=255)
                                                                        cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCt                                                                    >  3:75719/1‑149 (MQ=255)
                                                                                  gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                <  5:190694/111‑1 (MQ=255)
                                                                                  gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                >  6:190694/1‑111 (MQ=255)
                                                                                    gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAg                                                                                                                                 >  6:366351/1‑77 (MQ=255)
                                                                                    gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAg                                                                                                                                 <  5:366351/77‑1 (MQ=255)
                                                                                      tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                         >  4:338775/1‑114 (MQ=255)
                                                                                      tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                         <  3:338775/114‑1 (MQ=255)
                                                                                      tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGCCCGCCAACTCATGACGGGACTATTCACTGTACTCAGCCCTGTCTAAGCCCGCct                                                      >  6:398958/1‑148 (MQ=255)
                                                                                        cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCctct                                                    >  6:75031/1‑149 (MQ=255)
                                                                                           gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                  <  1:356166/148‑1 (MQ=255)
                                                                                           gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                  <  4:75719/148‑1 (MQ=255)
                                                                                                   tCCTGGTTGGCATCAAAAAATGCCTCAACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCAGCCATCTCATGACGGGTGTATTCAGTGTGATGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCt                                         >  5:393432/1‑149 (MQ=255)
                                                                                                        gTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                      <  5:75031/147‑1 (MQ=255)
                                                                                                        gTTGGCATCAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATATCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                   <  5:398958/149‑1 (MQ=255)
                                                                                                             cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                                  <  6:344227/98‑1 (MQ=255)
                                                                                                             cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                                  >  5:344227/1‑98 (MQ=255)
                                                                                                             cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAAc                               >  1:186731/1‑149 (MQ=255)
                                                                                                               tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                             >  6:194724/1‑149 (MQ=255)
                                                                                                               tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                             >  8:364989/1‑149 (MQ=255)
                                                                                                                      aTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                      <  8:326682/133‑1 (MQ=255)
                                                                                                                      aTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                      >  7:326682/1‑133 (MQ=255)
                                                                                                                      aTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                       >  3:238933/1‑148 (MQ=255)
                                                                                                                      aTGCCTCCACTTCGCCCGAGCTCGGCGCGTCGACTTCGCTGAGCAGGGAGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                       <  4:238933/148‑1 (MQ=255)
                                                                                                                          cTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCa                   >  5:16328/1‑148 (MQ=255)
                                                                                                                              aCTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGCTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCACTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGc              >  8:157153/1‑149 (MQ=255)
                                                                                                                                ttCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCg            >  4:375205/1‑149 (MQ=255)
                                                                                                                                   gCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                   >  4:222925/1‑123 (MQ=255)
                                                                                                                                   gCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                   <  3:222925/123‑1 (MQ=255)
                                                                                                                                    cccTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc        <  2:186731/149‑1 (MQ=255)
                                                                                                                                     ccctgctccGCGCAGTCGACGCGCCGAGCGAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc        >  7:270665/5‑149 (MQ=255)
                                                                                                                                      ctctctcgGCGCGTCTACTGCGCTGAGCAGGGCGTACACTTGCTGTGCCGCCATCTCATGCAGGGTGTATTCAGTGTACTGCGCCCttt                                                                  <  8:355699/89‑3 (MQ=255)
                                                                                                                                      ctcgctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                  >  7:355693/1‑89 (MQ=255)
                                                                                                                                      ctcgctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                  <  8:355693/89‑1 (MQ=255)
                                                                                                                                      atcgatcgTCGCGTCGACTGCGCTGAGAAGTGCGTACACTTGCTGGGCCGCCATATCATGCCGGGTGTATTCAGTGTACTGAGCCCTGt                                                                  >  7:355699/2‑89 (MQ=255)
                                                                                                                                         gctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtct   >  8:300164/1‑149 (MQ=255)
                                                                                                                                         gctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAt     <  5:340439/147‑1 (MQ=255)
                                                                                                                                          ctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtctc  <  4:292790/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTC  >  NC_002947/3948956‑3949242

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