Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,399,267 G→A 5.1% P223S (CCG→TCG)  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,2670GA 5.1% 237.2 / 2.3 79P223S (CCG→TCG) PP_3878minor capsid protein C
Reads supporting (aligned to +/- strand):  ref base G (35/40);  new base A (2/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

ATGGCGGAGGAGCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAG  >  NC_002947/4399119‑4399416
                                                                                                                                                     |                                                                                                                                                     
aTGGCGGAGGAGCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGc                                                                                                                                                      <  2:35428/150‑1 (MQ=255)
aTGGCGGAGGAGCCGCGATAACAAGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGGCGTGAACAGCTGATAGCTTTCTGGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGc                                                                                                                                                      >  2:1685662/1‑150 (MQ=255)
  ggCGGAGGAGCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTc                                                                                                                                                    <  2:1303946/150‑1 (MQ=255)
         gagCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCggg                                                                                                                                             >  1:1050098/1‑150 (MQ=255)
         gagCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCggg                                                                                                                                             >  2:519876/1‑150 (MQ=255)
             cgcgATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCa                                                                                                                                         <  2:1404776/150‑1 (MQ=255)
             cgcgATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCAGCTGGAGTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCa                                                                                                                                         <  1:1310524/150‑1 (MQ=255)
                   aaCAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCAt                                                                                                                                   >  1:1820962/1‑150 (MQ=255)
                      aGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCgg                                                                                                                                              >  2:2315761/1‑136 (MQ=255)
                      aGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCgg                                                                                                                                              <  1:2315761/136‑1 (MQ=255)
                       ggCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCAttttt                                                                                                                               <  2:1326700/150‑1 (MQ=255)
                          aGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTtgtg                                                                                                                            <  2:342809/150‑1 (MQ=255)
                            ccTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTtgtgtg                                                                                                                          >  1:608876/1‑150 (MQ=255)
                                 gTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCg                                                                                                                     >  2:2029819/1‑150 (MQ=255)
                                      gatgatGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCaaa                                                                                                                >  2:674079/1‑150 (MQ=255)
                                      gatgatGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCaaa                                                                                                                <  1:240979/150‑1 (MQ=255)
                                        tgatgTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGa                                                                                                              <  1:1065500/150‑1 (MQ=255)
                                          atgTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGAcc                                                                                                            >  1:1534837/1‑150 (MQ=255)
                                             tCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGTTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGtt                                                                                                         <  1:519876/150‑1 (MQ=255)
                                             tCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAGGACCGtg                                                                                                         >  2:1334021/1‑149 (MQ=255)
                                             tCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGtt                                                                                                         <  2:608876/150‑1 (MQ=255)
                                               cGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTtgt                                                                                                       <  1:674079/150‑1 (MQ=255)
                                               cGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACGTGCAGAGATTGCCAGGCATGCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTtgt                                                                                                       <  1:1334021/150‑1 (MQ=255)
                                                 cAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGa                                                                                                     <  2:1534837/150‑1 (MQ=255)
                                                  aaCATCACGGCCGCGTTATTCTGAAACTGCGGTCGTCAACCGCTATTAGCTTTCTTTCACAACTTCCTTAGTCACTTGCAGATATTTCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGGCCGTTGTGAc                                                                                                    <  2:722636/150‑1 (MQ=255)
                                                     aTCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTTGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCtt                                                                                                 <  2:903896/150‑1 (MQ=255)
                                                             cgcgGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCa                                                                                         >  1:1978940/1‑150 (MQ=255)
                                                              gcgGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCAt                                                                                        <  1:764368/150‑1 (MQ=255)
                                                              gcgGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGCTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCAt                                                                                        <  1:407314/150‑1 (MQ=255)
                                                                ggTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGc                                                                                      <  1:2306912/150‑1 (MQ=255)
                                                                ggTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGc                                                                                      >  2:269887/1‑150 (MQ=255)
                                                                    ttCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTc                                                                                  <  1:140124/150‑1 (MQ=255)
                                                                     tCGGCAACGTGGGTCGTGAACAGCTGATCGCTTTCATCTTGGTGTTCTTGTGCGACCTGCAGAGATTGCTCCGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCg                                                                                 <  1:1685662/150‑1 (MQ=255)
                                                                     tCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCg                                                                                 <  2:1820962/150‑1 (MQ=255)
                                                                        gCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGatc                                                                              <  2:1978940/150‑1 (MQ=255)
                                                                             gTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGa                                                                                                              >  1:1153077/1‑113 (MQ=255)
                                                                             gTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGa                                                                                                              <  2:1153077/113‑1 (MQ=255)
                                                                              tGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTc                                                                        >  2:1405912/1‑150 (MQ=255)
                                                                                 gTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCATACCCATGCCTTCGATCATCTTCGAg                                                                     <  2:1481863/150‑1 (MQ=255)
                                                                                  tCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGc                                                                    >  2:2026666/1‑150 (MQ=255)
                                                                                     tGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAgg                                                                                                                                       <  1:1242879/80‑1 (MQ=255)
                                                                                     tGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAgg                                                                                                                                       >  2:1242879/1‑80 (MQ=255)
                                                                                        aCAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCa                                                              >  1:2083662/1‑150 (MQ=255)
                                                                                        aCAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCa                                                              >  1:2299217/1‑150 (MQ=255)
                                                                                          aGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTAGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCAtg                                                            <  2:388546/150‑1 (MQ=255)
                                                                                                    tttCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTgc                                                                                                                         >  1:2004767/1‑79 (MQ=255)
                                                                                                    tttCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTgc                                                                                                                         <  2:2004767/79‑1 (MQ=255)
                                                                                                        ttGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCcaccac                                              >  1:1357446/1‑150 (MQ=255)
                                                                                                         tGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACcaca                                             <  2:2299217/150‑1 (MQ=255)
                                                                                                            aCGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGCTCATCTTCGAGCGATCCATGTGACTGGCCACCACACCg                                          >  1:839146/1‑150 (MQ=255)
                                                                                                            aCGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCg                                          >  2:2087333/1‑150 (MQ=255)
                                                                                                            aCGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCg                                          >  2:636846/1‑150 (MQ=255)
                                                                                                             cGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGa                                         >  1:393431/1‑150 (MQ=255)
                                                                                                              gACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTg                                                      <  1:2235236/136‑1 (MQ=255)
                                                                                                              gACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTg                                                      >  2:2235236/1‑136 (MQ=255)
                                                                                                              gACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGAt                                        <  2:2083662/150‑1 (MQ=255)
                                                                                                              gACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGACTCGTTCGGGCTCAGGTCAttttt                                                                                                                               >  1:514921/1‑63 (MQ=255)
                                                                                                              gACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGACTCGTTCGGGCTCAGGTCAttttt                                                                                                                               <  2:514921/63‑1 (MQ=255)
                                                                                                                cTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGtt                                                                                                         >  1:2169198/1‑83 (MQ=255)
                                                                                                                cTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGtt                                                                                                         <  2:2169198/83‑1 (MQ=255)
                                                                                                                 ttCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGAGATCCATGTGACTGGCCACCACACCGATGGa                                     >  2:759781/1‑150 (MQ=255)
                                                                                                                   cGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGAcc                                   <  2:1050098/150‑1 (MQ=255)
                                                                                                                   cGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGAcc                                   >  2:1255279/1‑150 (MQ=255)
                                                                                                                   cGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGAcc                                   <  1:1405912/150‑1 (MQ=255)
                                                                                                                              ttGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGa                         <  1:2029819/150‑1 (MQ=255)
                                                                                                                                gCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCa                                                                                         <  2:522022/84‑1 (MQ=255)
                                                                                                                                gCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCa                                                                                         >  1:522022/1‑84 (MQ=255)
                                                                                                                                gCAGAGATTGCTCGGTAATCGACTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGAcc                                                                                                            >  2:1326694/1‑65 (MQ=255)
                                                                                                                                gCAGAGATTGCTCGGTAATCGACTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGAcc                                                                                                            <  1:1326694/65‑1 (MQ=255)
                                                                                                                                  agagATTGCTCGGTAATCGGCTCGTTCGGGCTCCGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTc                     >  1:1337025/1‑150 (MQ=255)
                                                                                                                                  agagATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCAGTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTc                     >  1:931107/1‑150 (MQ=255)
                                                                                                                                  agagATTGCTCGGTAATAGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTc                     >  2:877083/1‑150 (MQ=255)
                                                                                                                                   gagaTTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTTGCCACCACACCGATGGACCCTACACCGGAGGTCa                    <  1:1353947/150‑1 (MQ=255)
                                                                                                                                   gagaTTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCa                    >  2:1335279/1‑150 (MQ=255)
                                                                                                                                     gaTTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAgg                  <  2:1337025/150‑1 (MQ=255)
                                                                                                                                          cTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTgc                                                                                                                         <  1:1721129/42‑1 (MQ=255)
                                                                                                                                          cTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTgc                                                                                                                         >  2:1721129/1‑42 (MQ=255)
                                                                                                                                           tCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAAcc            >  1:1208045/1‑150 (MQ=255)
                                                                                                                                                aaTCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGAGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGAc       <  1:759781/150‑1 (MQ=255)
                                                                                                                                                     gCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAg  >  2:620491/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGGCGGAGGAGCCGCGATAACAGGCAGCCTCGGTGGCGATGATGTCCGCAACATCACGGCCGCGGTATTCGGCAACGTGGGTCGTGAACAGCTGATAGCTTTCTTGCACGACTTCGTTGAGCACTTGCAGAGATTGCTCGGTAATCGGCTCGTTCGGGCTCAGGTCATTTTTGTGTGCGCCGGCAAAGACCGTTGTGACCTTCACACCCATGCCTTCGATCATCTTCGAGCGATCCATGTGACTGGCCACCACACCGATGGACCCTACACCGGAGGTCAGGCTAACCACGACTTCAG  >  NC_002947/4399119‑4399416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: