Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,385,412 C→T 6.8% A42T (GCC→ACC)  PP_3860 ← FluMu protein gp47

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,385,4120CT 6.8% 173.6 / 2.0 59A42T (GCC→ACC) PP_3860FluMu protein gp47
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

GACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGACTCCAGGTATCT  >  NC_002947/4385266‑4385555
                                                                                                                                                   |                                                                                                                                               
gacAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGgcgc                                                                                                                                               <  2:242021/149‑1 (MQ=255)
       gggCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTaa                                                                                                                                         <  5:179752/148‑1 (MQ=255)
        ggCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCTGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGa                                                                                                                                       <  2:316291/149‑1 (MQ=255)
        ggCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGa                                                                                                                                       <  5:15879/149‑1 (MQ=255)
          aaggCGGCTGTTCCAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAc                                                                                                                                     >  7:5355/4‑149 (MQ=255)
          aaggCGGCTGTTCCAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACCGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAc                                                                                                                                     >  3:265133/4‑149 (MQ=255)
           cTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAAcc                                                                                                                                    <  7:242453/149‑1 (MQ=255)
                    gAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGcc                                                                                                                           <  6:185784/149‑1 (MQ=255)
                    gAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGcc                                                                                                                           >  6:268270/1‑149 (MQ=255)
                    gAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGcc                                                                                                                           <  7:225742/149‑1 (MQ=255)
                    gAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGc                                                                                                                            <  5:114313/148‑1 (MQ=255)
                      aCAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGGTCCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCCCTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCa                                                                                                                          >  1:62972/1‑149 (MQ=255)
                                  cgcCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCag                                                                                                              <  6:213211/148‑1 (MQ=255)
                                   gcCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCcg                                                                                                                                  <  1:151415/127‑1 (MQ=255)
                                   gcCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCcg                                                                                                                                  >  2:151415/1‑127 (MQ=255)
                                   gcCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCagag                                                                                                            >  4:83444/1‑149 (MQ=255)
                                      cGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAgcgg                                                                                                         <  2:114961/149‑1 (MQ=255)
                                          cTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAgcggcg                                                                                                       >  3:105587/1‑147 (MQ=255)
                                                gcgcCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGt                                                                                               >  8:230909/1‑149 (MQ=255)
                                                gcgcCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGt                                                                                               >  3:195663/1‑149 (MQ=255)
                                                    cATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCGGGGTGGTACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCgg                                                                                           <  4:68551/149‑1 (MQ=255)
                                                         gg‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGCCCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCAt                                                                                      >  3:279861/1‑149 (MQ=255)
                                                              cATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTgcg                                                                                  <  6:217260/149‑1 (MQ=255)
                                                               aTTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTgcgc                                                                                 >  6:350091/1‑149 (MQ=255)
                                                               aTTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTgcgc                                                                                 <  5:368158/149‑1 (MQ=255)
                                                                     gcagcaTCTCTTCATCGCAGGTCGCCGGGTGGGACCGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCaaaa                                                                           <  4:265133/149‑1 (MQ=255)
                                                                      cagcaTCTCTGCATCGCAGGTCGACGGGTGTGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAAt                                                                          <  8:303856/149‑1 (MQ=255)
                                                                       agcaTCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATc                                                                         <  3:83444/149‑1 (MQ=255)
                                                                         caTCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATcc                                                                        >  3:339942/1‑148 (MQ=255)
                                                                          aTCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGATGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGa                                                                      >  8:274814/1‑149 (MQ=255)
                                                                             tctTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCtt                                                                                                                      <  4:304596/98‑1 (MQ=255)
                                                                             tctTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCtt                                                                                                                      >  3:304596/1‑98 (MQ=255)
                                                                                           cGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAAACCGATTGGGTACGCGTGAt                                                       >  8:309537/1‑147 (MQ=255)
                                                                                                 gTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCgg                                               <  4:105587/149‑1 (MQ=255)
                                                                                                   gggACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCggg                                              >  5:117947/1‑148 (MQ=255)
                                                                                                    ggACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCTAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCTGGc                                             <  7:230909/148‑1 (MQ=255)
                                                                                                    ggACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCggg                                              >  6:276354/1‑147 (MQ=255)
                                                                                                        tGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGt                                        >  2:85754/1‑149 (MQ=255)
                                                                                                          ccGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCg                                      >  5:127671/1‑149 (MQ=255)
                                                                                                          ccGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCg                                      <  4:195663/149‑1 (MQ=255)
                                                                                                             cgcgATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCa                                   <  5:276354/149‑1 (MQ=255)
                                                                                                                  tCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCa                              >  1:187816/1‑149 (MQ=255)
                                                                                                                     aCTCCTGATGGCCATACAGACCGTAGGCGGTGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAgcg                                                  >  4:135120/1‑126 (MQ=255)
                                                                                                                     aCTCCTGATGGCCATACAGACCGTAGGCGGTGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAgcg                                                  <  5:366940/126‑1 (MQ=255)
                                                                                                                     aCTCCTGATGGCCATACAGACCGTAGGCGGTGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAgcg                                                  >  6:366940/1‑126 (MQ=255)
                                                                                                                     aCTCCTGATGGCCATACAGACCGTAGGCGGTGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAgcg                                                  <  3:135120/126‑1 (MQ=255)
                                                                                                                       tCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCgggg                                                                                         >  4:95263/1‑85 (MQ=255)
                                                                                                                       tCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCgggg                                                                                         <  3:95263/85‑1 (MQ=255)
                                                                                                                          tGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCAt                      >  5:216694/1‑149 (MQ=255)
                                                                                                                            aTGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCAtt                     >  1:307838/1‑148 (MQ=255)
                                                                                                                               gCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGAt                  >  2:351883/1‑148 (MQ=255)
                                                                                                                               gCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATc                 <  2:307838/149‑1 (MQ=255)
                                                                                                                                ccATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGAt                  >  4:283877/1‑147 (MQ=255)
                                                                                                                                ccATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCa                <  3:283877/149‑1 (MQ=255)
                                                                                                                                ccATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCa                <  5:268270/149‑1 (MQ=255)
                                                                                                                                ccATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCa                <  6:117947/149‑1 (MQ=255)
                                                                                                                                ccATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGAGCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCa                <  2:62972/149‑1 (MQ=255)
                                                                                                                                  aTACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGa              >  5:318112/1‑149 (MQ=255)
                                                                                                                                       gACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGACTCCa         <  6:216694/149‑1 (MQ=255)
                                                                                                                                         ccGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGACTCCAgg       <  4:279861/149‑1 (MQ=255)
                                                                                                                                              ggcggcGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGACTCCAGGTATCt  <  4:339942/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GACAATGGGGCCTGCGGCTGGAACAGCCGCCTTGCGCCCGTCCTTAAGGCGCCATTCGG‑CCCATTTCAGCAGCATCTCTTCATCGCAGGTCGCCGGGTGGGACTGCCGCGCGATCCACTCCTGATGGCCATACAGACCGTAGGCGGCGCCACTAAGAACCCGCGAAGCCACCTTGGCGTCAGAGCGGCGCAGCGCGTCCGGGGCATTGCGCTCAAAATCCGATTGGGTACGCGTGATCAGCGCGGGCAGCGTCGGCATCTCATACGGCATTGATCAGACTCCAGGTATCT  >  NC_002947/4385266‑4385555

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