Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,410,822 G→C 7.1% H121Q (CAC→CAG PP_3894 ← replication protein O

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,410,8220GC 7.1% 219.1 / 7.5 70H121Q (CAC→CAGPP_3894replication protein O
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (2/3);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

GAGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCT  >  NC_002947/4410674‑4410967
                                                                                                                                                    |                                                                                                                                                 
gAGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgt                                                                                                                                                  >  2:700709/1‑150 (MQ=255)
  ggggTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgt                                                                                                                                                <  1:920013/150‑1 (MQ=255)
   gggTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgtg                                                                                                                                               <  1:363168/150‑1 (MQ=255)
     gTTGCCTTGTCGTTGCCGTGACCGCCTTCCGCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTa                                                                                                                                             >  2:715222/1‑150 (MQ=255)
         ccTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGAtt                                                                                                                                         >  2:1082235/1‑150 (MQ=255)
           ttGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgtg                                                                                                                                               >  2:1688243/1‑142 (MQ=255)
           ttGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtgtg                                                                                                                                               <  1:1688243/142‑1 (MQ=255)
           ttGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGCGTAGATTTc                                                                                                                                       >  2:1492206/1‑150 (MQ=255)
             gTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACACTGTGTAGaaaacac                                                                                                                                     <  2:1321845/150‑7 (MQ=255)
                 ttGCCGTGACCGCCTTCCTCGGCGCCAGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACg                                                                                                                                 <  2:1639000/150‑1 (MQ=255)
                    ccGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGCGTAGATTTCAACACGGtt                                                                                                                              <  1:1492206/150‑1 (MQ=255)
                     cGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTCTGTGTGTCAGTGTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTg                                                                                                                             <  2:2133241/150‑1 (MQ=255)
                         aCCGCCTTCCTCGGCGCCGGTTGGTGAGGAGGCGTCGACAACCAATGCGTTAGCCTTACTAGATTTCCGTCTTTGGTTTAGGCGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCtt                                                                                                                         <  1:655939/150‑1 (MQ=255)
                             ccTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACATAGTGTGTAGATTTCAACACGGTTGGCTTGCTg                                                                                                                     >  1:1453057/1‑150 (MQ=255)
                             ccTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTg                                                                                                                     >  1:2352547/1‑150 (MQ=255)
                              cTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCt                                                                                                                      >  2:15040/1‑148 (MQ=255)
                              cTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCt                                                                                                                      <  1:15040/148‑1 (MQ=255)
                                        gcCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCa                                                                                                          <  1:290858/150‑1 (MQ=255)
                                        gcCGGTCGATGCGGATGCGTCGACAACCAATGCGTCAGCATTGGTATTTGTATTTCTGTCTATTATGTGTGTACTTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTAGCTTGCTGTCGATTCTccc                                                                                                          >  2:180177/1‑149 (MQ=255)
                                            gTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCa                                                                                                      >  2:2088589/1‑150 (MQ=255)
                                               gatgcggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACAGAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGt                                                                                                   >  1:530937/1‑150 (MQ=255)
                                               gatgcggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGt                                                                                                   >  1:988143/1‑150 (MQ=255)
                                               gatgcggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGt                                                                                                   <  1:1082235/150‑1 (MQ=255)
                                                 tgcggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTAc                                                                                                 <  1:2088589/150‑1 (MQ=255)
                                                   cggatgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCg                                                                                               <  2:2128256/150‑1 (MQ=255)
                                                      atgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCt                                                                                                                      <  1:1598854/124‑1 (MQ=255)
                                                      atgcgTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCt                                                                                                                      >  2:1598854/1‑124 (MQ=255)
                                                           tCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGaaa                                                                                       <  2:48655/150‑1 (MQ=255)
                                                             gACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTCATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTGGAGCACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAgg                                                                                     <  2:944093/150‑2 (MQ=255)
                                                                aaCCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGaaa                                                                                       <  1:2096330/145‑1 (MQ=255)
                                                                aaCCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGaaa                                                                                       >  2:2096330/1‑145 (MQ=255)
                                                                 aCCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                                 <  2:280518/150‑1 (MQ=255)
                                                                 aCCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGAtt                                                                                 >  2:1374711/1‑150 (MQ=255)
                                                                   cAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgt                                                                                                                                                  >  2:2279624/1‑83 (MQ=255)
                                                                   cAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgt                                                                                                                                                  <  1:2279624/83‑1 (MQ=255)
                                                                    aaTGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGAttt                                                                                                                                        <  1:1362158/92‑1 (MQ=255)
                                                                    aaTGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGAttt                                                                                                                                        >  2:1362158/1‑92 (MQ=255)
                                                                         gTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtg                                                                                                                                                 >  1:1328080/1‑78 (MQ=255)
                                                                         gTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAgtg                                                                                                                                                 <  2:1328080/78‑1 (MQ=255)
                                                                          aacgCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTg                                                                                                                             <  1:1886979/94‑1 (MQ=255)
                                                                          aacgCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTg                                                                                                                             >  2:1886979/4‑97 (MQ=255)
                                                                                     atttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGt                                                                                                                    >  1:1202544/1‑95 (MQ=255)
                                                                                     atttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGt                                                                                                                    <  2:1202544/95‑1 (MQ=255)
                                                                                       ttgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgc                                                           <  1:700709/150‑1 (MQ=255)
                                                                                        tgtatttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATcgcg                                                          >  2:1622991/1‑150 (MQ=255)
                                                                                             ttCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTACCGCCATCGCGGaaaa                                                     <  1:715222/150‑1 (MQ=255)
                                                                                                 gTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACAAAGTGTGTAGATTTAAACAGGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGGAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAGGCAc                                                 >  2:516504/1‑150 (MQ=255)
                                                                                                  gcTTTTTTGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCAcg                                                >  2:1611041/2‑150 (MQ=255)
                                                                                                        aTGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGa                                          <  1:1374711/150‑1 (MQ=255)
                                                                                                             tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGa                                                                                                                 >  2:2179216/1‑74 (MQ=255)
                                                                                                             tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGa                                                                                                                 <  1:2179216/74‑1 (MQ=255)
                                                                                                             tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGAt                                     >  2:734852/1‑150 (MQ=255)
                                                                                                             tgtAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGAt                                     >  2:68224/1‑150 (MQ=255)
                                                                                                               tAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATCTCGACAGTGGTGGTTTCCTTAAAAATCTCAACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCa                                   >  2:295717/1‑150 (MQ=255)
                                                                                                                aaTTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGATTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCaa                                  <  1:1651106/150‑1 (MQ=255)
                                                                                                                 aTTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGGAAAGCACGCGCTGACGGATCAAc                                 >  1:985103/1‑150 (MQ=255)
                                                                                                                  ttttCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACg                                >  1:2079448/1‑150 (MQ=255)
                                                                                                                    ttCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGaa                              <  2:540681/150‑1 (MQ=255)
                                                                                                                    gtCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGaa                              >  2:1845221/2‑150 (MQ=255)
                                                                                                                        acacAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGa                          <  1:734852/150‑1 (MQ=255)
                                                                                                                          acaGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATc                        <  2:2352547/150‑1 (MQ=255)
                                                                                                                          acaGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGGTGGCTTCCTGTCAATTCTCAACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATc                        >  2:233925/1‑150 (MQ=255)
                                                                                                                            aGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGc                      <  1:1611041/150‑1 (MQ=255)
                                                                                                                                      ttgtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACg            >  1:1552383/2‑150 (MQ=255)
                                                                                                                                      gtgtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACg            >  1:702086/1‑150 (MQ=255)
                                                                                                                                      gggtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACg            >  1:1383251/3‑150 (MQ=255)
                                                                                                                                       tgtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGacac             >  2:889774/1‑148 (MQ=255)
                                                                                                                                       tgtTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGacac             <  1:889774/148‑1 (MQ=255)
                                                                                                                                           atCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTg                                           >  1:562361/2‑114 (MQ=255)
                                                                                                                                           atCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTg                                           <  2:562361/113‑1 (MQ=255)
                                                                                                                                           atCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAAc                                 <  1:1982115/123‑1 (MQ=255)
                                                                                                                                           atCTACACACTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAAc                                 >  2:1982115/2‑124 (MQ=255)
                                                                                                                                            tCTACACAGTGTGTAGATCTCAGCACTGGCGGCTTCCTGTCAATTCTCCACTCAGGTACAGGTCGGAAAGACATTGAATCACGACTGCCGCCATCCCCGAGAAGCACGCGCTGACGTAGCAACAACTTCATCTCGCGAGACAAGTTGgcg      >  2:746167/1‑150 (MQ=255)
                                                                                                                                              tACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGgcgcg    >  1:1272991/1‑150 (MQ=255)
                                                                                                                                                cacaGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCt  >  2:2010973/1‑150 (MQ=255)
                                                                                                                                                 acaGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCg                                                                                               <  1:1756979/56‑1 (MQ=255)
                                                                                                                                                 acaGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCg                                                                                               >  2:1756979/1‑56 (MQ=255)
                                                                                                                                                 acaGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGac               >  1:976232/1‑136 (MQ=255)
                                                                                                                                                 acaGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGac               <  2:976232/136‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAGGGGTTGCCTTGTCGTTGCCGTGACCGCCTTCCTCGGCGCCGGTCGATGCGGATGCGTCGACAACCAATGCGTTAGCATTGGTATTTGTATTTCTGTCTTTTATGTGTGTAATTTTCAACACAGTGGCAGGGGTGTTTTCTACACAGTGTGTAGATTTCAACACGGTTGGCTTGCTGTCGATTCTCCACTCAGGTACCGGTAGGAAAGTGATTGGGTCGCGACTGCCGCCATCGCGGAAAAGCACGCGCTGACGGATCAACGAATTGATCGCCCGAGACACGTTGGCGCGCT  >  NC_002947/4410674‑4410967

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