Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,385,286 G→T 6.3% P84T (CCA→ACA)  PP_3860 ← FluMu protein gp47

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,385,2860GT 6.3% 200.4 / 3.5 63P84T (CCA→ACA) PP_3860FluMu protein gp47
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base T (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

CGGCGCGCACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGC  >  NC_002947/4385139‑4385433
                                                                                                                                                   |                                                                                                                                                   
cGGCGCGCACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTgg                                                                                                                                                     <  7:209055/148‑1 (MQ=255)
    gcgcACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAg                                                                                                                                                >  8:346663/1‑149 (MQ=255)
    gcgcACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAg                                                                                                                                                >  1:389239/1‑149 (MQ=255)
      gcACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACGATACCGGCATCCCCCAAGGCGGTATTAGCGCCGCGGACGACAATGGGGCCTGCGGCGGGACCAgcc                                                                                                                                              >  2:156644/1‑149 (MQ=255)
       cACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAgccg                                                                                                                                             <  4:398079/149‑1 (MQ=255)
          cGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGTTCCAgccgcag                                                                                                                                          <  8:219096/149‑3 (MQ=255)
            gCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACTGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGTTCCAgccgcaggc                                                                                                                                        <  5:359004/149‑5 (MQ=255)
                   cgcCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCCGCGGCTGGAACAGCCGCCTTGCGccc                                                                                                                                   <  1:217314/147‑1 (MQ=255)
                    gcCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTc                                                                                                                                >  3:380643/1‑149 (MQ=255)
                          ccATCCACCAGTCGCGTCGTTACGGTCACGACGTAGCGGCTGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTaa                                                                                                                           <  3:87029/148‑1 (MQ=255)
                                 ccaGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCAt                                                                                                                   >  4:324252/1‑149 (MQ=255)
                                    gACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACa                                                                                                                                                 >  7:329457/1‑116 (MQ=255)
                                    gACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACa                                                                                                                                                 <  8:329457/116‑1 (MQ=255)
                                     aCGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCgg                                                                                                               <  6:168779/149‑1 (MQ=255)
                                      cgcgTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGc                                                                                                              <  1:115045/149‑1 (MQ=255)
                                            gtTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCAttt                                                                                                        <  2:195853/149‑1 (MQ=255)
                                             ttACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCAttt                                                                                                        >  4:204934/1‑148 (MQ=255)
                                               aCGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGTTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGACCTTAAGGCGCCATTCGGCCCATTTca                                                                                                      <  8:282913/148‑1 (MQ=255)
                                                    cACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAgcagca                                                                                                >  1:297501/1‑149 (MQ=255)
                                                     acgacgTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAgcagca                                                                                                >  8:386963/1‑148 (MQ=255)
                                                       gacgTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGCGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTcagc                                                                                                    >  6:286746/1‑142 (MQ=255)
                                                       gacgTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGAGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTcagc                                                                                                    <  5:286746/142‑1 (MQ=255)
                                                        acgTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCAtctc                                                                                            <  1:296347/149‑1 (MQ=255)
                                                               ggcggcCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATctct                                                                                           >  6:393487/1‑143 (MQ=255)
                                                               ggcggcCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATctct                                                                                           <  5:393487/143‑1 (MQ=255)
                                                                  ggcCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCTTCTCTTCATCGCAg                                                                                  <  7:346663/149‑1 (MQ=255)
                                                                         cTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCg                                                                           <  1:156644/149‑1 (MQ=255)
                                                                             tGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTg                                                                       >  4:150883/1‑149 (MQ=255)
                                                                               cTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTgg                                                                      <  4:185706/148‑1 (MQ=255)
                                                                                  gTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCa                                                                                        >  6:256516/1‑127 (MQ=255)
                                                                                  gTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCa                                                                                        <  5:256516/127‑1 (MQ=255)
                                                                                  gTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGa                                                                    <  3:324252/147‑1 (MQ=255)
                                                                                       ccATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTggg                                                                     <  4:284088/141‑1 (MQ=255)
                                                                                       ccATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTggg                                                                     >  3:284088/1‑141 (MQ=255)
                                                                                       ccATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCcgc                                                             <  2:297501/149‑1 (MQ=255)
                                                                                        cATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCcgcg                                                            <  2:389239/149‑1 (MQ=255)
                                                                                             cGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCa                                                                                        >  2:400972/1‑116 (MQ=255)
                                                                                             cGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCa                                                                                        <  1:400972/116‑1 (MQ=255)
                                                                                             cGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGc                                                                             <  5:238488/127‑1 (MQ=255)
                                                                                             cGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGc                                                                             >  6:238488/1‑127 (MQ=255)
                                                                                                  tCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTc                                                  <  4:380643/149‑1 (MQ=255)
                                                                                                   ccaccaAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATAGCAGCTCGCCGGGTGGGACTGCCGCGCGATCCACtac                                                 >  8:223532/1‑147 (MQ=255)
                                                                                                     accaAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCt                                                >  1:145543/1‑148 (MQ=255)
                                                                                                        aaGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTc                                                  >  4:297831/1‑143 (MQ=255)
                                                                                                        aaGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTc                                                  <  3:297831/143‑1 (MQ=255)
                                                                                                        aaGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATg                                            >  3:332479/1‑149 (MQ=255)
                                                                                                             gTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTgg                                                                      >  5:387821/1‑118 (MQ=255)
                                                                                                             gTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTgg                                                                      <  6:387821/118‑1 (MQ=255)
                                                                                                              ttattaGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTGCTGATGGCCATa                                      >  4:277870/1‑149 (MQ=255)
                                                                                                              ttattaGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATa                                      >  5:226405/1‑149 (MQ=255)
                                                                                                               tattaGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATa                                      <  6:226405/148‑1 (MQ=255)
                                                                                                                    gcgcCGCTGACGACAATGGGGCCTGCGGCTGGGACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGAcc                                >  1:404710/1‑149 (MQ=255)
                                                                                                                      gccgcTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATAc                                     <  7:21607/142‑1 (MQ=255)
                                                                                                                      gccgcTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATAc                                     >  8:21607/1‑142 (MQ=255)
                                                                                                                           tGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAggcg                         >  2:14456/1‑149 (MQ=255)
                                                                                                                               gacAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGgcgc                     >  3:23357/1‑149 (MQ=255)
                                                                                                                                  aaTGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCAc                  <  3:204934/149‑1 (MQ=255)
                                                                                                                                       ggCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAg              >  2:17956/1‑148 (MQ=255)
                                                                                                                                       ggCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAg              >  2:24276/1‑148 (MQ=255)
                                                                                                                                           tGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAccc         >  3:48188/1‑149 (MQ=255)
                                                                                                                                           ggCGGCTGTTCCAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAccc         >  5:6249/2‑149 (MQ=255)
                                                                                                                                           ggCGGCTGTTCCAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCACAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAccc         >  6:232405/2‑149 (MQ=255)
                                                                                                                                             cGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGTTGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCcg        <  7:223532/148‑1 (MQ=255)
                                                                                                                                                 tGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAg   >  3:342210/1‑149 (MQ=255)
                                                                                                                                                  ggAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGc  >  5:244633/1‑149 (MQ=255)
                                                                                                                                                   gAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAg   <  7:386963/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGCGCGCACCGGCACCTGCGCCGTACCATCCACCAGACGCGTCGTTACGGTCACGACGTAGCGGCGGCCGTCCTGATGCTGGTAGACCATACCGGCATCCACCAAGGCGTTATTAGCGCCGCTGACGACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGC  >  NC_002947/4385139‑4385433

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