Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,886 A→C 5.7% R333R (CGA→CGC lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,8860AC 5.7% 204.8 / 5.7 71R333R (CGA→CGClonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

GCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGG  >  NC_000913/459743‑460034
                                                                                                                                               |                                                                                                                                                    
gCAGAAGCGTCCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCt                                                                                                                                                 <  6:64241/149‑1 (MQ=255)
 cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                <  7:308031/149‑1 (MQ=255)
 cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                >  8:259973/1‑149 (MQ=255)
 cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                >  3:209305/1‑149 (MQ=255)
 cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                <  4:138298/149‑1 (MQ=255)
   gAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                                                                                              <  6:24736/149‑1 (MQ=255)
      gCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTa                                                                                                                                           >  3:109590/1‑149 (MQ=255)
        gACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATc                                                                                                                                         >  3:291678/1‑149 (MQ=255)
                          tatCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCgtg                                                                                                                       <  2:394554/149‑1 (MQ=255)
                            tCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCgtgt                                                                                                                      >  7:232267/1‑148 (MQ=255)
                            tCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTAAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTc                                                                                                                     <  1:279687/149‑1 (MQ=255)
                                   gATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACaaaa                                                                                                              <  4:209305/149‑1 (MQ=255)
                                       gTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCaa                                                                                                          <  2:278003/149‑1 (MQ=255)
                                       gTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCaa                                                                                                          <  6:302634/149‑1 (MQ=255)
                                        tACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAg                                                                                                         >  5:71109/1‑149 (MQ=255)
                                         aCAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAgg                                                                                                        >  3:159022/1‑149 (MQ=255)
                                            ggTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGAc                                                                                                     <  8:327482/149‑1 (MQ=255)
                                             gTGCCGTGGAATGCGCGTAGCGAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGAcc                                                                                                    >  7:315362/1‑149 (MQ=255)
                                              tGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGAcc                                                                                                    >  7:140505/1‑148 (MQ=255)
                                               gCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                <  6:210667/103‑1 (MQ=255)
                                               gCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                                >  5:210667/1‑103 (MQ=255)
                                                 cGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATc                                                                                                <  7:349007/149‑1 (MQ=255)
                                                   tGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCct                                                                                              <  1:209311/149‑1 (MQ=255)
                                                   tGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCct                                                                                              >  3:235780/1‑149 (MQ=255)
                                                   tGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCct                                                                                              >  2:240642/1‑149 (MQ=255)
                                                     gAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCctct                                                                                            >  3:19338/1‑149 (MQ=255)
                                                     gAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCctc                                                                                             <  3:297510/148‑1 (MQ=255)
                                                      aaTGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTg                                                                                           >  5:47817/1‑149 (MQ=255)
                                                          cgcgTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTg                                                                                       <  7:259973/149‑1 (MQ=255)
                                                          cgcgTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTg                                                                                       <  8:154264/149‑1 (MQ=255)
                                                              tAGCAAGGGCAAAAACGACCTGCGTCAGGCGCAGGAAATCCTTGATAACGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAg                                                                                   <  6:17223/149‑1 (MQ=255)
                                                                tcAATGTCAAATAAGAGCTGCGTCTGGTTCATGAAATCGTTGATACCGAGCATTATTGTGTGGAGCGCGTGAAAGATAGATTCCTTGAGTATCTTGCGGTTCAAAGCAGTGTCAACAAAATCAAGGGACCGATCCTCTGACTGGTAgt                                                                                  <  8:40781/147‑2 (MQ=255)
                                                                gCAAGGTCAAAATAGAGCTGCGTCAGGCGCAGGAAATCCTTGTTACCGACCATTATGGTCTGGTGCGCTTGAAAGATGGAATCCTTGAGTATGTTGCGGTTCAAAGCCGTGTGAACAAAATCAAGGGACCGATCCTCTGCCTGGTAgg                                                                                  <  4:291678/148‑1 (MQ=255)
                                                                gCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAgg                                                                                  <  1:240642/148‑1 (MQ=255)
                                                                gCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATACTTGATACCGACCATTATGGTCTGGAGCTCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAgg                                                                                  <  6:16448/148‑1 (MQ=255)
                                                                     gTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgccgcc                                                                            >  7:373352/1‑149 (MQ=255)
                                                                        aaaaaaGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCggg                                                                         >  3:370669/1‑149 (MQ=255)
                                                                          aaaaGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCggggg                                                                       >  8:294735/1‑149 (MQ=255)
                                                                             aGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtag                                                                    <  1:253266/149‑1 (MQ=255)
                                                                             aGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtag                                                                    <  6:71109/149‑1 (MQ=255)
                                                                              gACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATAGAATCCTTGAGTATCTTGCGGTTCAACGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGGAGGGCCGCCGGGggtagg                                                                   >  1:322928/1‑149 (MQ=255)
                                                                                 cTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtaggta                                                                 >  3:123769/1‑148 (MQ=255)
                                                                                      tCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACct                                                           <  4:123769/149‑1 (MQ=255)
                                                                                        aGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACctct                                                         <  7:294735/149‑1 (MQ=255)
                                                                                           cgcAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCtt                                                      <  8:315362/149‑1 (MQ=255)
                                                                                                gAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCa                                                 <  8:232267/149‑1 (MQ=255)
                                                                                                   aTCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGGAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTc                                              >  7:315638/1‑149 (MQ=255)
                                                                                                    tCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGATCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGCGGTAGGTAAAACCTCTCTTGGTCAGTc                                              >  3:217787/1‑148 (MQ=255)
                                                                                                     ccTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgg                                                                          >  6:300354/1‑119 (MQ=255)
                                                                                                     ccTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgg                                                                          <  5:300354/119‑1 (MQ=255)
                                                                                                     ccTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATATTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTCCATGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCa                                            >  8:159078/1‑149 (MQ=255)
                                                                                                      cTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAt                                           >  1:311193/1‑149 (MQ=255)
                                                                                                       ttGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTTGTCAGTCCAtt                                          >  1:190471/1‑149 (MQ=255)
                                                                                                       ttGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                          >  5:227121/1‑149 (MQ=255)
                                                                                                          aTACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTTTGTAAAACCTCTCTTGGTCAGTCCATTga                                        >  4:65462/1‑147 (MQ=255)
                                                                                                             ccGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCaaa                                    <  5:238360/149‑1 (MQ=255)
                                                                                                              cGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgcc                                                                               <  8:292849/105‑1 (MQ=255)
                                                                                                                    ttATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGTTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAACACCTCTCTTGGTCAGTCAATTGCCAAAGCCACCg                             >  4:110236/1‑149 (MQ=255)
                                                                                                                      aTGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCgg                            <  6:227121/148‑1 (MQ=255)
                                                                                                                       tGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGAGCGTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAACACGTCTCTTGGTCGGTGCATTGCCAAAGCCCCCGGGc                          >  6:385637/1‑149 (MQ=255)
                                                                                                                          tCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGCGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGt                        >  4:364288/1‑148 (MQ=255)
                                                                                                                              gAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGGGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAat                    <  6:47817/148‑1 (MQ=255)
                                                                                                                                  gcgTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGAGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGTGTAGGTAAAACCTCTCTTGGTCTGTCCATTGCCAAAGCCACCGGGCGTAAATATGTc               <  6:312448/149‑1 (MQ=255)
                                                                                                                                   cgTGAAAGATCGAATCCTTGAGTATCTTGCGGGTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTcc              <  5:97267/149‑1 (MQ=255)
                                                                                                                                    gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                          <  7:255534/120‑1 (MQ=255)
                                                                                                                                    gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                          >  8:249357/1‑120 (MQ=255)
                                                                                                                                    gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                          >  8:255534/1‑120 (MQ=255)
                                                                                                                                    gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                          <  7:249357/120‑1 (MQ=255)
                                                                                                                                        aaGATCGAATCCTTGAGTATCTTGCGGTTCAAATCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATg         >  4:73284/1‑149 (MQ=255)
                                                                                                                                         aGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCAATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATgg        >  3:78431/1‑149 (MQ=255)
                                                                                                                                           aTCGAATCCTTGAGTATCTTGCGGTTAAAAGCCTTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGgcg      >  6:145497/1‑149 (MQ=255)
                                                                                                                                               aaTCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTgg  <  4:109590/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGG  >  NC_000913/459743‑460034

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