Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,949,092 Δ1 bp 8.8% coding (442/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,0920C. 8.8% 172.8 / 10.4 57coding (442/642 nt)PP_3485hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base . (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTC  >  NC_002947/3948946‑3949236
                                                                                                                                                  |                                                                                                                                                   
gcgcgcCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcgc                                                                                                                                                   >  6:269542/1‑149 (MQ=255)
  gcgcctgcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcgctc                                                                                                                                                 >  3:191939/1‑149 (MQ=255)
        gcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcgctc                                                                                                                                                 <  4:328154/143‑1 (MQ=255)
        gcgcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCctcgctc                                                                                                                                                 >  3:328154/1‑143 (MQ=255)
          gcctCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTc                                                                                                                                       <  6:275242/149‑1 (MQ=255)
               ccGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTg                                                                                                                                  >  4:73349/1‑149 (MQ=255)
                            aGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgaggg                                                                                                                     <  3:101655/148‑6 (MQ=38)
                             gCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCg                                                                                                                   >  7:194278/1‑149 (MQ=255)
                                   cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                    <  1:171101/142‑1 (MQ=38)
                                   cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGCGAGGGc                                                                                                                    >  2:171101/1‑142 (MQ=38)
                                    ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACtt                                                                                                            >  7:235790/1‑149 (MQ=255)
                                     cGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgaggg                                                                                                                     <  6:238533/139‑6 (MQ=38)
                                     cGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTGCTCAGCGCAGTCGACGCGCCGAGcgaggg                                                                                                                     >  5:238533/1‑134 (MQ=38)
                                      gAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGc                                                                                                          <  3:73349/149‑1 (MQ=255)
                                        gCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCACTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTg                                                                                                        >  8:29941/1‑149 (MQ=255)
                                         ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTgg                                                                                                       <  7:29941/149‑1 (MQ=255)
                                         ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTgg                                                                                                       <  4:191939/149‑1 (MQ=255)
                                           aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgc                                                                                                     >  6:233854/1‑149 (MQ=255)
                                            aggccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                    <  3:283365/149‑1 (MQ=255)
                                                 aggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCAt                                                                                               <  8:235790/149‑1 (MQ=255)
                                                   gccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCAtct                                                                                             <  5:269542/149‑1 (MQ=255)
                                                     cGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAggg                                                                                                                     >  6:342880/1‑123 (MQ=255)
                                                     cGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAggg                                                                                                                     <  5:342880/123‑1 (MQ=255)
                                                             cgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                   <  5:233854/149‑1 (MQ=255)
                                                                  gACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                              >  1:42611/1‑149 (MQ=255)
                                                                  gACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                              >  3:235439/1‑149 (MQ=255)
                                                                    ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                            <  8:194278/149‑1 (MQ=255)
                                                                    ccAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                            >  8:142592/1‑149 (MQ=255)
                                                                     cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                           >  3:47614/1‑149 (MQ=255)
                                                                     cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                           >  2:199174/1‑149 (MQ=255)
                                                                       aTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                         <  7:142592/149‑1 (MQ=255)
                                                                             tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                          <  7:151333/142‑1 (MQ=255)
                                                                             tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                          >  8:151333/1‑142 (MQ=255)
                                                                                             cgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCg                                                   >  5:208385/1‑149 (MQ=255)
                                                                                              gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGc                                                  <  4:235439/149‑1 (MQ=255)
                                                                                              gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGc                                                  >  2:22292/1‑149 (MQ=255)
                                                                                              gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGc                                                  <  2:42611/149‑1 (MQ=255)
                                                                                              gcTCCTTGCAGTCTGTCCTGGCTGCCATCAAAACATGCCTCCACTTCGCCCGCGCTCGGCGCGGCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGc                                                  <  4:281502/149‑1 (MQ=255)
                                                                                                  cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCt                                                                                                         <  5:211324/90‑1 (MQ=255)
                                                                                                  cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCt                                                                                                         >  6:211324/1‑90 (MQ=255)
                                                                                                   ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                    <  6:347814/94‑1 (MQ=255)
                                                                                                   ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGgcc                                                                                                    >  5:347814/1‑94 (MQ=255)
                                                                                                   ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGccc                                                    <  6:199281/142‑1 (MQ=255)
                                                                                                   ttGAAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGccc                                                    >  5:199281/1‑142 (MQ=255)
                                                                                                        gTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGGCGTCGACTGCGCTGAGCAGGGCGAACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCATGTCGAAGCCCGCCACTGCCAg                                         >  1:261550/1‑149 (MQ=255)
                                                                                                              ccTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCtt                                  >  8:43636/1‑149 (MQ=255)
                                                                                                                  gTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCTACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGa                              <  7:316795/149‑1 (MQ=255)
                                                                                                                  gTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCg                               >  8:316795/1‑148 (MQ=255)
                                                                                                                     ggCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAAt                            <  6:208385/148‑1 (MQ=255)
                                                                                                                      gCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATaa                          >  6:331175/1‑149 (MQ=255)
                                                                                                                      gCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATaa                          >  8:328652/1‑149 (MQ=255)
                                                                                                                       cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGTGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAAc                         <  5:331175/149‑1 (MQ=255)
                                                                                                                           aaaaaaTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCg                     >  2:108590/1‑149 (MQ=255)
                                                                                                                             aaaaTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTg                                                                                                        >  5:42251/1‑64 (MQ=255)
                                                                                                                             aaaaTGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTg                                                                                                        <  6:42251/64‑1 (MQ=255)
                                                                                                                                    cTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                            <  4:52243/117‑1 (MQ=255)
                                                                                                                                     tCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                            >  8:34842/1‑116 (MQ=255)
                                                                                                                                      ccACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCTGGTGTATTCAg                                                                           <  5:252250/84‑1 (MQ=255)
                                                                                                                                      ccACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAg                                                                           >  6:252250/1‑84 (MQ=255)
                                                                                                                                       cACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATg         >  3:21597/1‑149 (MQ=255)
                                                                                                                                              cccTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc  <  2:261550/149‑1 (MQ=255)
                                                                                                                                                 tcgctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAAt          >  8:133164/1‑138 (MQ=255)
                                                                                                                                                 tcgctcgGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAAt          <  7:133164/138‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGCGCCTGCGCCTCCCGCTGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGCGTCGACTGCGCTGAGCAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTC  >  NC_002947/3948946‑3949236

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