Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,429 T→G 27.9% 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,4290TG27.9% 89.7 / 7.0 63intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (13/31);  new base G (17/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.60e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACT  >  NC_000913/422287‑422576
                                                                                                                                              |                                                                                                                                                   
gTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtccgcctc                                                                                                                                              >  2:21470/1‑150 (MQ=255)
   gATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCAGCCCCCGGCCTCCTcccc                                                                                                                                           >  2:469373/1‑150 (MQ=255)
       gttgtGCTGTTGATTGGCTGGATGTTTTAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGGAAGTTATCCGCCCCCGGGCGCCGCCCGCaaa                                                                                                                                       >  2:23250/1‑150 (MQ=255)
          gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATACGCCCCCGTACTGCTCCCCCAAatat                                                                                                                                    >  2:824990/1‑150 (MQ=255)
             cTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCATCGTCCCCCAACTATCCt                                                                                                                                 >  2:1117035/1‑150 (MQ=255)
                 tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGACGCCGTCCTCCTCCCCCAAAGATAATTCat                                                                                                                             >  1:189036/1‑149 (MQ=255)
                   aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTAGTCCGAAAGTTATCCGCCGGCGGCCTCCGGCCCCAAATAAGATTCAgat                                                                                                                           >  2:219173/1‑150 (MQ=255)
                    ttGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCTCCCGCGGCCTCAAGACCCAAAAAAACTTCAgaat                                                                                                                          >  1:202469/1‑148 (MQ=255)
                      ggCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCGGTCGGCCGCCCCCAAAGATACTTCAgatgtt                                                                                                                        >  2:1085182/1‑148 (MQ=255)
                      ggCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                                                                                                        <  1:1010025/150‑1 (MQ=255)
                      ggCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                                                                                                        <  1:730845/150‑1 (MQ=255)
                         tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCGACCTACCCCTCATGTACTTAAGGAGATGAg                                                                                                                     >  1:341601/1‑150 (MQ=255)
                         tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAg                                                                                                                     <  1:410355/150‑1 (MQ=255)
                             tGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                 <  1:522443/150‑1 (MQ=255)
                               tttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATcc                                                                                                               <  1:1053492/150‑1 (MQ=255)
                               tgtAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCATGATCCGGGTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATcc                                                                                                               <  1:828668/148‑1 (MQ=255)
                                 tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTACTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTg                                                                                                             <  1:134359/150‑1 (MQ=255)
                                       ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                                                                                                       <  1:779171/150‑1 (MQ=255)
                                         gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCGTCCTCCCCCAAATATCATTAAGAAGATGAGTGTTCCTGGAGTAggg                                                                                                     >  2:733275/1‑149 (MQ=255)
                                         gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                     <  2:546122/150‑1 (MQ=255)
                                         gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCACAAACACCAAAAATCCTTAAGATGATGAGTGATCCTGCAATAGGc                                                                                                     >  1:655173/1‑150 (MQ=255)
                                          ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCt                                                                                                    <  1:309418/150‑1 (MQ=255)
                                           cACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCATGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCACCGGCCTGCTCCCCCAAATAACCATCAGAACATGATGGTTCCTGCATTAGGCtt                                                                                                   >  2:1346443/1‑149 (MQ=255)
                                             cGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                 <  1:353501/150‑1 (MQ=255)
                                              gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATgg                                                                                                <  2:43368/150‑1 (MQ=255)
                                               aTCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                               <  2:731084/150‑1 (MQ=255)
                                                tCGTCAGCTTGCTGAAAACCAGTAATACTTGCGTTCTGTAAGCCGGTTCAGGAGTAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCGCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCa                                                                                              <  1:671811/150‑1 (MQ=255)
                                                   tCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTGAGGAGGAAATCCGGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATTAGTGATCCTGCATTAGGCTATGGCAAgg                                                                                           <  1:101655/150‑1 (MQ=255)
                                                    cAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt                                                                                          <  1:22596/150‑1 (MQ=255)
                                                     aGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTg                                                                                         <  2:1234112/150‑1 (MQ=255)
                                                      gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGa                                                                                        <  2:362344/150‑1 (MQ=255)
                                                           ctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGAGCGTGCATTAGGGTATGGCCAGGTGGTCAgt                                                                                   >  2:306100/1‑149 (MQ=255)
                                                             gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGCCGGCCGCCGCCCCCAAATATCCTACAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                 >  2:3479/1‑150 (MQ=255)
                                                                        aTCCTTGCTTTCTGGAAGCCGGTTCAGGAGGAAAGCCTGCTCCGGATATTCCGAAAGTAATCCGCCCCCGTCCTCCTCCCCCAAATATCATTCAGATGATGAGGGATCCAGCAATAGGCTATGGCAAGGTGATCAGATTTTCATCACAgg                                                                      <  2:204208/150‑1 (MQ=255)
                                                                          ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGTTGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAgggg                                                                    <  1:977189/150‑1 (MQ=255)
                                                                           cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGGCGCCGCCCCCACATAGCCTTCAGCTGAGGAGTGAGCCGGGAGTTGGCGATGGCCAGGTGATCAGCTTTTCATAACTCGGGa                                                                   >  1:679656/1‑150 (MQ=255)
                                                                            ttGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATACGGCTATTCCGAAAGTTATCCGCAACCGTCCTCCGCCCCCAAAAATAGTTGAGATGATGTGAGATCGTGCATGAGTCTAAGGGTACCAGCTGAGAATATCATCACAGGGGaa                                                                  >  2:135366/1‑150 (MQ=255)
                                                                              gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAtt                                                                <  1:1152517/150‑1 (MQ=255)
                                                                               cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGTAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTa                                                               <  1:469373/150‑1 (MQ=255)
                                                                               cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCTCCCCCGAAGAGCCGGCAGATGATGCGTGATCGGGGTTTAGGCGATGGCTAGGGGATCTGCTGTTCATCACAGGGGAtttt                                                               >  2:1118174/1‑149 (MQ=255)
                                                                                    tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAACGCTCCAGCAGAAGATGAGTGATCCGGCATTAGGGTATGGCAAGGTGATCAGATTTTCAAGACATGGGAGTTatgatt                                                          >  2:429799/1‑149 (MQ=255)
                                                                                     gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCAATGAGATGTTGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAGTATGATGt                                                         >  2:773433/1‑150 (MQ=255)
                                                                                      tAATCCGTTTCAGGTTGAAATCCTGATCCGGATATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATTTCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                        <  1:167539/150‑1 (MQ=255)
                                                                                           cGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATg                                                   <  2:280116/150‑1 (MQ=255)
                                                                                           cGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAAGAACCTGCAGATGATGAGTGATCCTGCATTATGCTATGGCAAGGTGATGAGGTTTTGATAACATGGGAATTATGATTTTAAATg                                                   >  2:1287160/1‑150 (MQ=255)
                                                                                               tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATg                                               <  1:117984/150‑1 (MQ=255)
                                                                                                  ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCCAATATCATTAAGATGATGAGTGCTAGGGGAATAGGCTTTGGGAAGGTGATCAGATTTTGATCACATGGGAATTTTTATTTTTAATTCAtctaa                                            >  2:848178/1‑146 (MQ=255)
                                                                                                  ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCa                                            <  1:1206697/150‑1 (MQ=255)
                                                                                                   gaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCAc                                           <  2:1134647/150‑1 (MQ=255)
                                                                                                              gATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCATCAGATGATGAGTGATCCTGCATTAGGCTATGGTAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcc                                <  1:1027985/150‑1 (MQ=255)
                                                                                                                 ccGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCAATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccatg                             <  1:787952/150‑4 (MQ=255)
                                                                                                                 ccGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCGGCAGGTGccatg                             <  2:216462/150‑4 (MQ=255)
                                                                                                                  cGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCGCCCAAAACTCCGTCCGATGAGGAGGGGTGCAGCGGGAGGGTATGGAAGGGGGAGCAGATTTTCATACCGGGGGTAATATGATGGTTCATGGATGGCCCCTGCCGGGGcccccc                            >  2:59799/1‑149 (MQ=255)
                                                                                                                    gCTATTCCGAAAGTTATCCGCCGCCGGCGGCCGCCCCCCAAAATCCTGAAGATGATGAGTGATCATGCATTAGGAAATGGCTAGGTGATAAGAGTTTAATATGAGGGGAATTATGATGTTAAATGAATGGCACGTGCCCGTTccccccat                          >  2:1001642/1‑147 (MQ=255)
                                                                                                                     cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                         <  2:854452/150‑1 (MQ=255)
                                                                                                                          ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                    <  2:427700/150‑1 (MQ=255)
                                                                                                                          ccGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCAATATCCATCAGATGATGAGTGATCGTGCATTAGGGTATGGGAAGTTGATGAGATTTTATTAAAAGGGGAATTTTGATGTTAAATGAATGGCACATGCCGGTGCCCCCATTTGTTat                    >  2:409560/1‑149 (MQ=255)
                                                                                                                            gAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTAAGCTGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGGTCAGATTGTCATAAACGGGGAATTATGATGTTAAATTAATGGCACCAGCCGGTGCCCCCCATTGTTaaaa                  >  2:685718/1‑149 (MQ=255)
                                                                                                                             aaaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATAAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAca                 <  1:11549/150‑1 (MQ=255)
                                                                                                                                gTTATCCGCCCCCGTCCTCCTCCCCCAAAAATACTTCAGATGATTAGTGATACTGCATTAGGCTATGGCAAGGTGATCAGATTTTAATAACAGGGGAATTATGATGTTAAATGCATGGAACCTGCCGGTGCCCCCATTTGTTAaaaaaaa              >  1:310057/1‑149 (MQ=255)
                                                                                                                                gTTATCCGCCCCCGTCCTACTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaac              >  1:886122/1‑149 (MQ=255)
                                                                                                                                  tATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCATGCAATAGGCTATGGCAAGGTGATGAGAATTTCATCACAGTGGAATTATGATGTTAAATTAATGGCACCTGCCGGTGCCCCCATTTGTTAaaaaaaaaa            >  2:295210/1‑147 (MQ=255)
                                                                                                                                  tATCCGCCCCGTCCTCCTCCCCCAAATATCCTTGAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGGGCCCGCATTTGATAAAcaaagcaa           >  2:213182/1‑150 (MQ=255)
                                                                                                                                     ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAGCAGATTTTCATCAAAGGGGAATAATGATGTTAAATGCATGTAACCTGCCGGGGCCCCCAATTGTTAAAcaaagaaaac         >  2:874142/1‑149 (MQ=255)
                                                                                                                                            cttcctcctcaGCCATATATCCTTCAGATGATGATTGATCCTGCTTTAGGCTTTGGCAAGTTGATCAGATTTTCATCACACGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAGAATGCAAAGATCAACt  <  2:38018/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACT  >  NC_000913/422287‑422576

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: