Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 161,942 T→C 5.9% P30P (CCA→CCG ligT ← 2'‑5' RNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913161,9420TC 5.9% 300.2 / 6.4 99P30P (CCA→CCGligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base C (3/3);  total (51/48)
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.45e-01

Only 100 of 102 total aligned reads displayed.

CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACC  >  NC_000913/161794‑162090
                                                                                                                                                    |                                                                                                                                                    
cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                      >  4:707974/1‑149 (MQ=255)
cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                      <  2:237450/149‑1 (MQ=255)
  gCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTgg                                                                                                                                                    >  3:382963/1‑149 (MQ=255)
   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                                   >  3:477417/1‑149 (MQ=255)
   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                                   <  1:504526/149‑1 (MQ=255)
     gcatcaAGGGTGACGGTGAATCAAGGTGGAAGAACTCGCCCGGCTAAAAGAGAAAGCCCCCTGCTTTTATTGGCGGTCGCTGAGCTTAAAAATGCCCGAGTCAGATGCACATTCTCGGCGAAGACCCGAAGTTCCGCTTCAGGTGTGaa                                                                                                                                                 <  2:354410/148‑1 (MQ=255)
      catcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGacg                                                                                                                                                <  1:318046/149‑3 (MQ=255)
       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                               >  3:567591/1‑149 (MQ=255)
       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                               >  3:510265/1‑149 (MQ=255)
         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                             >  3:743656/1‑149 (MQ=255)
         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                             <  1:433150/149‑1 (MQ=255)
           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                           <  1:787522/149‑1 (MQ=255)
           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                           <  1:507602/149‑1 (MQ=255)
            gCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGtgtgtg                                                                                                                                          >  1:184385/1‑149 (MQ=255)
               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                                  <  4:7153/138‑1 (MQ=255)
               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                                  >  3:7153/1‑138 (MQ=255)
                   tgtgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGcc                                                                                                                                   <  4:510265/149‑1 (MQ=255)
                     tgAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCaa                                                                                                                                 >  1:59428/1‑149 (MQ=255)
                      gAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                >  4:333889/1‑149 (MQ=255)
                        aaCCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGacgt                                                                                                                                               >  4:247396/1‑129 (MQ=255)
                        aaCCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGACgt                                                                                                                                               <  3:247396/132‑1 (MQ=255)
                          ccAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGAt                                                                                                                            >  1:494851/1‑149 (MQ=255)
                          ccAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGAt                                                                                                                            <  3:560214/149‑1 (MQ=255)
                           cAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATa                                                                                                                           <  3:707974/149‑1 (MQ=255)
                            aGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATaa                                                                                                                          <  4:567591/149‑1 (MQ=255)
                             ggTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAAt                                                                                                                         <  1:625231/149‑1 (MQ=255)
                                 gACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCTCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGt                                                                                                                     >  1:654257/1‑149 (MQ=255)
                                   cGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTc                                                                                                                   <  3:68742/149‑1 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                >  3:246425/1‑129 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAAt                                                                                                                                <  4:246425/129‑1 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  2:249428/1‑149 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  4:420590/1‑149 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGACGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                            >  2:249767/1‑149 (MQ=255)
                                             cgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCt                                                                                                         <  1:310504/149‑1 (MQ=255)
                                                cTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      <  2:157086/149‑1 (MQ=255)
                                                    aaGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTa                                                                                                  >  2:268437/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  3:36169/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  4:787750/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  2:344129/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGCGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                                >  4:544490/1‑149 (MQ=255)
                                                           aGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGAt                                                                                           <  1:590494/149‑1 (MQ=255)
                                                            gcgcCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATa                                                                                          >  3:527919/1‑149 (MQ=255)
                                                                   ctctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      <  2:301016/130‑1 (MQ=255)
                                                                   ctctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                      >  1:301016/1‑130 (MQ=255)
                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATgg                                                                                                                              >  3:696000/1‑105 (MQ=255)
                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATgg                                                                                                                              <  4:696000/105‑1 (MQ=255)
                                                                    tctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCAGAGACCGCACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCAATTCGCGTATTTATGCATGTAAGTCATTAGCAAAGaa                                                                                  <  2:584481/149‑1 (MQ=255)
                                                                      tttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACa                                                                                <  3:787750/149‑1 (MQ=255)
                                                                        tctcGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGa                                                                              <  2:59428/149‑1 (MQ=255)
                                                                            ggCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGtt                                                                          <  1:344129/149‑1 (MQ=255)
                                                                                cTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTc                                                                                                                   <  3:9959/104‑1 (MQ=255)
                                                                                cTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCgg                                                                      <  1:684640/149‑1 (MQ=255)
                                                                                   aCTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTc                                                                   >  1:376532/1‑149 (MQ=255)
                                                                                     ttCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAg                                                                 >  1:608972/1‑149 (MQ=255)
                                                                                      tCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGa                                                                <  3:333889/149‑1 (MQ=255)
                                                                                         ccTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAt                                                             >  4:621999/1‑149 (MQ=255)
                                                                                          cTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAtg                                                            >  1:545857/1‑149 (MQ=255)
                                                                                           tAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGAGGGAAGTGCGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACAtgt                                                           >  3:201665/1‑149 (MQ=255)
                                                                                              aaaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        <  4:743656/149‑1 (MQ=255)
                                                                                               aaTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAgg                                                       >  2:222614/1‑149 (MQ=255)
                                                                                                aTGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGc                                                      >  3:681433/1‑149 (MQ=255)
                                                                                                   ccAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGa                                                                                            <  3:50989/108‑1 (MQ=255)
                                                                                                   ccAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGa                                                                                            >  4:50989/1‑108 (MQ=255)
                                                                                                     agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGt                                                                                               >  3:567858/1‑103 (MQ=255)
                                                                                                     agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGt                                                                                               <  4:567858/103‑1 (MQ=255)
                                                                                                         tCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTa                                             <  2:641883/149‑1 (MQ=255)
                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                            <  1:268437/149‑1 (MQ=255)
                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                            <  1:25982/149‑1 (MQ=255)
                                                                                                             aTGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGaa                                         >  1:84981/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        >  4:427685/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        <  3:420590/149‑1 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        >  4:51093/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGACGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                        <  1:249767/149‑1 (MQ=255)
                                                                                                                 aaaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAc                                     >  2:670714/1‑149 (MQ=255)
                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                                    >  3:88748/1‑149 (MQ=255)
                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                                    >  3:309850/1‑149 (MQ=255)
                                                                                                                    ttATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                  >  3:684914/1‑149 (MQ=255)
                                                                                                                    ttATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACcgc                                  <  2:376532/149‑1 (MQ=255)
                                                                                                                       tCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGAt                               <  4:201665/149‑1 (MQ=255)
                                                                                                                                  cGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                    <  1:670714/149‑1 (MQ=255)
                                                                                                                                  cGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAAtgtg                    <  3:621999/149‑1 (MQ=255)
                                                                                                                                   ggACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTTGGTTATGAATTACCGCGATGCTACAATGTgg                   <  2:545857/149‑1 (MQ=255)
                                                                                                                                    gACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGgc                  >  1:453296/1‑149 (MQ=255)
                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                <  2:624150/116‑1 (MQ=255)
                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                                >  1:624150/1‑116 (MQ=255)
                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGcgcg               >  2:209556/1‑149 (MQ=255)
                                                                                                                                         cccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaa             >  4:507464/1‑149 (MQ=255)
                                                                                                                                         cccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaa             >  1:498776/1‑149 (MQ=255)
                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        <  2:496373/105‑1 (MQ=255)
                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAg                                                        >  1:496373/1‑105 (MQ=255)
                                                                                                                                          ccGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaag            >  4:32663/1‑149 (MQ=255)
                                                                                                                                            gCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  3:544490/149‑1 (MQ=255)
                                                                                                                                             ccTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAAt         <  4:82966/149‑1 (MQ=255)
                                                                                                                                              cTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg        <  2:216519/149‑1 (MQ=255)
                                                                                                                                              cTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATg        <  1:249428/149‑1 (MQ=255)
                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          >  2:604012/1‑142 (MQ=255)
                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          >  2:310483/1‑142 (MQ=255)
                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  1:604012/142‑1 (MQ=255)
                                                                                                                                                   gTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCgaagaa          <  1:310483/142‑1 (MQ=255)
                                                                                                                                                    tGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAAcc  >  3:79723/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATGTTAACC  >  NC_000913/161794‑162090

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: