Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 991,864 G→A 5.9% A415T (GCA→ACA)  pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913991,8640GA 5.9% 218.6 / 2.6 68A415T (GCA→ACA) pepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base A (2/2);  total (33/35)
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.89e-01

CCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCT  >  NC_000913/991717‑992009
                                                                                                                                                   |                                                                                                                                                 
ccGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTg                                                                                                                                                   >  5:61898/1‑148 (MQ=255)
 cGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCa                                                                                                                                                 <  5:25169/149‑1 (MQ=255)
   gATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTAATGATGGTAGTGCAGc                                                                                                                                               >  6:135105/1‑149 (MQ=255)
       tGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCTTCATGATGGTAGTGCAGCGAcc                                                                                                                                           >  3:67178/1‑149 (MQ=255)
         gTCATTGAGATGAACAACTTCTACACCCTGCCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCt                                                                                                                                          <  2:189555/148‑1 (MQ=255)
            aTTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                                                                                                                                      >  1:77691/1‑149 (MQ=255)
             ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgac                                                                                                                                     >  5:271425/1‑149 (MQ=255)
             ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgac                                                                                                                                     >  7:105808/1‑149 (MQ=255)
                agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgac                                                                                                                                  >  8:183997/1‑149 (MQ=255)
                 gaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACt                                                                                                                                 <  7:39092/149‑1 (MQ=255)
                   tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACttt                                                                                                                               <  4:109043/149‑1 (MQ=255)
                   tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACttt                                                                                                                               >  2:245920/1‑149 (MQ=255)
                   tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCCTCATGATGGTAGTGCAGCGACCTGTGACGACttt                                                                                                                               <  4:67178/149‑1 (MQ=255)
                       caacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtgtg                                                                                                                            >  1:53251/1‑148 (MQ=255)
                        aacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtgtg                                                                                                                            >  2:304763/1‑147 (MQ=255)
                          cTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                                                                                                                        >  6:265598/1‑149 (MQ=255)
                               acacCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATg                                                                                                                   >  8:279447/1‑149 (MQ=255)
                                cacCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGTATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTAGTGCAGGCGATgg                                                                                                                  <  1:324002/149‑1 (MQ=255)
                                cacCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATgg                                                                                                                  <  6:61898/149‑1 (MQ=255)
                                cacCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATg                                                                                                                   <  8:105808/148‑1 (MQ=255)
                                  cccTGACAGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCATGCTTTGCGAAGAAAACTTCCAGAAAGTGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa                                                                                                                <  8:390317/149‑1 (MQ=255)
                                     tGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGAt                                                                                                             >  8:372808/1‑149 (MQ=255)
                                       aCCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGAt                                                                                                             <  2:77691/147‑1 (MQ=255)
                                        ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                                                                                                                           >  2:108532/1‑132 (MQ=255)
                                        ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                                                                                                                           <  1:108532/132‑1 (MQ=255)
                                         cGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGt                                                                                                         <  1:304763/149‑1 (MQ=255)
                                         cGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGt                                                                                                         <  5:265598/149‑1 (MQ=255)
                                                 agaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTc                                                                                                 <  2:53251/149‑1 (MQ=255)
                                                     gggCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGCCCTCTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGAtc                                                                                             >  6:163810/1‑149 (MQ=255)
                                                     gggCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGAtc                                                                                             <  4:392926/149‑1 (MQ=255)
                                                           ggAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTcc                                                                                        >  1:316294/1‑148 (MQ=255)
                                                               gTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCAttt                                                                                   <  7:183997/149‑1 (MQ=255)
                                                                tGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTc                                                                                  >  4:404374/1‑149 (MQ=255)
                                                                  aTTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCAAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTccg                                                                                >  7:390179/1‑149 (MQ=255)
                                                                          gATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCTTTgg                                                                         >  4:188651/1‑148 (MQ=255)
                                                                           aTCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTa                                                                       >  1:414173/1‑149 (MQ=255)
                                                                               acacCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTACAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGc                                                                   >  1:144234/1‑149 (MQ=255)
                                                                               acacCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTACAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGc                                                                   >  1:431918/1‑149 (MQ=255)
                                                                                       cTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTcc                                                             <  7:279447/147‑1 (MQ=255)
                                                                                       cTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCgg                                                           <  2:316294/149‑1 (MQ=255)
                                                                                       cTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTACAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCgg                                                           <  2:431918/149‑1 (MQ=255)
                                                                                       cTTGGCGAAGAAAAATTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTACAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCgg                                                           <  2:144234/149‑1 (MQ=255)
                                                                                           gCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTaca                                                       >  5:19586/1‑149 (MQ=255)
                                                                                              aagaaAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACAcc                                                     >  8:285421/1‑148 (MQ=255)
                                                                                               agaaAACTTCCAGAAAGGGATGCAGCTTTATTTCGCGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCTTTGGTACAGCCAGTCCGGTACACCGa                                                   >  4:244015/1‑149 (MQ=255)
                                                                                                 aaaaCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTTCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAtt                                                 >  8:421816/1‑149 (MQ=255)
                                                                                                 aaaaCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGAtt                                                 <  3:404374/149‑1 (MQ=255)
                                                                                                             aaGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCaaa                                     <  3:97708/149‑1 (MQ=255)
                                                                                                              aGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAg                                    <  6:36893/149‑1 (MQ=255)
                                                                                                               gggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAga                                   <  6:271425/149‑1 (MQ=255)
                                                                                                                ggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgac                                  <  3:188651/149‑1 (MQ=255)
                                                                                                                ggATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGGCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgac                                  <  5:135105/149‑1 (MQ=255)
                                                                                                                    gCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga                                >  1:233945/1‑147 (MQ=255)
                                                                                                                    gCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAgacga                                <  2:233945/147‑1 (MQ=255)
                                                                                                                     cAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTa                             <  2:414173/149‑1 (MQ=255)
                                                                                                                              ttCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCgg                     >  2:365557/1‑148 (MQ=255)
                                                                                                                                  aGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGCTGGAAGATGCTTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAccc                <  7:372808/149‑1 (MQ=255)
                                                                                                                                   gCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTTGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAAcc                >  3:132507/1‑148 (MQ=255)
                                                                                                                                    cGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGa              <  7:285421/149‑1 (MQ=255)
                                                                                                                                        atgatgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCa           <  1:245920/148‑1 (MQ=255)
                                                                                                                                        atgatgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg          >  6:318484/1‑149 (MQ=255)
                                                                                                                                        atgatgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAg          <  5:318484/149‑1 (MQ=255)
                                                                                                                                         tgatgGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGt         <  2:153834/149‑1 (MQ=255)
                                                                                                                                             ggTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTAGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAACCCGAGCTGTacac     >  1:413286/1‑149 (MQ=255)
                                                                                                                                              gTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACAcc    <  6:19586/149‑1 (MQ=255)
                                                                                                                                                aGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCt  >  4:553/1‑149 (MQ=255)
                                                                                                                                                 gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCg                      <  6:29885/128‑1 (MQ=255)
                                                                                                                                                 gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCt  >  1:187105/1‑148 (MQ=255)
                                                                                                                                                 gTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATACg                      >  5:29885/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTCCCATTTCCGCCGTTGGTACAGCCAGTCCGGTACACCGATTGTGACCGTCAAAGACGACTACAATCCGGAAACCGAGCAGTACACCCT  >  NC_000913/991717‑992009

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