Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,883 T→A 5.5% D332E (GAT→GAA 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,8830TA 5.5% 231.7 / 3.6 73D332E (GAT→GAAlonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

CGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCG  >  NC_000913/459735‑460030
                                                                                                                                                    |                                                                                                                                                   
cGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGAt                                                                                                                                                     <  1:274421/149‑1 (MQ=255)
cGATGTCGGCAGAAGCGACCGGAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGAt                                                                                                                                                     <  2:323588/149‑1 (MQ=255)
 gacgTCGGCAGAAGCGACCGTAGTGCGTGGCTATATCGACAGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTAGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                                                                                    <  6:320916/146‑1 (MQ=255)
  aTGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                                                                                    >  2:209311/1‑148 (MQ=255)
        gCAGAAGCGTCCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCt                                                                                                                                             <  6:64241/149‑1 (MQ=255)
         cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                            <  4:138298/149‑1 (MQ=255)
         cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                            >  8:259973/1‑149 (MQ=255)
         cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                            >  3:209305/1‑149 (MQ=255)
         cAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCtt                                                                                                                                            <  7:308031/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                                                                                          >  2:240642/1‑149 (MQ=255)
                                                           tGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCct                                                                                          >  3:235780/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                                                                                   <  8:154264/149‑1 (MQ=255)
                                                                  cgcgTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTg                                                                                   <  7:259973/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                                                                <  6:71109/149‑1 (MQ=255)
                                                                                     aGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtag                                                                <  1:253266/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:249357/120‑1 (MQ=255)
                                                                                                                                            gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                      >  8:255534/1‑120 (MQ=255)
                                                                                                                                            gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                      >  8:249357/1‑120 (MQ=255)
                                                                                                                                            gTGAAAGAACGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCAtt                                      <  7:255534/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)
                                                                                                                                                    |                                                                                                                                                   
CGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCG  >  NC_000913/459735‑460030

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: