Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 327,999 G→A S245S (TCC→TCT)  betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913327,9990GA96.6% 1176.9 / 18.0 351S245S (TCC→TCT) betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base G (8/4);  new base A (176/163);  total (184/167)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 359 total aligned reads displayed.

GAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGG  >  NC_000913/327851‑328146
                                                                                                                                                    |                                                                                                                                                   
gAATGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAgaga                                                                                                                                                   <  3:1342619/151‑1 (MQ=255)
    gAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAagag                                                                                                                                               <  3:277645/151‑1 (MQ=255)
     aaGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAgg                                                                                                                                              <  5:303499/151‑1 (MQ=255)
       gACGCGGCTGCCATTTGTACACTCCTGAGCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAACACGATCAGCGGTGATTTAGCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGgc                                                                                                                                             <  8:833150/150‑1 (MQ=255)
           cgGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGccgccg                                                                                                                                        <  7:1299643/151‑1 (MQ=255)
            gggTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGa                                                                                                                                       <  7:1582278/151‑1 (MQ=255)
             ggTGCCATTGGTACACACCTGACCGGTGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCGAGTTCCATGGTCACTTCTTTCAGAGAAGAGgccggggag                                                                                                                                      <  4:1303148/151‑6 (MQ=255)
             ggTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGccgccg                                                                                                                                        >  8:1531548/1‑149 (MQ=255)
               ggCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGtt                                                                                                                                    <  3:333051/150‑1 (MQ=255)
                ggcATTGGTAGACAGCTGACGGGAGCTGAAGTAGTTTGCGATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTTCCGCCCAGTTCCATGGTCACTTCTTTGAGAGAAGAGGCCGCCGAGtt                                                                                                                                    <  2:570711/148‑1 (MQ=255)
                gCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGtt                                                                                                                                    <  6:994021/150‑1 (MQ=255)
                  cATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGc                                                                                                                                 >  3:164312/1‑151 (MQ=255)
                  cATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGc                                                                                                                                 >  3:793047/1‑151 (MQ=255)
                      ggTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCAt                                                                                                                              <  2:792182/150‑1 (MQ=255)
                       gTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCa                                                                                                                            <  2:19591/151‑1 (MQ=255)
                           acacCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACttt                                                                                                                        >  3:94984/1‑151 (MQ=255)
                              ccTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACtttttt                                                                                                                     <  3:648415/151‑1 (MQ=255)
                                tGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGc                                                                                                                   >  4:275665/1‑151 (MQ=255)
                                tGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGc                                                                                                                   >  2:751471/1‑151 (MQ=255)
                                   ccGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGc                                                                                                                >  5:1137296/1‑151 (MQ=255)
                                   ccGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGc                                                                                                                <  3:489208/151‑1 (MQ=255)
                                     ggAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTg                                                                                                              >  6:489108/1‑151 (MQ=255)
                                      gAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTgg                                                                                                             >  8:1017957/1‑151 (MQ=255)
                                       aGCTGAAGCAGTTTGCCCTCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGGGATTTACCGCCCAGTTCCATGGTCACTGCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGc                                                                                                            <  2:224103/151‑1 (MQ=255)
                                           gaagaaGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCGATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGac                                                                                                         >  8:708391/1‑150 (MQ=255)
                                              gaacgtggCCATCATGGCGATCTCGGCGGCGAGATCGAGATCCGCATCAGCGAAAACGATCAGCGCTGATGTCCCGCCCAGTTCCAGGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccg                                                                                                     <  6:1248857/144‑1 (MQ=255)
                                               aagTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCCCTTTTTTGCCGCTGGCGACAccgc                                                                                                    >  8:624185/1‑151 (MQ=255)
                                               aagTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAACGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgc                                                                                                    <  8:1308775/151‑1 (MQ=255)
                                                agTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCGGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgcc                                                                                                   >  2:137496/1‑151 (MQ=255)
                                                agTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgcc                                                                                                   <  1:1140613/151‑1 (MQ=255)
                                                 gTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACAccgccg                                                                                                  <  6:904519/151‑1 (MQ=255)
                                                   ttGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCgg                                                                                                 <  5:1601890/150‑1 (MQ=255)
                                                   ttGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCgg                                                                                                 <  3:153276/150‑1 (MQ=255)
                                                       catcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAAt                                                                                            <  5:122847/151‑1 (MQ=255)
                                                       catcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAAt                                                                                            >  1:1018180/1‑151 (MQ=255)
                                                       catcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAAt                                                                                            >  2:1533725/1‑151 (MQ=255)
                                                        atcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATg                                                                                           >  5:773676/1‑151 (MQ=255)
                                                         tcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGCCACCGCCGGTAAATg                                                                                           >  8:753124/1‑150 (MQ=255)
                                                         tcatGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATg                                                                                           >  5:850039/1‑150 (MQ=255)
                                                              gCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACtt                                                                                     <  4:1117888/151‑1 (MQ=255)
                                                               cGATATCGGCGGCGAGATCGAGATCCGCATCATCGCAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACttt                                                                                    >  7:1332068/1‑151 (MQ=255)
                                                               cGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACttt                                                                                    >  7:186618/1‑151 (MQ=255)
                                                                    tCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCa                                                                                >  2:260203/1‑150 (MQ=255)
                                                                     gggcggcGAGATGGAGATCCGCATGATCGTAAAGGATGTGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                              <  7:523409/150‑1 (MQ=255)
                                                                     cggcggCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAAt                                                                              <  7:404843/151‑1 (MQ=255)
                                                                      ggcggcGAGATCGAGATCCGCATCATCGACAACGAGCCGCGGGGATTTACCGCCCACTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTTGCGACACCGCCGGTAAATGACACTTTGGCAATg                                                                             <  1:644958/151‑1 (MQ=255)
                                                                        cggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGcc                                                                           >  3:1249111/1‑151 (MQ=255)
                                                                         ggcgAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGccc                                                                          >  4:1017010/1‑151 (MQ=255)
                                                                          gcgAGCTCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGACCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCAGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCg                                                                         <  7:331688/151‑1 (MQ=255)
                                                                           cgagatcgagatCCGCATCATCGAAAACGATCAGCGCTGATTTACCGCCCAGTTCCATCGACACTTCTTTAAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCAGGCGACACCGCCCGTAAATGACACTTTGCCAATCCCCgg                                                                        >  4:421109/1‑151 (MQ=255)
                                                                             agatcgagatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGAt                                                                      >  4:149445/1‑151 (MQ=255)
                                                                              gatcgagatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTGCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATg                                                                     <  8:133412/151‑1 (MQ=255)
                                                                               atcgagatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGc                                                                    <  6:408688/151‑1 (MQ=255)
                                                                                   agatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCgg                                                                <  1:993242/151‑1 (MQ=255)
                                                                                    gatcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGt                                                               <  5:994301/151‑1 (MQ=255)
                                                                                     atcCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTc                                                              <  5:924625/151‑1 (MQ=255)
                                                                                       ccGCATCATCGAAAACGATCAGCGGTCATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCAGTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAg                                                            <  6:1521788/151‑1 (MQ=255)
                                                                                        cGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTc                                                              >  1:389580/1‑148 (MQ=255)
                                                                                             catcGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATAtt                                                       <  3:1302594/150‑1 (MQ=255)
                                                                                             catcGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTg                                                      >  1:739287/1‑151 (MQ=255)
                                                                                             catcCAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATAtt                                                       <  7:688363/150‑1 (MQ=255)
                                                                                              atcGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGc                                                     <  1:1582508/151‑1 (MQ=255)
                                                                                                cGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGccc                                                   >  5:382122/1‑151 (MQ=255)
                                                                                                   aaaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGt                                                <  6:696985/151‑1 (MQ=255)
                                                                                                    aaCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtc                                               >  7:883389/1‑151 (MQ=255)
                                                                                                     aCGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGtct                                              <  5:138322/151‑1 (MQ=255)
                                                                                                          cAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcgcg                                         <  7:1072816/151‑1 (MQ=255)
                                                                                                          cAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCcgc                                          <  5:1395915/150‑1 (MQ=255)
                                                                                                             cGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGccc                                      >  3:858005/1‑151 (MQ=255)
                                                                                                             cGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATCCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGccc                                      >  6:1480189/1‑151 (MQ=255)
                                                                                                              ggTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCa                                     <  7:81148/151‑1 (MQ=255)
                                                                                                              ggTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCa                                     <  8:487186/151‑1 (MQ=255)
                                                                                                               gTGATTTACCGCCCAGTTCGATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTATTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCAc                                    <  8:184717/151‑1 (MQ=255)
                                                                                                                  aTTTACCGCCCAGTTCCATGGTCACTTCTTTCAGATAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGcc                                 <  5:737232/151‑1 (MQ=255)
                                                                                                                   tttACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGccc                                >  4:1552516/1‑151 (MQ=255)
                                                                                                                    ttACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGccc                                <  2:1662681/150‑1 (MQ=255)
                                                                                                                    ttACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCg                               >  3:1137985/1‑151 (MQ=255)
                                                                                                                     tACCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCgg                              >  2:345743/1‑151 (MQ=255)
                                                                                                                      aCCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCgg                              <  4:904902/150‑1 (MQ=255)
                                                                                                                      aCCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGc                             >  4:797907/1‑151 (MQ=255)
                                                                                                                      aCCGCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGACTCCGCGCCCACGCCCGGc                             >  6:1116147/1‑151 (MQ=255)
                                                                                                                         gCCCAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCaa                           >  8:1128361/1‑150 (MQ=255)
                                                                                                                          cccAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGCCACTTTGGCAATGACCGGATGCTCGGTCAGATATTGCACGGTCTCCGCGCCCACCCCCGCCAacc                         >  5:1017713/1‑150 (MQ=255)
                                                                                                                           ccAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAacac                        >  2:663187/1‑151 (MQ=255)
                                                                                                                           ccAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAacac                        >  4:526576/1‑151 (MQ=255)
                                                                                                                            cAGTTCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACg                       >  7:612208/1‑151 (MQ=255)
                                                                                                                               ttCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGtt                     >  8:540387/1‑150 (MQ=255)
                                                                                                                                tCCATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTa                    <  8:1061505/150‑1 (MQ=255)
                                                                                                                                 ccATGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTaa                   >  1:1458851/1‑150 (MQ=255)
                                                                                                                                   aTGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAAt                 >  3:174440/1‑150 (MQ=255)
                                                                                                                                   aTGGTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATa                <  6:835909/151‑1 (MQ=255)
                                                                                                                                      gTCACTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGc             <  2:1280489/151‑1 (MQ=255)
                                                                                                                                         aCTTCTTTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGt          >  6:298055/1‑151 (MQ=255)
                                                                                                                                           ttcttTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTcc        >  1:323495/1‑151 (MQ=255)
                                                                                                                                           ttcttTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTcc        >  5:1490456/1‑151 (MQ=255)
                                                                                                                                            tcttTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCg       >  2:452671/1‑151 (MQ=255)
                                                                                                                                             cttTCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGCTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCgg      >  7:56828/1‑151 (MQ=255)
                                                                                                                                                tCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGCCACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGc     >  2:1395450/1‑149 (MQ=255)
                                                                                                                                                tCAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGc     >  5:197538/1‑149 (MQ=255)
                                                                                                                                                 cAGAGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAgg  >  3:1042350/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAACGAAGACGCGGGTGCCATTGGTACACACCTGACCGGAGCTGAAGAAGTTTGCCATCATGGCGATATCGGCGGCGAGATCGAGATCCGCATCATCGAAAACGATCAGCGGTGATTTACCGCCCAGTTCCATGGTCACTTCTTTCAGGGAAGAGGCCGCCGAGTTAGCCATCACTTTTTTGCCGCTGGCGACACCGCCGGTAAATGACACTTTGGCAATGCCCGGATGCTCGGTCAGATATTGCCCGGTCTCCGCGCCCACGCCCGGCAACACGTTAAATACGCCGTCCGGCAGG  >  NC_000913/327851‑328146

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