Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,700,374 T→A 5.2% I284I (ATT→ATA) ‡ phaA → poly(3‑hydroxyalkanoate) polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,700,3740TA 5.2% 242.4 / 2.5 77I284I (ATT→ATA) ‡phaApoly(3‑hydroxyalkanoate) polymerase
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGG  >  NC_002947/5700227‑5700516
                                                                                                                                                    |                                                                                                                                              
ggAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTa                                                                                                                                               <  4:251747/149‑1 (MQ=255)
 gAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTAc                                                                                                                                              >  2:295854/1‑149 (MQ=255)
     agagCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgc                                                                                                                                          <  7:253384/149‑1 (MQ=255)
         ccTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgcagc                                                                                                                                       >  7:15000/1‑148 (MQ=255)
         ccTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagca                                                                                                                                      <  4:294158/149‑1 (MQ=255)
          cTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCaa                                                                                                                                     <  4:134965/149‑1 (MQ=255)
             gcgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagcacc                                                                                                                                    <  5:69471/147‑3 (MQ=255)
              cgcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGAc                                                                                                                                 <  1:103813/149‑1 (MQ=255)
               gcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGAc                                                                                                                                 >  1:179710/1‑148 (MQ=255)
               gcgcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGAc                                                                                                                                 <  2:179710/148‑1 (MQ=255)
                 gcTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTg                                                                                                                              <  6:215682/149‑1 (MQ=255)
                   ttCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgca                                                                                                                                         <  4:127856/136‑1 (MQ=255)
                   ttCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGgca                                                                                                                                         >  3:127856/1‑136 (MQ=255)
                   ttCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGCCGGTAATCGacagcacc                                                                                                                                    >  2:31269/1‑139 (MQ=255)
                        cctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCTCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACAt                                                                                                                         >  8:265739/1‑147 (MQ=255)
                        cctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATg                                                                                                                        <  1:264218/148‑1 (MQ=255)
                        cctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGc                                                                                                                       >  2:304723/1‑149 (MQ=255)
                        cctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGc                                                                                                                       >  2:151476/1‑149 (MQ=255)
                         ctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCt                                                                                                                      >  5:138055/1‑149 (MQ=255)
                         ctgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCt                                                                                                                      <  3:33252/149‑1 (MQ=255)
                          tgcGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTc                                                                                                                     >  5:138482/1‑149 (MQ=255)
                                         agACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGCTCGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCgg                                                                                                       <  1:238103/149‑1 (MQ=255)
                                           aCCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTgg                                                                                                    <  1:304723/149‑1 (MQ=255)
                                           aCCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTgg                                                                                                    >  5:215104/1‑149 (MQ=255)
                                                 atcatcaGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCa                                                                                               <  8:15000/148‑1 (MQ=255)
                                                  tcatcaGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAg                                                                                                                                    >  4:9955/1‑110 (MQ=255)
                                                  tcatcaGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAg                                                                                                                                    <  3:9955/110‑1 (MQ=255)
                                                    atcaGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTg                                                                                           >  2:78298/1‑149 (MQ=255)
                                                        gCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACGTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTAcc                                                                                       >  5:261865/1‑149 (MQ=255)
                                                            gcgcAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCAc                                                                                   >  4:228144/1‑149 (MQ=255)
                                                            gcgcAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGCGGCATCACTTGTACCGCa                                                                                    <  2:61765/148‑1 (MQ=255)
                                                              gcAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtg                                                                                 >  2:308528/1‑149 (MQ=255)
                                                               cAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGTAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtgg                                                                                <  8:211534/149‑1 (MQ=255)
                                                                  cccGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtggtgg                                                                             <  5:67468/149‑1 (MQ=255)
                                                                  cccGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACtggtgg                                                                             >  5:104481/1‑149 (MQ=255)
                                                                    cGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGc                                                                           >  4:275189/1‑149 (MQ=255)
                                                                             cccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgcc                                                                  <  8:185300/149‑1 (MQ=255)
                                                                              ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCg                                                                                                        >  4:235509/1‑110 (MQ=255)
                                                                              ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCg                                                                                                        <  3:235509/110‑1 (MQ=255)
                                                                              ccAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATgccg                                                                 >  5:15050/1‑149 (MQ=255)
                                                                                    tGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTg                                                           <  1:78298/149‑1 (MQ=255)
                                                                                      aaTGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTgg                                                          <  6:215104/148‑1 (MQ=255)
                                                                                       aTGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGc                                                         >  6:47355/1‑148 (MQ=255)
                                                                                            gCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGaa                                                    >  4:225127/1‑148 (MQ=255)
                                                                                                gTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTg                                                                                                     >  4:10510/1‑95 (MQ=255)
                                                                                                gTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAgg                                               >  1:58370/1‑149 (MQ=255)
                                                                                                 tCCACCTACATCGATGCTCTTAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAATGACCTGAACATGCTCGGCGCCTTCTCCGGTGGCATCACTTGTACCGCACTGGTGGGACACTATGCCGCCATTGGCGAGAACAAg                                                <  7:146845/147‑1 (MQ=255)
                                                                                                 tCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAg                                                <  6:104481/147‑1 (MQ=255)
                                                                                                 tCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGt                                              >  4:238055/1‑149 (MQ=255)
                                                                                                   cACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCAc                                                                        >  6:312024/1‑121 (MQ=255)
                                                                                                   cACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCAc                                                                        <  5:312024/121‑1 (MQ=255)
                                                                                                    aCCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACATGAACATGCTCGGCGCCTGCTCCGTTGGCATACCTTGTACCGCAATGGTGGGCCACTATGCCGCCATGGGCGAGAACATGGTCaa                                           >  4:208489/1‑149 (MQ=255)
                                                                                                     ccTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAAc                                          <  8:188434/149‑1 (MQ=255)
                                                                                                     ccTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAAc                                          <  3:275189/149‑1 (MQ=255)
                                                                                                        aCATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAgg                                               <  2:57325/141‑1 (MQ=255)
                                                                                                        aCATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAgg                                               >  1:57325/1‑141 (MQ=255)
                                                                                                        aCATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGcc                                       >  5:233758/1‑149 (MQ=255)
                                                                                                         cATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGc                                        <  7:265739/147‑1 (MQ=255)
                                                                                                              aTGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGAcc                                 <  3:238055/149‑1 (MQ=255)
                                                                                                                     gAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTgg                          >  8:143336/1‑149 (MQ=255)
                                                                                                                      aaaGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCctgctg                           <  1:308528/147‑1 (MQ=255)
                                                                                                                      aaaGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTgg                          <  3:228144/148‑1 (MQ=255)
                                                                                                                       aagaagCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCctgc                             >  6:160064/1‑144 (MQ=255)
                                                                                                                       aagaagCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCctgc                             <  5:160064/144‑1 (MQ=255)
                                                                                                                        agaagCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTTCTCCGTTGGCATCACTTGTACCGCACTTGTGGGCCACTATGCCGCCATTGTCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCa                       <  6:15050/149‑1 (MQ=255)
                                                                                                                           agCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGgcgc                                                                                                                >  1:225108/1‑57 (MQ=255)
                                                                                                                           agCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGgcgc                                                                                                                <  2:225108/57‑1 (MQ=255)
                                                                                                                              cGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCg                                                        <  6:17604/110‑1 (MQ=255)
                                                                                                                              cGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCg                                                        >  5:17604/1‑110 (MQ=255)
                                                                                                                              cGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGc                 <  1:295854/149‑1 (MQ=255)
                                                                                                                                        ggTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGAcacca        <  5:47355/149‑1 (MQ=255)
                                                                                                                                         ttgctgcCGGTAATCGACAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCCCTATGCCGCCATTGGGGAGAACAAGGTCAACGACCTGACCCTGCTGGTCAGCGTGATGGCcaccac       >  4:197661/2‑149 (MQ=255)
                                                                                                                                          tGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGAcaccac       >  7:292064/1‑148 (MQ=255)
                                                                                                                                             tGTCGATTACCGGCAGCAATGACCTGAACATGCGCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGTGCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATg   <  8:140952/149‑1 (MQ=255)
                                                                                                                                               tCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATgg  <  6:138482/148‑1 (MQ=255)
                                                                                                                                                  aTTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGt                                              >  1:232356/1‑101 (MQ=255)
                                                                                                                                                  aTTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGt                                              <  2:232356/101‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGAAAAGAGCCTGGCGCGCTTCTGCCTGCGCTCGCAGCAGCAGACCTTCATCATCAGCTGGCGCAACCCGACCAAGGCCCAGCGTGAATGGGGCCTGTCCACCTACATCGATGCGCTGAAAGAAGCCGTCGACGCGGTGCTGTCGATTACCGGCAGCAAGGACCTGAACATGCTCGGCGCCTGCTCCGGTGGCATCACTTGTACCGCACTGGTGGGCCACTATGCCGCCATTGGCGAGAACAAGGTCAACGCCCTGACCCTGCTGGTCAGCGTGCTGGACACCACCATGG  >  NC_002947/5700227‑5700516

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: