Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,429 T→G 23.8% intergenic (+70/‑89) proY → / → malZ proline‑specific permease/maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4290TG23.8% 91.7 / 3.9 62intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (13/30);  new base G (14/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.01e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAC  >  NC_000913/422295‑422575
                                                                                                                                      |                                                                                                                                                  
ttgtGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATTCGCCCCCGTCCTCCTCCCCCaaa                                                                                                                                      <  4:46224/149‑1 (MQ=255)
   tgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGTCGGCCGCCCCCaaagaa                                                                                                                                   >  8:262075/1‑146 (MQ=255)
     cTGTTGATTGGCTGGTTTTTTTAACGCCGGCACGATGGTGAGGTGGCTGAAAACGAGTTATCGTTGCGTTCTGTAAGCCGGTTCTGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATcc                                                                                                                                 <  8:399702/149‑1 (MQ=255)
       gTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCTCCTCCCCCCAATATcgt                                                                                                                                >  4:138646/1‑146 (MQ=255)
        ttGATGGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTc                                                                                                                              <  7:189746/149‑1 (MQ=255)
                cTGGATGTTTAAACGCCGCCACGATCGTCAACTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGatgatg                                                                                                                      <  1:38678/149‑1 (MQ=255)
                   gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtc                                                                                                                                                   >  6:295537/1‑117 (MQ=255)
                   gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtc                                                                                                                                                   <  5:295537/117‑1 (MQ=255)
                   gATGTTTAAACGCCGCCACGATCGTAAGCTTGCTTAAAACCAGTAATCCTAGCGTTCTGTAAGCCGGTTCAGGAGGAAATACTGATACGGAAATTCCGAAAGTTATCCGCCCCCGAGCTCCTAACCCATAAATCAATAAGATGATGAGt                                                                                                                   >  4:222425/1‑149 (MQ=255)
                       tttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc                                                                                                               <  4:100267/149‑1 (MQ=255)
                               ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGTTCCTGCATTa                                                                                                       <  2:181558/149‑1 (MQ=255)
                                 gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                     >  3:128045/1‑149 (MQ=255)
                                 gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                     <  8:400284/149‑1 (MQ=255)
                                  ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCGCCTCCTCCCCCAAATAACCTAAAGATGATGAGTTATCATGAATTAgga                                                                                                    >  6:362972/1‑148 (MQ=255)
                                     cGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGTCCGCCTCCCCCCAAAATCCTTCCGATGATGAGTGATCCTGCATGAGGCTAt                                                                                                 >  1:138649/1‑149 (MQ=255)
                                      gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCAAAAACCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                >  2:195691/1‑149 (MQ=255)
                                                  gctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATc                                                                                     <  7:262075/148‑1 (MQ=255)
                                                  gctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCATAAGATGATGAGTGATACTGCATTAGGCAATGGCAAGGTGAt                                                                                      >  6:309790/1‑147 (MQ=255)
                                                   ctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAg                                                                                   <  1:403007/149‑1 (MQ=255)
                                                     gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGTCCTGCGCCCCCAAATATACTTCAGATGATGTGTGATGCTGCATTAGGCTAGGGGAAGGTGGTCAGAt                                                                                 >  1:129635/1‑149 (MQ=255)
                                                     gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg                                                                                          <  5:146347/140‑1 (MQ=255)
                                                     gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCAAAAATCCTTAAGATGATGAGTGATCCGGCATTAGGCTATGGCAAGGttgatcgg                                                                                  >  6:146347/1‑141 (MQ=255)
                                                       aaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                <  8:83414/148‑1 (MQ=255)
                                                        aaCCAGTAATCCGTGCGTTCTGTAAGCCGGGTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtttt                                                                              <  2:171371/149‑1 (MQ=255)
                                                        aaCCAGTAATCATTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATACTTCAGAGTATGAGTGATCCTGCTTTAGGATTTGTCTAGGTTATCAGAtttt                                                                              >  8:389249/1‑149 (MQ=255)
                                                          ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATACTTCAGATGATGAGTGATCATGCATTAGGATATGGCAAGGTGATCAGATTTtct                                                                            >  8:8865/1‑148 (MQ=255)
                                                                aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCACAAATCCTGCAGATGATGAGGGTGCAGGGCTTAGGGTATGGGAAGGTGATCAGATTTTCAGCagag                                                                      >  2:133013/1‑149 (MQ=255)
                                                                 tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAgg                                                                     <  4:128045/149‑1 (MQ=255)
                                                                 tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTACCCAAAATATAAGTCAGAGTATGAGTGCTACGGAAATGGAGTATGGCAAGGGGAGCAGATTTTCATAACAgg                                                                     >  6:14070/1‑149 (MQ=255)
                                                                 tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCACCGCCACCCCCCAATATCCGACCTATGATCAGTGAACCTGCATTAGGCTATGGCAATGTGATCAGATTTTCATcacacg                                                                     >  6:18908/1‑147 (MQ=255)
                                                                  ccTTGCTTTCTGTAAGCCGGTTCAGGAGTAAATCCTGATCCGGCTATTGCGAAAGTTTTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATTATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATGAGATTTTCATCACAgg                                                                     <  4:62410/148‑1 (MQ=255)
                                                                   cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCAGCCCCCCTCCTCATGCCCCAAATATACTTAAGATGATCAGTTATACTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATGACAggc                                                                    >  8:10232/1‑147 (MQ=255)
                                                                     tGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                <  7:204793/134‑1 (MQ=255)
                                                                       cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAtt                                                               <  2:246160/149‑1 (MQ=255)
                                                                       cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAt                                                                <  1:32093/148‑1 (MQ=255)
                                                                        gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCTCCCCCCCATATCCTACAGATGATGAGTGATCCTGGAGTAGGGGATGGCAAGGTGGTCAGATTTTGATCAAAGGGGAATTa                                                              >  4:34367/1‑149 (MQ=255)
                                                                         ttCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTat                                                             >  2:8777/1‑149 (MQ=255)
                                                                            tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCGGCCCCCAACACTCGTGCAGAGGAAGAGTGATAGGGAATAAGGCTATGGCAAGGTGATCAGATTGTCAAGACACGGGCATtattat                                                          >  6:378688/1‑149 (MQ=255)
                                                                              tAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGt                                                        <  1:195691/149‑1 (MQ=255)
                                                                                     gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCACCTCCCCCAAATAAACTTCAGATGATGAGTGATGCTGCATTAGGCTATGGGAAGGGGATCAGATTTTCATGACAGGGGAATTATGATGTTAAAtt                                                  >  8:160898/1‑147 (MQ=255)
                                                                                       tCAGGTTGAAATCCGGATCCGGCTAATGCGTAATTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                               <  2:92889/149‑1 (MQ=255)
                                                                                        gagcaggaAATCCTGATCCGGGAAATCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                               <  2:87519/144‑1 (MQ=255)
                                                                                          ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGTAATTATGATGTTAAATGCATg                                              <  3:239433/147‑1 (MQ=255)
                                                                                          ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAAGATGCTTAAGATGATGAGTGATCGTGAATTAGGATATTGAAAGGGGATCTGTTTATGATCTCATGGGAATTTTGATGTTAATTGAATgtg                                            >  4:239433/1‑147 (MQ=255)
                                                                                            aggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCAc                                          <  6:299928/149‑1 (MQ=255)
                                                                                             ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCAAAAATCCGTAAGTTGATGAGTGATCGTGGATTAGGGGATGGGAAGGTGATAAGATTTTCGTCACAGGGGAATTATGATGTTAAATTGATGGCaa                                          >  6:328765/1‑147 (MQ=255)
                                                                                                 aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCTGGCCTCCTCCCCCAAATATCCTTCAGCTGCTGAGGGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                     <  2:287460/149‑1 (MQ=255)
                                                                                                 aTCCTGATCCGGCTATTCCGAAAGTTATCCCCCCCCCACCTCCTCCCCCAAATATACTTAAAATGAATAGTGCACCAGCATTAGGCTATCTCAAACGGATCAGATTTTCATCACAGGGGAAATATTATTTTAAATGAATGGCAcatacc                                     >  2:160094/1‑144 (MQ=255)
                                                                                                  tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCTCCTCATCCCCCAAATATACTTCAGATGATGAGTGATCATGCATTAGGATATGGCAAGGTGATCATATTTTCATGACAGGGGAATTATGATGTTTAATGCATTTCACATGcag                                    >  2:120908/1‑147 (MQ=255)
                                                                                                    cTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                  <  1:133013/149‑1 (MQ=255)
                                                                                                      gATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCGGCATTAGGCTATGGGAAGGTGTTGAGTTTTTGATCACAGGGGAGTTATGATTTTAAATGAATGGAGCATGCCGGTgg                                >  4:349619/1‑148 (MQ=255)
                                                                                                          cGGCTATTCCGAAAGTTATCCGCCCCCGTCCGACGCACCCAAATATCCTTAAGAAGATGAGTGATACTGCATTAGGCTTAGGAAAGGTGCGCAGATTATCATAAAAAGGGAATTATGATAATTAAGTCATGGCCCCTACCGGAGcgccc                            >  4:312467/1‑149 (MQ=255)
                                                                                                                ttCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATAACCTTCAGATGATGAGTGATACTGCATTAGGGTATGGCAAGGTGATAAGATTTTAATAACAGGGGAATTATTATGTTATATACATGGCACAGGCCGGTGCCCGCAGTTGt                      >  6:234970/1‑149 (MQ=255)
                                                                                                                   cGAAAGTTATCCGCGCCCGGCCGCCGCCCCCAAATATCCTTCAGAAGATGAGTGATGCTGCATTTGGCTATGGCAAGGTGATCAGATTTTAATCTAAGGGGAATTATTATGTTTAATGCATGGCGGATGCAGGGGCCCCCAGTTGTTa                    >  8:221936/1‑148 (MQ=255)
                                                                                                                    gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGc                                                 <  1:66724/118‑1 (MQ=255)
                                                                                                                        gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGTGTGATCCTGCATTAGGCTATGGCATGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa              <  5:18908/149‑1 (MQ=255)
                                                                                                                        gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa              <  3:34367/149‑1 (MQ=255)
                                                                                                                        gTTATCCGCCCCCGTCCTCCTCACCCAAATATCATTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACATGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa              <  5:14070/149‑1 (MQ=255)
                                                                                                                         ttATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaag             <  5:362972/149‑1 (MQ=255)
                                                                                                                           aTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                     <  5:353507/139‑1 (MQ=255)
                                                                                                                               gcacccGTCCTCCTCACCCAAATATCATTCAGATGATGAGTGATCCTTCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa       <  3:222425/146‑1 (MQ=255)
                                                                                                                               gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaag        >  4:371744/1‑148 (MQ=255)
                                                                                                                                    cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc  <  7:53361/149‑1 (MQ=255)
                                                                                                                                    cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc  >  1:1352/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
TTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAC  >  NC_000913/422295‑422575

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