Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,399,522 A→T 6.8% I123N (ATT→AAT)  PP_3878 ← minor capsid protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,399,5220AT 6.8% 216.5 / 4.7 73I123N (ATT→AAT) PP_3878minor capsid protein C
Reads supporting (aligned to +/- strand):  ref base A (29/39);  new base T (3/2);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGAT  >  NC_002947/4399374‑4399670
                                                                                                                                                    |                                                                                                                                                    
gCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACaa                                                                                                                                                      >  8:87552/1‑149 (MQ=255)
 cAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAAt                                                                                                                                                     <  1:110871/149‑1 (MQ=255)
 cAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAAt                                                                                                                                                     <  1:25851/149‑1 (MQ=255)
  aaGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAACACATtg                                                                                                                                                    >  1:138849/1‑149 (MQ=255)
  aaGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAACACATtg                                                                                                                                                    >  5:337278/1‑149 (MQ=255)
   aGCAGAGGCATTCCGGTACCCCCCTGAGTACGCCATGATGTTAACCAGGCCGGTGATTGGTTTGATCTTGGGCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAAtgt                                                                                                                                                   <  4:304993/149‑1 (MQ=255)
    gCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAAtgtg                                                                                                                                                  >  6:229220/1‑149 (MQ=255)
     cAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAACACATtgtgc                                                                                                                                                 <  6:337278/149‑2 (MQ=255)
         cGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGa                                                                                                                                             <  3:338510/149‑1 (MQ=255)
         cGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAACACATTGTGCTCGa                                                                                                                                             <  2:138849/149‑1 (MQ=255)
          gCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTc                                                                                                                                               >  4:87622/1‑146 (MQ=255)
          gCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTc                                                                                                                                               <  3:87622/146‑1 (MQ=255)
             aTCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCAt                                                                                                                                         >  8:137428/1‑149 (MQ=255)
             aTCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCAt                                                                                                                                         <  4:218033/149‑1 (MQ=255)
              tCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATg                                                                                                                                        <  2:101347/149‑1 (MQ=255)
                aGGTAACCCCCTGAGTACGACATGAAGTTAACCAGGCCGTTGATTGGTTTTATCTTGGTCGCTGCCCGAATGTGCGCCGCCAGTGCAAAAGCGCCGAGGGCGCTGCCGCCTGGTCTGTCAATGTCGAGCACAATGTGTTCGACGATggg                                                                                                                                      <  6:295591/149‑1 (MQ=255)
                   tAATCCGCTGAGTACGGCATGAAGTTAACGCGGCCGGTGAGTTGGTTGTTCTTGGTCGCGGCACGAATGTCCGCCGCCAGTTCAAAATCGCCTACGGCGCTGCCTCCTGGGCTGTCAATGTCGAGCACAATGTATTCGACCATGGGGTc                                                                                                                                   <  2:319983/149‑1 (MQ=255)
                   tAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTc                                                                                                                                   <  2:246126/149‑1 (MQ=255)
                      cccccTGAGTACGCAATGAAGTTAACCTGCCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGCCGCCGCTGCCGCCCGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGt                                                                                                                                <  1:110700/149‑1 (MQ=255)
                       ccccTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGt                                                                                                                                >  6:131402/1‑148 (MQ=255)
                         ccTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGAt                                                                                                                             >  7:51693/1‑149 (MQ=255)
                             aGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGc                                                                                                                          <  8:298920/148‑1 (MQ=255)
                                  gCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTAtt                                                                                                                    >  3:61873/1‑149 (MQ=255)
                                  gCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTAt                                                                                                                     >  3:103174/1‑148 (MQ=255)
                                    cATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTAtt                                                                                                                    >  4:138565/1‑147 (MQ=255)
                                     aTGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAg                                                                                                                 >  6:3722/1‑149 (MQ=255)
                                      tGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAgg                                                                                                                <  5:192728/149‑1 (MQ=255)
                                             aaCCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGtt                                                                                                                                                >  5:211525/1‑110 (MQ=255)
                                             aaCCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGtt                                                                                                                                                <  6:211525/110‑1 (MQ=255)
                                             aaCCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTc                                                                                                         >  3:193047/1‑149 (MQ=255)
                                                cAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGgctgcgag                                                                                                        >  3:277201/1‑144 (MQ=255)
                                                   gCCGGTTATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATg                                                                                                                            <  7:212185/124‑1 (MQ=255)
                                                   gCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATg                                                                                                                            >  8:212185/1‑124 (MQ=255)
                                                    ccGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGccc                                                                                                  >  3:91666/1‑149 (MQ=255)
                                                    ccGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGccc                                                                                                  <  5:207055/149‑1 (MQ=255)
                                                      ggTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCtt                                                                                                <  7:87552/149‑1 (MQ=255)
                                                       gTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTc                                                                                               <  4:193047/149‑1 (MQ=255)
                                                        tGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCa                                                                                              >  8:321266/1‑149 (MQ=255)
                                                          aTTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCa                                                                                              >  2:172012/1‑147 (MQ=255)
                                                            tGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAATGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGc                                                                                          >  1:267874/1‑149 (MQ=255)
                                                             ggTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATg                                                                                                                            >  1:124723/1‑114 (MQ=255)
                                                             ggTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATg                                                                                                                            <  2:124723/114‑1 (MQ=255)
                                                             ggTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCGAGCAGCCCtgaataag                                                                                          <  4:277201/148‑8 (MQ=255)
                                                              gTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTg                                                                                        >  4:354513/1‑149 (MQ=255)
                                                                  gATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTc                                                                                     <  5:229220/148‑1 (MQ=255)
                                                                        ggTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATg                                                                                  <  3:138565/145‑1 (MQ=255)
                                                                         gTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGtctc                                                                             >  4:217536/1‑149 (MQ=255)
                                                                            gcgGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCg                                                                            >  2:296112/1‑147 (MQ=255)
                                                                             cgGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAg                                                                         >  2:252657/1‑149 (MQ=255)
                                                                               gCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGc                                                                       <  4:103174/149‑1 (MQ=255)
                                                                                     aTGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCaa                                                                 <  5:131402/149‑1 (MQ=255)
                                                                                                cAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcgc                                                      <  7:137428/149‑1 (MQ=255)
                                                                                                cAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCcgc                                                      <  4:91666/149‑1 (MQ=255)
                                                                                                  gTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGa                                                    <  3:354513/149‑1 (MQ=255)
                                                                                                      aaaaGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTca                                                <  8:51693/149‑1 (MQ=255)
                                                                                                       gcaGCGCCGTCGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTc                                                                                                         <  8:25700/89‑1 (MQ=255)
                                                                                                       gcaGCGCCGTCGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTc                                                                                                         >  7:25700/3‑91 (MQ=255)
                                                                                                            gcCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGaa                                          <  2:267874/149‑1 (MQ=255)
                                                                                                             ccGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAAc                                         <  5:3722/149‑1 (MQ=255)
                                                                                                             ccGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAAc                                         <  7:321266/149‑1 (MQ=255)
                                                                                                                  ggCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCAt                                    <  3:217536/149‑1 (MQ=255)
                                                                                                                     gcTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGa                                 <  7:280073/149‑1 (MQ=255)
                                                                                                                             cTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATc                         <  1:296112/149‑1 (MQ=255)
                                                                                                                               gggCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTgg                                                                                       >  4:135276/1‑85 (MQ=255)
                                                                                                                               gggCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTgg                                                                                       <  3:135276/85‑1 (MQ=255)
                                                                                                                                    gTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCca                  >  1:324406/1‑149 (MQ=255)
                                                                                                                                     ccatggTCGAACACATTGTGCTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCcac                 >  8:250671/6‑149 (MQ=255)
                                                                                                                                      cAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTAATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCcaca                >  7:128175/1‑149 (MQ=255)
                                                                                                                                        aTGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACAcc              <  4:61873/149‑1 (MQ=255)
                                                                                                                                            cGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCCGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTc          >  8:231813/1‑149 (MQ=255)
                                                                                                                                             gAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCcacgc               <  2:30923/143‑3 (MQ=255)
                                                                                                                                             gAGCACAATGTGTTCGACCATGGGGTCAGTCATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCcacgc               >  1:30923/1‑141 (MQ=255)
                                                                                                                                                  caATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCacac               >  6:130101/1‑138 (MQ=255)
                                                                                                                                                  caATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCacac               <  5:130101/138‑1 (MQ=255)
                                                                                                                                                   aaTGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGa   >  5:222483/1‑149 (MQ=255)
                                                                                                                                                    aTGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGAt  <  8:335776/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAAGCAGACGCAATCAGGTAACCCCCTGAGTACGCCATGAAGTTAACCAGGCCGGTGATTGGTTTGATCTTGGTCGCGGCCCGAATGTCCGCCGCCAGTTCAAAAGCGCCGACGGCGCTGCCGCCTGGGCTGTCAATGTCGAGCACAATGTGTTCGACCATGGGGTCAGTGATGGCCTTATTCAGGGCTGCTCGCAGCCCTTCATAGCTGGTCATGGTCTCGCAGGCGTTCAAGTGGCTGCCGCGACTCACAAGAACGCCATGGACGGGGATCACCTCCACACCCGTCTGGGCGAT  >  NC_002947/4399374‑4399670

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