Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 161,938 C→G 5.2% E32Q (GAG→CAG)  ligT ← 2'‑5' RNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913161,9380CG 5.2% 346.3 / 7.9 97E32Q (GAG→CAG) ligT2'‑5' RNA ligase
Reads supporting (aligned to +/- strand):  ref base C (48/44);  new base G (2/3);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATG  >  NC_000913/161794‑162084
                                                                                                                                                |                                                                                                                                                  
cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                <  2:237450/149‑1 (MQ=255)
cGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGt                                                                                                                                                >  4:707974/1‑149 (MQ=255)
  gCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTgg                                                                                                                                              >  3:382963/1‑149 (MQ=255)
   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                             <  1:504526/149‑1 (MQ=255)
   cGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTggg                                                                                                                                             >  3:477417/1‑149 (MQ=255)
     gcatcaAGGGTGACGGTGAATCAAGGTGGAAGAACTCGCCCGGCTAAAAGAGAAAGCCCCCTGCTTTTATTGGCGGTCGCTGAGCTTAAAAATGCCCGAGTCAGATGCACATTCTCGGCGAAGACCCGAAGTTCCGCTTCAGGTGTGaa                                                                                                                                           <  2:354410/148‑1 (MQ=255)
      catcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGacg                                                                                                                                          <  1:318046/149‑3 (MQ=255)
       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                         >  3:510265/1‑149 (MQ=255)
       atcaAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgt                                                                                                                                         >  3:567591/1‑149 (MQ=255)
         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                       >  3:743656/1‑149 (MQ=255)
         caAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgt                                                                                                                                       <  1:433150/149‑1 (MQ=255)
           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                     <  1:507602/149‑1 (MQ=255)
           aGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAgtgtgt                                                                                                                                     <  1:787522/149‑1 (MQ=255)
            gCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGtgtgtg                                                                                                                                    >  1:184385/1‑149 (MQ=255)
               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                            >  3:7153/1‑138 (MQ=255)
               tGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCACCTGAGGCGGGa                                                                                                                                            <  4:7153/138‑1 (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                                                                                                      >  4:420590/1‑149 (MQ=255)
                                          gtccgGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGAtt                                                                                                      >  2:249428/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                                                                                          >  2:344129/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                          >  4:787750/1‑149 (MQ=255)
                                                      gagaAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAg                                                                                          >  3:36169/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                                                                                                >  1:301016/1‑130 (MQ=255)
                                                                   ctctTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCa                                                                                                <  2:301016/130‑1 (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)
                                                                      tttCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACa                                                                          >  2:667762/1‑149 (MQ=255)
                                                                        tctcGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGa                                                                        <  2:59428/149‑1 (MQ=255)
                                                                            ggCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGtt                                                                    <  1:344129/149‑1 (MQ=255)
                                                                                cTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTc                                                                                                             >  4:9959/1‑104 (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                                                                                         <  4:567858/103‑1 (MQ=255)
                                                                                                     agagTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGt                                                                                         >  3:567858/1‑103 (MQ=255)
                                                                                                         tCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTa                                       <  2:641883/149‑1 (MQ=255)
                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                      <  1:25982/149‑1 (MQ=255)
                                                                                                          cAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTAt                                      <  1:268437/149‑1 (MQ=255)
                                                                                                             aTGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGaa                                   >  1:84981/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                  >  4:51093/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                  >  4:427685/1‑149 (MQ=255)
                                                                                                              tGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                  <  3:420590/149‑1 (MQ=255)
                                                                                                              tGCAAATTATCGACGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAAt                                  <  1:249767/149‑1 (MQ=255)
                                                                                                                 aaaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAc                               >  2:670714/1‑149 (MQ=255)
                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                              >  3:309850/1‑149 (MQ=255)
                                                                                                                  aaTTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTAcc                              >  3:88748/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                                          >  1:624150/1‑116 (MQ=255)
                                                                                                                                       gTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCgg                                          <  2:624150/116‑1 (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)
                                                                                                                                                |                                                                                                                                                  
CGGCGTCATCAAGCGTGAGTGTGAAACCAGGTTGACGAATCCGTCCGGCTAAAAGAGAAAGCGCCTTCTCTTTCTCGGCGCTCACTTCGCCTAAAAATGCCAGAGTCAGATGCAAATTATCGGCGGCGACCGGACGTCCCGCCTCAGGTGGGAAGTGTGTGGCGCGCCAATGGATAATCTGTTCGCGGATTTCTGCAGGTAAGTCGATAGCAAAGAACAGACGTTGCGGTTCAGACATGTGAGGCACTCGGTTATGAATTACCGCGATGCTACAATGTGGCGCGAAGAATG  >  NC_000913/161794‑162084

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: ‑