Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,768,092 A→G 5.5% I711V (ATC→GTC)  PP_1577 → lambda family tail tape meausure protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,768,0920AG 5.5% 220.2 / 2.1 73I711V (ATC→GTC) PP_1577lambda family tail tape meausure protein
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (38/35)
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.32e-01

GCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTAGAGGTCC  >  NC_002947/1767948‑1768239
                                                                                                                                                |                                                                                                                                                   
gCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATccc                                                                                                                                                 >  7:283944/1‑149 (MQ=255)
 cATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCcg                                                                                                                                                <  8:329444/149‑1 (MQ=255)
 cATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCcg                                                                                                                                                >  8:88609/1‑149 (MQ=255)
  aTTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCcgc                                                                                                                                               >  2:54315/1‑149 (MQ=255)
          cgATGATCGGGGTTCCCTTCATGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAggg                                                                                                                                        >  6:217715/1‑148 (MQ=255)
                      ttCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGctgctg                                                                                                                           <  3:164680/149‑1 (MQ=255)
                         ccTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGa                                                                                                                         <  6:359103/148‑1 (MQ=255)
                              gTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGATAACATCCCGCGAGAGGGTACCTGGCTGCTGGATggcg                                                                                                                    <  1:330730/148‑1 (MQ=255)
                               tGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATggcggc                                                                                                                  <  4:213123/149‑1 (MQ=255)
                                   gCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagc                                                                                                              <  8:223807/149‑1 (MQ=255)
                                   gCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGCGTTGGCTCGGCCCTTTCTGGGCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGgcgagc                                                                                                              <  7:277396/149‑1 (MQ=255)
                                      ggCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGcgagcgag                                                                                                           >  4:72281/1‑149 (MQ=255)
                                        cAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTc                                                                                                         <  5:51945/149‑1 (MQ=255)
                                         aaTTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTcc                                                                                                        >  3:206198/1‑149 (MQ=255)
                                              cACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGc                                                                                                   >  4:109309/1‑149 (MQ=255)
                                              cACCTGGCGCTGCCATTGCAGCCGCAGCCGCAACAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGc                                                                                                   >  1:39296/1‑149 (MQ=255)
                                                 cTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCt                                                                                                >  8:25644/1‑149 (MQ=255)
                                                 cTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCt                                                                                                >  1:52339/1‑149 (MQ=255)
                                                 cTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGACCCTATCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCt                                                                                                >  1:71579/1‑149 (MQ=255)
                                                      gcTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTgcg                                                                                             >  6:198429/1‑147 (MQ=255)
                                                            aTTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCg                                                                                     >  6:127708/1‑149 (MQ=255)
                                                                cagccgcagccgCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGa                                                                                  <  7:189563/148‑1 (MQ=255)
                                                                cagccgcagccgCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGAc                                                                                 <  2:9881/149‑1 (MQ=255)
                                                                cagccgcagccgCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGAc                                                                                 <  3:72281/149‑1 (MQ=255)
                                                                          cgcagcAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCt                                                                       <  5:198429/149‑1 (MQ=255)
                                                                           gcagcaGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTa                                                                      >  3:151369/1‑149 (MQ=255)
                                                                             agcagGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTAcc                                                                    >  3:342007/1‑149 (MQ=255)
                                                                                 ggCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTg                                                                  >  6:329393/1‑147 (MQ=255)
                                                                                       cTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGc                                                          <  7:88609/149‑1 (MQ=255)
                                                                                        tgatgaCTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGc                                                          >  5:251252/1‑148 (MQ=255)
                                                                                                 gggTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGtt                                                <  5:329393/149‑1 (MQ=255)
                                                                                                   gTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGa                                              <  4:151369/149‑1 (MQ=255)
                                                                                                   gTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGa                                              <  1:115128/149‑1 (MQ=255)
                                                                                                    ttGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             <  6:30030/149‑1 (MQ=255)
                                                                                                    ttGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             <  2:39296/149‑1 (MQ=255)
                                                                                                      ggCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGc                                           <  4:204706/149‑1 (MQ=255)
                                                                                                        cTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCa                                         >  7:32129/1‑149 (MQ=255)
                                                                                                         tCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGATGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCAc                                        >  1:279225/1‑149 (MQ=255)
                                                                                                         tCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCAc                                        >  7:6533/1‑149 (MQ=255)
                                                                                                         tCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCAc                                        >  6:329102/1‑149 (MQ=255)
                                                                                                             cccTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTc                                     >  3:324153/1‑148 (MQ=255)
                                                                                                              ccTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCTGGCCAACAGCGTTGATGCCACGTcc                                    >  6:212961/1‑148 (MQ=255)
                                                                                                              ccTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTcc                                    >  5:207542/1‑148 (MQ=255)
                                                                                                              ccTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTcc                                    >  5:354518/1‑148 (MQ=255)
                                                                                                               cTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCaa                                  <  6:354518/149‑1 (MQ=255)
                                                                                                               cTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACATCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTcc                                    >  4:252552/1‑147 (MQ=255)
                                                                                                               cTTTCTGGCCAGGCTCACGGGGGTATGGAAAACATCCCGCGAGAGGGGACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCaa                                  <  8:59079/149‑1 (MQ=255)
                                                                                                                tttCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCAAACAGCGTTGATGCCACGTCCAAg                                 >  6:184648/1‑149 (MQ=255)
                                                                                                                       ccAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCtggtg                          <  8:32129/149‑1 (MQ=255)
                                                                                                                       ccAGGCTCACGGTGGTATGGAAAACATCACGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAAGCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCtggtg                          >  6:356045/1‑149 (MQ=255)
                                                                                                                           gCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCtggtggtgg                      <  4:324153/149‑1 (MQ=255)
                                                                                                                              cACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACCACCAGGCTTGGACGTGgcat                   <  2:52339/149‑3 (MQ=255)
                                                                                                                                     gtttGGAAAACATCCCGCGAGAGGGTACGTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAAGCGGGAGCTCACCGGCTACCTGCAGCTGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAAc            <  2:71579/146‑1 (MQ=255)
                                                                                                                                     gtATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTcc                                    >  6:106270/1‑125 (MQ=255)
                                                                                                                                     gtATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTcc                                    <  5:106270/125‑1 (MQ=255)
                                                                                                                                       aTGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACg           >  3:234129/1‑148 (MQ=255)
                                                                                                                                        tGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACg           >  1:207732/1‑147 (MQ=255)
                                                                                                                                        tGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGGCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTa         <  6:78292/149‑1 (MQ=255)
                                                                                                                                         ggAAAACGTCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             >  2:95038/1‑112 (MQ=255)
                                                                                                                                         ggAAAACGTCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             >  2:90058/1‑112 (MQ=255)
                                                                                                                                         ggAAAACGTCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             <  1:95038/112‑1 (MQ=255)
                                                                                                                                         ggAAAACGTCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGAt                                             <  1:90058/112‑1 (MQ=255)
                                                                                                                                          gAAAACATGCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGgcgc                   <  2:44443/137‑1 (MQ=255)
                                                                                                                                          gAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTag        <  1:54315/148‑1 (MQ=255)
                                                                                                                                          gAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTa         >  3:343325/1‑147 (MQ=255)
                                                                                                                                          gAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTa         >  3:337707/1‑147 (MQ=255)
                                                                                                                                          gAAAACATCCCGCGAGAGGGTACCTGGCTGCTGAATCGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATTCCACGACCAAGCCTGGTGGTGgcgc                   >  1:44443/1‑137 (MQ=255)
                                                                                                                                            aaaCATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAAc            >  4:337016/1‑142 (MQ=255)
                                                                                                                                            aaaCATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAAc            <  3:337016/142‑1 (MQ=255)
                                                                                                                                             aaCATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGcc                                                         <  6:270074/96‑1 (MQ=255)
                                                                                                                                             aaCATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGcc                                                         >  5:270074/1‑96 (MQ=255)
                                                                                                                                              aCATCCCGCGAGAGGGTACCTGCCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCGGCGCACAACCGGGACCTCCCCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCCAGCCTGGTGGCGGCGCCATGAACGTAGAggcc   <  4:206198/149‑3 (MQ=255)
                                                                                                                                               cATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTAGAGGTcc  <  5:127708/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCATTCACCGCGATGATCGGGGTTCCCTTCGTGGGGCCGGCAATTGCACCTGGCGCTGCCATTGCAGCCGCAGCCGCAGCAGGCGTGCTGATGACTGGGGTTGGCTCGGCCCTTTCTGGCCAGGCTCACGGTGGTATGGAAAACATCCCGCGAGAGGGTACCTGGCTGCTGGATGGCGGCGAGCGAGTCCTCAGCCCTGCGCAGAACCGGGACCTCACCGGCTACCTGCAGCGGGCCAACAGCGTTGATGCCACGTCCAAGCCTGGTGGTGGCGCCATGAACGTAGAGGTCC  >  NC_002947/1767948‑1768239

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: