Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,623 A→G 8.3% D191G (GAC→GGC)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,6230AG 8.3% 165.9 / 4.8 60D191G (GAC→GGC) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base G (3/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

CTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGC  >  NC_002947/3908476‑3908771
                                                                                                                                                   |                                                                                                                                                    
cTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGAc                                                                                                                                                     <  5:154834/149‑1 (MQ=255)
   aTCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTg                                                                                                                                                  >  3:58647/1‑149 (MQ=255)
    tCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTgg                                                                                                                                                 >  8:48070/1‑149 (MQ=255)
       aCTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGATCTGGccc                                                                                                                                              <  3:247573/149‑1 (MQ=255)
       aCTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGccc                                                                                                                                              >  2:143748/1‑149 (MQ=255)
        cTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCg                                                                                                                                             >  3:322832/1‑149 (MQ=255)
           cagcaACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                                                                                                            <  4:58647/147‑1 (MQ=255)
           cagcaACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                                                                                                            <  3:93751/147‑1 (MQ=255)
              caACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACa                                                                                                                                       <  1:143748/149‑1 (MQ=255)
                aCCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGt                                                                                                                                     >  3:228344/1‑149 (MQ=255)
                 cccGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGtt                                                                                                                                    <  7:48070/149‑1 (MQ=255)
                  ccGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGtt                                                                                                                                    >  4:305324/1‑148 (MQ=255)
                       cagcGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGc                                                                                                                              <  6:35586/149‑1 (MQ=255)
                           gcCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGc                                                                                                                          <  5:100343/149‑1 (MQ=255)
                              aGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGcccc                                                                                                                       <  3:192405/149‑1 (MQ=255)
                               gCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCg                                                                                                                      >  7:212206/1‑149 (MQ=255)
                                            cTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGGCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGAccc                                                                                                         >  7:119881/1‑149 (MQ=255)
                                               gCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTgg                                                                                                                                                 <  5:259778/106‑1 (MQ=255)
                                               gCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTgg                                                                                                                                                 >  6:259778/1‑106 (MQ=255)
                                               gCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCt                                                                                                        <  4:228344/147‑1 (MQ=255)
                                                cccGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAg                                                                                                     >  6:262176/1‑149 (MQ=255)
                                                 ccGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCggg                                                                                                                                           <  7:132751/110‑1 (MQ=255)
                                                 ccGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCggg                                                                                                                                           >  8:132751/1‑110 (MQ=255)
                                                 ccGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAgg                                                                                                    <  8:291037/149‑1 (MQ=255)
                                                  cgcgAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGc                                                                                                   >  2:155699/1‑149 (MQ=255)
                                                    cgAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCt                                                                                                 <  3:305324/149‑1 (MQ=255)
                                                          cgacgaTGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCg                                                                                           >  6:77599/1‑149 (MQ=255)
                                                           gacgaTGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGt                                                                                          <  5:77599/149‑1 (MQ=255)
                                                                  cAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAag                                                                                   <  1:155699/149‑1 (MQ=255)
                                                                      gcTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAagcagc                                                                               >  6:194344/1‑149 (MQ=255)
                                                                          gTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                                                                                                            >  8:165693/1‑84 (MQ=255)
                                                                          gTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCgg                                                                                                                                            <  7:165693/84‑1 (MQ=255)
                                                                          gTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGcc                                                                                                              <  6:319580/114‑1 (MQ=255)
                                                                            gcggcgCCTGAGCAACGCCCTGCAGCGCCCCTGCCTTGAAACGGGCGACTGGCTGTCACATCAACACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTga                                                                         <  7:71347/149‑2 (MQ=255)
                                                                              ggcgCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTAt                                                                       <  6:133112/149‑1 (MQ=255)
                                                                               gcgcCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATg                                                                      >  7:26987/1‑149 (MQ=255)
                                                                                      aGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGccc                                                                <  5:262176/148‑1 (MQ=255)
                                                                                         aaCGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAACACCCCGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCACAGGCCTATGCCGAGCTCACGCAGCGCCAGTATGGCGCCccgcc                                                            >  4:264734/1‑149 (MQ=255)
                                                                                                tGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGccc                                                     >  3:128060/1‑149 (MQ=255)
                                                                                                 gCAGCGCCCCTGCCTTGAAACGGTCGACTCGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGtgt                                                                       >  4:28414/1‑128 (MQ=255)
                                                                                                  cAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGccccc                                                   <  8:119881/149‑1 (MQ=255)
                                                                                                    gcgcCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCg                                                  >  3:44623/1‑148 (MQ=255)
                                                                                                        cccTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGcc                                             >  6:213972/1‑149 (MQ=255)
                                                                                                         ccTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCcg                                            >  4:201029/1‑149 (MQ=255)
                                                                                                               ttGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCAc                                       >  1:291671/1‑148 (MQ=255)
                                                                                                                          gACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAAt                            <  8:26987/148‑1 (MQ=255)
                                                                                                                               gCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaa                       >  7:173118/1‑148 (MQ=255)
                                                                                                                               gCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaa                       >  7:44478/1‑148 (MQ=255)
                                                                                                                                 tGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaacaa                    >  6:1943/1‑149 (MQ=255)
                                                                                                                                   tCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGaacaa                    <  5:194344/147‑1 (MQ=255)
                                                                                                                                            ttcccGGGCCAGGTCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCCGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCa         >  1:284311/3‑149 (MQ=255)
                                                                                                                                            gACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCa         <  3:201029/149‑1 (MQ=255)
                                                                                                                                              cccGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGa       >  5:124856/1‑149 (MQ=255)
                                                                                                                                              cccGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGa       >  4:47839/1‑149 (MQ=255)
                                                                                                                                               ccGGGCCAGGTCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCcgc                                           >  2:264472/1‑112 (MQ=255)
                                                                                                                                               ccGGGCCAGGTCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCcgc                                           <  1:264472/112‑1 (MQ=255)
                                                                                                                                               ccGGGCCAGGTCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCg            >  3:12818/1‑143 (MQ=255)
                                                                                                                                               ccGGGCCAGGTCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCg            <  4:12818/143‑1 (MQ=255)
                                                                                                                                               ccGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGAc      >  4:60622/1‑149 (MQ=255)
                                                                                                                                                  gACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCgg   >  6:305100/1‑149 (MQ=255)
                                                                                                                                                   aCCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGc  <  8:212206/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAGCAACGACCTGCAGCGCCCCTGCCTTGAAACGGTCGACTGGCTGTCACATCAAGACCCGGACCTGGCCCGGGAACAGTTCAAGGCCTGCCCCGCCAACGCCGACCCTCAGGCCTATGCCGTGCTCAAGCAGCGCCTGTATGGCGCCCCGCCGCAGCCCCCGCAGCCGCGCACCGCCGCAGAATGGGAACAACGCTATCGGCAGACCGGC  >  NC_002947/3908476‑3908771

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: ‑