Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,604,864 C→A 9.7% intergenic (+26/‑36) PP_2282 → / → PP_2283 minor capsid protein 10/tail tubular protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,604,8640CA 9.7% 221.8 / 9.9 71intergenic (+26/‑36)PP_2282/PP_2283minor capsid protein 10/tail tubular protein A
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base A (2/5);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGGCGATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTCCCGTAGCGGATGCCGC  >  NC_002947/2604722‑2605012
                                                                                                                                               |                                                                                                                                                    
cTGGCGATGGAGTCGGACTACGACATCTGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCtt                                                                                                                                                <  2:75908/149‑1 (MQ=255)
cTGGCGATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCCCACACAAGGAAGTCCtt                                                                                                                                                >  7:59209/1‑149 (MQ=255)
 tGGCGATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTtg                                                                                                                                               >  5:122345/1‑149 (MQ=255)
  ggCGATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTtgt                                                                                                                                              >  5:304434/1‑149 (MQ=255)
     gATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTtgtgtg                                                                                                                                           >  2:40799/1‑149 (MQ=255)
     gATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTtgtgtg                                                                                                                                           >  6:179251/1‑149 (MQ=255)
         gAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCCCACACAAGGAAGTCCTTGAGTGAGGc                                                                                                                                       <  8:59209/149‑1 (MQ=255)
           gTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCCCACACAAGGAAGTCCTTGAGTGAGGCta                                                                                                                                     <  6:255064/149‑2 (MQ=255)
           gTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCCCACACAAGGAAGTCCTTGAGTGAGGCta                                                                                                                                     <  5:35948/149‑2 (MQ=255)
                aCTACGACATCCGCCGTCAAGGCACCCTCATGGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCtttttt                                                                                                                                 >  1:226947/1‑149 (MQ=255)
                  tACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTc                                                                                                                              >  2:3173/1‑149 (MQ=255)
                     gACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTg                                                                                                                           >  5:268423/1‑149 (MQ=255)
                          ccgccgTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTggg                                                                                                                                         >  2:265024/1‑130 (MQ=255)
                          ccgccgTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTggg                                                                                                                                         <  1:265024/130‑1 (MQ=255)
                               gTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAAt                                                                                                                 >  8:263025/1‑149 (MQ=255)
                                  aaGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCACCCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGa                                                                                                               <  7:219001/148‑1 (MQ=255)
                                  aaGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGa                                                                                                               <  7:263025/148‑1 (MQ=255)
                                   aGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAg                                                                                                             <  2:226947/149‑1 (MQ=255)
                                      cACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGaca                                                                                                          <  8:296182/149‑1 (MQ=255)
                                       aCCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGacac                                                                                                         >  7:214857/1‑149 (MQ=255)
                                       aCCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGaca                                                                                                          >  4:126658/1‑148 (MQ=255)
                                         ccTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTc                                                                                                       <  7:180424/149‑1 (MQ=255)
                                           tCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTggg                                                                                                                                         <  5:171949/113‑1 (MQ=255)
                                           tCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTggg                                                                                                                                         >  6:171949/1‑113 (MQ=255)
                                             aTGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGccc                                                                                                   <  4:259018/149‑1 (MQ=255)
                                              attctCGACAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGccc                                                                                                   <  4:263945/145‑1 (MQ=255)
                                                    gCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGa                                                                                             <  5:179251/149‑1 (MQ=255)
                                                     ccAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGAc                                                                                            <  2:162454/149‑1 (MQ=255)
                                                      cAAGAACGCCATGGGCCACGGCATACTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGAcc                                                                                           >  3:347752/1‑149 (MQ=255)
                                                         gTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAg                                                                                        <  6:304434/149‑1 (MQ=255)
                                                                aTGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATg                                                                                 <  3:3746/149‑1 (MQ=255)
                                                                aTGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGCACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATg                                                                                 >  1:118519/1‑149 (MQ=255)
                                                                    gCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCg                                                                              >  6:183489/1‑148 (MQ=255)
                                                                     ccACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCtt                                                                                                                                                >  3:235378/1‑81 (MQ=255)
                                                                     ccACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCtt                                                                                                                                                <  4:235378/81‑1 (MQ=255)
                                                                     ccACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTAGTGACTAGAAAATGAAGACAAACGCCCAAACGACCGAGCAAGATGCCGt                                                                             >  8:56787/1‑148 (MQ=255)
                                                                        cGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAg                                                                                    >  4:175580/1‑138 (MQ=255)
                                                                        cGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAg                                                                                    <  3:175580/138‑1 (MQ=255)
                                                                        cGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAAc                                                                         <  5:252499/149‑1 (MQ=255)
                                                                        cGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAAc                                                                         <  1:40799/149‑1 (MQ=255)
                                                                             tcctccGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAAcatcat                                                                    <  8:214857/149‑1 (MQ=255)
                                                                                tccGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCtttttttt                                                                                                                               >  7:16368/1‑87 (MQ=255)
                                                                                tccGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCtttttttt                                                                                                                               <  8:16368/87‑1 (MQ=255)
                                                                                tccGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCt                                                                 >  8:212068/1‑149 (MQ=255)
                                                                                  cGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCt                                                                 <  6:268423/147‑1 (MQ=255)
                                                                                     ccAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTaaa                                                                          <  8:123548/135‑1 (MQ=255)
                                                                                     ccAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTaaa                                                                          >  7:123548/1‑135 (MQ=255)
                                                                                      catcagcCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCTTAAACATCATGCTCAATAc                                                           <  4:347752/146‑1 (MQ=255)
                                                                                         cagcCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCAt                                                        <  2:118519/149‑1 (MQ=255)
                                                                                             cGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGc                                                    <  3:126658/149‑1 (MQ=255)
                                                                                                      gCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCa                                           <  1:3173/149‑1 (MQ=255)
                                                                                                            aaCCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCaa                                      >  3:194988/1‑148 (MQ=255)
                                                                                                            aaCCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTACTTGTGTGGGGATTTTTTTTCGTGCCTAGAAAATGAATACACTCGCCCAAACGACCAAGATAGATGCCGTAAACATCATGATCAATACCATTGGCGAGTAGCCAGTCAAc                                     >  3:105636/1‑149 (MQ=255)
                                                                                                                     accaccAACCTAGCCTCACTCAAGGACTTCCTTTTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCTAGTCGCCAGTCAACTCCCTGgat                            <  7:56787/149‑1 (MQ=255)
                                                                                                                       caccaACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTATAAAATGACGACACTCGCCCAAAGGACCGAGCTATATGCCGTAAACATCATGCTCAATACCATTGGCGACTCGCCaatg                                        <  5:263027/135‑4 (MQ=255)
                                                                                                                       caccaACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCt                                <  4:293403/143‑1 (MQ=255)
                                                                                                                       caccaACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCt                                >  3:293403/1‑143 (MQ=255)
                                                                                                                       caccaACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGAGTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGCCCGAGCTAGATGCCGTAAACATCAAGCTCAATACCATTGCCTACTCGCCaatg                                        >  6:263027/1‑132 (MQ=255)
                                                                                                                         ccaACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGatgatg                        >  2:334577/1‑149 (MQ=255)
                                                                                                                         aaaaaaaaGCCCCACACAAGGAAGTCCTTGAGTGAGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGatgatg                        <  2:223821/142‑1 (MQ=38)
                                                                                                                          caACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATgg                       <  5:183489/149‑1 (MQ=255)
                                                                                                                             ccTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATa                                                            <  6:172331/109‑1 (MQ=255)
                                                                                                                             ccTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATa                                                            >  5:172331/1‑109 (MQ=255)
                                                                                                                                 gCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTccc                <  8:21320/149‑1 (MQ=255)
                                                                                                                                   ctcactcaAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTc                  <  2:341390/145‑1 (MQ=255)
                                                                                                                                   ctcactcaAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTc                  >  1:341390/1‑145 (MQ=255)
                                                                                                                                   ctcactcaAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTCCCGt              <  6:122345/149‑1 (MQ=255)
                                                                                                                                    ccacacaAGGAAGTCCTTGAGTGGGGCTTTTTTTTCGTGCCTAGAAAATGaa                                                                                                              <  8:200673/51‑1 (MQ=18)
                                                                                                                                    ccacacaAGGAAGTCCTTGAGTGGGGCTTTTTTTTCGTGCATAGAAAATGaa                                                                                                              >  7:200673/2‑52 (MQ=18)
                                                                                                                                          aaGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTCCCGTAGCGGAt       >  8:187403/1‑149 (MQ=255)
                                                                                                                                          aaGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTCCCGTAGCGGAt       >  1:201781/1‑149 (MQ=255)
                                                                                                                                               cTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTCGCGAGTCGCCAGTCAACTCCCTGGATGATCGCCTTCCCGTAGCGGATgccgc  >  8:107833/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTGGCGATGGAGTCGGACTACGACATCCGCCGTCAAGGCACCCTCATGGTCGCCAAGTACGCCATGGGCCACGGCATCCTCCGCCCAGCAGCCGCTGTAGAGCTGAAAACCGCCTAACCACCAACCTAGCCTCACTCAAGGACTTCCTTGTGTGGGGCTTTTTTTTCGTGCCTAGAAAATGAAGACACTCGCCCAAACGACCGAGCTAGATGCCGTAAACATCATGCTCAATACCATTGGCGAGTCGCCAGTCAACTCCCTGGATGATGGCCTTCCCGTAGCGGATGCCGC  >  NC_002947/2604722‑2605012

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: