Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,884,024 A→G 6.6% Q348R (CAA→CGA)  PP_3428 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,884,0240AG 6.6% 170.6 / 2.2 60Q348R (CAA→CGA) PP_3428TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AACTGGAAACCGTGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGCACCGC  >  NC_002947/3883877‑3884164
                                                                                                                                                   |                                                                                                                                            
aaCTGGAAACCGTGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCaa                                                                                                                                             <  7:152806/149‑1 (MQ=255)
     gAAACCGTGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAggg                                                                                                                                        >  5:131662/1‑149 (MQ=255)
     gAAACCGTGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGAGCCAGAACCCGGACGACCCTGCGCTGCAACAggg                                                                                                                                        >  7:263697/1‑149 (MQ=255)
      aaaCCGTGCTCAATCGGCAGCCAGCCAATCGCCTGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGc                                                                                                                                       <  7:200605/149‑1 (MQ=255)
            tGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAg                                                                                                                                 <  5:43653/149‑1 (MQ=255)
               tCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCGACAGGGCCTGTAgg                                                                                                                                >  7:184905/1‑147 (MQ=255)
               tCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCGACAGGGCCTGTAgg                                                                                                                                >  7:190124/1‑147 (MQ=255)
               tCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAgg                                                                                                                                >  7:241413/1‑147 (MQ=255)
               tCAATCGACAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAgg                                                                                                                                >  6:179842/1‑147 (MQ=255)
                cAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAg                                                                                                                             >  4:61580/1‑149 (MQ=255)
                   tCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGc                                                                                                                          >  3:259114/1‑149 (MQ=255)
                       cagccagcCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTgg                                                                                                                      <  3:45953/149‑1 (MQ=255)
                         gccagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGAt                                                                                                                    <  6:110475/149‑1 (MQ=255)
                         gccagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGAt                                                                                                                    <  6:18651/149‑1 (MQ=255)
                          ccagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGAtt                                                                                                                   <  8:263697/149‑1 (MQ=255)
                           cagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGc                                                                                                                               <  3:25741/136‑1 (MQ=255)
                           cagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGc                                                                                                                               >  4:25741/1‑136 (MQ=255)
                           cagccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGAtt                                                                                                                   >  6:350122/1‑148 (MQ=255)
                            agccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACGCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCa                                                                                                                 <  8:241413/149‑1 (MQ=255)
                             gccaATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAg                                                                                                                <  5:179842/149‑1 (MQ=255)
                               caATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGtt                                                                                                              <  1:158456/149‑1 (MQ=255)
                                     ccGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGTGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAAccc                                                                                                        >  4:252484/1‑149 (MQ=255)
                                      cGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAAcc                                                                                                         >  8:178277/1‑147 (MQ=255)
                                      cGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAAcc                                                                                                         <  7:178277/147‑1 (MQ=255)
                                         gCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAg                                                                                                    <  5:350122/149‑1 (MQ=255)
                                          cACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCGACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGc                                                                                                   <  8:184905/149‑1 (MQ=255)
                                          cACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCGACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGc                                                                                                   <  8:190124/149‑1 (MQ=255)
                                                     ctgATTCAGTACTGGAAAGTGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAgg                                                                                         <  3:252484/148‑1 (MQ=255)
                                                     ctgATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGGGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAgg                                                                                         <  2:281407/148‑1 (MQ=255)
                                                     ctgATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAgg                                                                                         <  2:189271/148‑1 (MQ=255)
                                                      tgATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGaa                                                                                                           <  1:33154/129‑1 (MQ=255)
                                                      tgATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGaa                                                                                                           >  2:33154/1‑129 (MQ=255)
                                                      tgATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGt                                                                                       >  1:41995/1‑149 (MQ=255)
                                                          tCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAg                                                                                   >  6:124322/1‑149 (MQ=255)
                                                            aGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGt                                                                                  <  6:248947/148‑1 (MQ=255)
                                                             gTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTgg                                                                                <  4:324679/149‑1 (MQ=255)
                                                             cagcagcaaaGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTgg                                                                                <  3:105483/142‑1 (MQ=255)
                                                                   gAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTAt                                                                            >  5:209219/1‑147 (MQ=255)
                                                                      agagTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGc                                                                       >  8:125660/1‑149 (MQ=255)
                                                                      agagTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGc                                                                       <  1:235232/149‑1 (MQ=255)
                                                                          tCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCAcc                                                                   <  7:125660/149‑1 (MQ=255)
                                                                              gCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTg                                                                >  4:36911/1‑148 (MQ=255)
                                                                                cAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAg                                                             <  8:294574/149‑1 (MQ=255)
                                                                                cAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAg                                                             <  3:36911/149‑1 (MQ=255)
                                                                                 aGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAg                                                             >  1:48568/1‑148 (MQ=255)
                                                                                  gTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCTCCAGGGTCCAGTGTGTATCAGGCCACCCTGAAGga                                                           <  6:131662/149‑1 (MQ=255)
                                                                                     gCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGgaga                                                         <  6:209219/148‑1 (MQ=255)
                                                                                     gCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGgaga                                                         <  3:161851/148‑1 (MQ=255)
                                                                                     gCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGagag                                                        >  5:263023/1‑149 (MQ=255)
                                                                                        gAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTc                                                      >  6:287751/1‑148 (MQ=255)
                                                                                        gAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTc                                                      <  5:287751/148‑1 (MQ=255)
                                                                                              gCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTgg                                               >  4:235444/1‑149 (MQ=255)
                                                                                               cAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTgg                                               <  4:259114/148‑1 (MQ=255)
                                                                                                aaGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGc                                              >  1:169387/1‑148 (MQ=255)
                                                                                                        ctcgcCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGtt                                     >  1:307826/1‑149 (MQ=255)
                                                                                                           gcCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTg                                    <  3:61580/147‑1 (MQ=255)
                                                                                                                         ccAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTg                    >  7:53993/1‑149 (MQ=255)
                                                                                                                                      gacCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGc       >  2:248900/1‑149 (MQ=255)
                                                                                                                                           tGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGCACCGc  >  6:128210/1‑149 (MQ=255)
                                                                                                                                                tGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGCAcc    <  3:90671/142‑1 (MQ=255)
                                                                                                                                                tGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGCAcc    >  4:90671/1‑142 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AACTGGAAACCGTGCTCAATCGGCAGCCAGCCAATCGCCGGGCACGGGTGCTGCTGATTCAGTACTGGAAAGAGTCCGGCCAGTTGCAGAACGTGCAAGTATCGCTCGCCGAGCTGGAGCGCCAGAACCCGGACGACCCTGCGCTGCAACAGGGCCTGTAGGCAGTGCCTGGATTCAGTTGAACCCTGAGCCTGCGCCAGGGTCCAGTGGGTATCAGGCCACCCTGAAGGAGAGTCGTTTTGGCCAGTTCGTTGGCGCTGGATGAACGTGCCTTGATCCTCGCACCGC  >  NC_002947/3883877‑3884164

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