Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 504,920 T→C 6.2% L3L (TTA→CTA)  ushA → bifunctional UDP‑sugar hydrolase/5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913504,9200TC 6.2% 2533.3 / 101.9 784L3L (TTA→CTA) ushAbifunctional UDP‑sugar hydrolase/5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (382/353);  new base C (22/27);  total (404/380)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-02

Only 100 of 831 total aligned reads displayed.

ATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAATATGG  >  NC_000913/504772‑505068
                                                                                                                                                    |                                                                                                                                                    
  gCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAatta                                                                                                                                                    >  2:278263/1‑149 (MQ=255)
  gCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAatta                                                                                                                                                    <  2:1181848/149‑1 (MQ=255)
   cgtcgtCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattat                                                                                                                                                   <  4:6219751/149‑1 (MQ=255)
    gtcgtcTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattatt                                                                                                                                                  <  4:402422/149‑1 (MQ=255)
     tcgtcTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattatt                                                                                                                                                  >  1:591767/1‑148 (MQ=255)
          tGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCg                                                                                                                                            >  3:582215/1‑149 (MQ=255)
           gaaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCgg                                                                                                                                           <  3:4310770/149‑1 (MQ=255)
            aaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCTCCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTg                                                                                                                                                 >  2:6238632/1‑142 (MQ=255)
            aaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCggg                                                                                                                                          <  3:5997281/149‑1 (MQ=255)
                  gTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGggcgtgg                                                                                                                                    <  4:5957925/149‑1 (MQ=255)
                   tctcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGggcgtggc                                                                                                                                   >  3:3744421/1‑149 (MQ=255)
                     tcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattatt                                                                                                                                                  <  2:2543931/132‑1 (MQ=255)
                     tcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGggcgtggcgt                                                                                                                                 <  3:3203282/149‑1 (MQ=255)
                     tcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGggcgtggcgt                                                                                                                                 <  4:5825167/149‑1 (MQ=255)
                      cAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGtt                                                                                                                                >  2:6825973/1‑149 (MQ=255)
                       aaCGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTa                                                                                                                               >  1:4230638/1‑149 (MQ=255)
                       aaCGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTa                                                                                                                               <  4:3517351/149‑1 (MQ=255)
                        aCGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTGGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAg                                                                                                                              >  2:5939179/1‑149 (MQ=255)
                           ccGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAgcgc                                                                                                                           <  2:7078125/149‑1 (MQ=255)
                            cGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCt                                                                                                                          <  4:3060603/149‑1 (MQ=255)
                             gAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTg                                                                                                                         >  4:4589124/1‑149 (MQ=255)
                             gAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTg                                                                                                                         >  1:5203420/1‑149 (MQ=255)
                               aTACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGtt                                                                                                                       <  1:7017436/149‑1 (MQ=255)
                               aTACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGtt                                                                                                                       <  3:4154873/149‑1 (MQ=255)
                                 aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTaa                                                                                                                     <  4:1248491/149‑1 (MQ=255)
                                 aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTCTTaa                                                                                                                     <  1:6177604/149‑1 (MQ=255)
                                  cAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAactat                                                                                                                                                   <  4:673718/118‑1 (MQ=255)
                                  cAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAAc                                                                                                                    <  4:4830846/149‑1 (MQ=255)
                                      aTTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACcaca                                                                                                                <  1:3833048/149‑1 (MQ=255)
                                                 ggATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGc                                                                                                     >  3:4567801/1‑149 (MQ=255)
                                                  gATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCg                                                                                                    >  1:1554958/1‑149 (MQ=255)
                                                         aTTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             <  4:2107273/149‑1 (MQ=255)
                                                         aTTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             <  2:2024322/149‑1 (MQ=255)
                                                         aTTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCCGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             <  1:4524934/149‑1 (MQ=255)
                                                         aTTTACCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACTCTGGCGAGTGaaa                                                                                             <  1:3138027/149‑1 (MQ=255)
                                                                ttCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGctct                                                                                      >  2:796469/1‑149 (MQ=255)
                                                                ttCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGctct                                                                                      <  3:5238954/149‑1 (MQ=255)
                                                                 tCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTg                                                                                     <  3:6496165/149‑1 (MQ=255)
                                                                   aCCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGc                                                                                   <  2:4871886/149‑1 (MQ=255)
                                                                    ccGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGc                                                                                   >  2:4585054/1‑148 (MQ=255)
                                                                         cGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGa                                                                             <  4:274554/149‑1 (MQ=255)
                                                                            aGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCa                                                                          >  1:225844/1‑149 (MQ=255)
                                                                             gACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGggcgtgg                                                                                                                                    <  2:323402/90‑1 (MQ=255)
                                                                             gACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGggcgtgg                                                                                                                                    >  1:323402/1‑90 (MQ=255)
                                                                             gACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAg                                                                         >  1:180956/1‑149 (MQ=255)
                                                                             gACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAg                                                                         <  2:5174999/149‑1 (MQ=255)
                                                                                ttGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGAt                                                                      >  3:953239/1‑149 (MQ=255)
                                                                                   tCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGa                                                                                                   >  1:1242698/1‑117 (MQ=255)
                                                                                   tCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAAc                                                                                            >  4:15833/1‑124 (MQ=255)
                                                                                    cTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATaaaa                                                                  <  1:1003536/149‑1 (MQ=255)
                                                                                     tATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAAc                                                                 >  3:2454431/1‑149 (MQ=255)
                                                                                     tATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAAc                                                                 >  1:6736129/1‑149 (MQ=255)
                                                                                     tATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAAc                                                                 >  2:2210055/1‑149 (MQ=255)
                                                                                        gatgTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTa                                                              <  3:153219/149‑1 (MQ=255)
                                                                                        gatgTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTa                                                              <  2:1679916/149‑1 (MQ=255)
                                                                                         atgTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTAc                                                             <  3:787632/149‑1 (MQ=255)
                                                                                            tCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             >  2:1385757/1‑114 (MQ=255)
                                                                                             cgcgTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACaaaa                                                         <  4:5505975/149‑1 (MQ=255)
                                                                                              gcgTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAAc                                                                                                                    >  3:2623386/1‑89 (MQ=255)
                                                                                              gcgTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAAc                                                                                                                    <  4:2623386/89‑1 (MQ=255)
                                                                                              gcgTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACACAAt                                                        >  3:3450960/1‑149 (MQ=255)
                                                                                                  cATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACa                                                    <  4:7023248/149‑1 (MQ=255)
                                                                                                    tACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             >  1:560450/1‑106 (MQ=255)
                                                                                                     aCTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGtt                                                 >  3:2175871/1‑149 (MQ=255)
                                                                                                      cTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTAc                                                             <  2:5560540/136‑1 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGTGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCt                                               >  3:1686152/1‑149 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACt                                                                                           <  1:5412188/105‑1 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGa                                                                       >  4:6482879/1‑125 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGGTAAAACCTACAAAATTACAGTTCt                                               >  3:5584984/1‑149 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCt                                               >  3:3937749/1‑149 (MQ=255)
                                                                                                       tATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCt                                               >  3:680144/1‑149 (MQ=255)
                                                                                                         tttttCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGc                                             >  4:5101308/1‑149 (MQ=255)
                                                                                                         tttttCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGc                                             >  1:6330066/1‑149 (MQ=255)
                                                                                                            ttCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATa                                          <  1:2370583/149‑1 (MQ=255)
                                                                                                            ttCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGGGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATa                                          >  2:3761427/1‑149 (MQ=255)
                                                                                                               aaCACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCa                                       >  3:190444/1‑149 (MQ=255)
                                                                                                                 cacGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAAt                                     <  3:5483458/149‑1 (MQ=255)
                                                                                                                  acGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGaaa                                                                                             <  3:1093205/92‑1 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGatc                                 <  3:444287/149‑1 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGatc                                 <  2:4273759/149‑1 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGatc                                 <  4:3744421/149‑1 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGatc                                 <  2:3076772/149‑1 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTa                                                                                >  4:3263012/1‑102 (MQ=255)
                                                                                                                     ttGAAATCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGatc                                 >  1:2091694/1‑149 (MQ=255)
                                                                                                                          aTCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTa                                                                                                            <  1:4998837/69‑1 (MQ=255)
                                                                                                                          aTCAGGTCAGGGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATcatcat                            >  4:1379040/1‑149 (MQ=255)
                                                                                                                           tcaggtcaggGAGAGAAGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATg                           <  2:4975157/149‑1 (MQ=255)
                                                                                                                                  aggGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCAttt                    >  4:4343887/1‑149 (MQ=255)
                                                                                                                                      agagaAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAATCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTgg                <  2:5798/149‑1 (MQ=255)
                                                                                                                                        agaAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGgcg              >  3:3008719/1‑149 (MQ=255)
                                                                                                                                           aGTATGAAACTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCaa           <  1:734492/149‑1 (MQ=255)
                                                                                                                                            gTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAAt          <  4:5271454/149‑1 (MQ=255)
                                                                                                                                            gTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAAt          <  1:2798617/149‑1 (MQ=255)
                                                                                                                                             tATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATg         >  2:3451276/1‑149 (MQ=255)
                                                                                                                                               tGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAAcc                                                                <  3:1327457/92‑1 (MQ=255)
                                                                                                                                                  aaTTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAAc                                                                                            <  3:4970072/61‑1 (MQ=255)
                                                                                                                                                  aaTTATTGCAGCGGGGCGTGGCGTCAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAatat    <  2:6573562/149‑1 (MQ=255)
                                                                                                                                                   attattGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAATATg   <  2:4009058/149‑1 (MQ=255)
                                                                                                                                                   attattGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAATATg   <  1:4874089/149‑1 (MQ=255)
                                                                                                                                                    ttattGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAATATgg  >  1:3573094/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTACACTGGCGAGTGAAACTGCTCTGGCGTATGAGCAGGATAAAACCTACAAAATTACAGTTCTGCATACCAATGATCATCATGGGCATTTTTGGCGCAATGAATATGG  >  NC_000913/504772‑505068

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