Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,008,626 G→A 6.0% V373I (GTT→ATT) ‡ ygeX → 2,3‑diaminopropionate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,008,6260GA 6.0% 259.9 / 3.3 82V373I (GTT→ATT) ‡ygeX2,3‑diaminopropionate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (40/37);  new base A (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAG  >  NC_000913/3008477‑3008772
                                                                                                                                                     |                                                                                                                                                  
tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                                                                                    >  2:949240/1‑150 (MQ=255)
tGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTgg                                                                                                                                                    >  2:105973/1‑150 (MQ=255)
 gCTGGGTAATCCGTACGGCAACGACCCGCGGATCATCTGCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAATGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGt                                                                                                                                                   <  1:490771/150‑1 (MQ=255)
 gCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGt                                                                                                                                                   <  2:908009/150‑1 (MQ=255)
  cTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCACAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGtt                                                                                                                                                  <  1:140969/150‑1 (MQ=255)
     ggTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATc                                                                                                                                               <  1:211420/150‑1 (MQ=255)
        aaTCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGc                                                                                                                                            >  2:1351262/1‑150 (MQ=255)
          tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                                                                                          <  1:1031223/150‑1 (MQ=255)
          tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                                                                                          >  2:1230883/1‑150 (MQ=255)
          tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                                                                                          <  1:105973/150‑1 (MQ=255)
          tCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCAc                                                                                                                                          <  1:10083/150‑1 (MQ=255)
               aCGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCCTGGTGCTGGTTATCACCACTGaac                                                                                                                                     >  2:1297559/1‑149 (MQ=255)
               aCGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGGTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGATAACCAGCAccacgg                                                                                                                                     <  1:790112/150‑6 (MQ=255)
                 ggCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGc                                                                                                                                   <  1:1094564/150‑1 (MQ=255)
                  gCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCg                                                                                                                                  >  1:888517/1‑150 (MQ=255)
                  gCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCg                                                                                                                                  <  2:1203914/150‑1 (MQ=255)
                   cAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGa                                                                                                                                 >  2:1182486/1‑150 (MQ=255)
                    aaCGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGac                                                                                                                                <  2:319452/150‑1 (MQ=255)
                       gACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACAcc                                                                                                                             <  2:618480/150‑1 (MQ=255)
                           cgcgCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACg                                                                                                                         <  1:226560/150‑1 (MQ=255)
                              gcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGa                                                                                                                      <  2:1034993/150‑1 (MQ=255)
                              gcATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGa                                                                                                                      <  2:380806/150‑1 (MQ=255)
                                atcatcTCCGGTGAATCCGGCGCTGTCCGTTTGGGCGTTATCGCAGCGGATCATTAACACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAAAAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAg                                                                                                                    >  2:22312/1‑150 (MQ=255)
                                 tcatcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGc                                                                                                                   <  1:1529036/150‑1 (MQ=255)
                                   atcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGACGCAc                                                                                                                 >  1:1086385/1‑150 (MQ=255)
                                   atcTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCAc                                                                                                                 >  1:1233958/1‑150 (MQ=255)
                                    tctcCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACt                                                                                                                >  2:764742/1‑150 (MQ=255)
                                       ccGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTAcc                                                                                                             >  1:10809/1‑150 (MQ=255)
                                         ggTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTGATTATCAGCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAAGAAAGATGCCGTGGTGCTGGTTTTCAGCACTGAAGTCGACACCGACGTGAAGCACTACggc                                                                                                           <  2:1520332/150‑1 (MQ=255)
                                                  ggCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTc                                                                                                  <  1:433476/150‑1 (MQ=255)
                                                  ggCGCTGTCGGATTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGACGCACTACCGCGAAGTTGTc                                                                                                  >  2:436118/1‑150 (MQ=255)
                                                        gTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAATGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGaa                                                                                            >  1:1115925/1‑150 (MQ=255)
                                                          cgttttGGTCGTTCTCGCAGAGGTTAATTATAACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAgg                                                                                          >  2:750536/4‑150 (MQ=255)
                                                            ttttGGGCGTTCTCGCAGCGGTTCATTATCACCCTCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTTCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCa                                                                                        <  1:660412/149‑1 (MQ=255)
                                                            gTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCTCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCa                                                                                        <  1:656227/150‑1 (MQ=255)
                                                             tttGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCATAACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCaa                                                                                       <  2:962347/150‑1 (MQ=255)
                                                                gggCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGAGAAGTTGTCTGGGAAGGCAAaca                                                                                    >  2:1260304/1‑150 (MQ=255)
                                                                   cGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGc                                                                                 <  2:10809/150‑1 (MQ=255)
                                                                      gctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGt                                                                              <  1:330546/149‑1 (MQ=255)
                                                                       ctcGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTa                                                                             >  2:173898/1‑150 (MQ=255)
                                                                         cGCAGCGGTTCATTATCCCCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCAGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGc                                                                           >  2:863708/1‑150 (MQ=255)
                                                                         cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGc                                                                                                                   >  2:854921/1‑110 (MQ=255)
                                                                         cGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGc                                                                                                                   <  1:854921/110‑1 (MQ=255)
                                                                             gCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGTCAAACACGCTGTAGCACCt                                                                       <  2:888517/150‑1 (MQ=255)
                                                                              cGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCtt                                                                      <  1:1182486/150‑1 (MQ=255)
                                                                                gTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTaa                                                                    >  1:544756/1‑150 (MQ=255)
                                                                                gTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTaa                                                                    >  1:584769/1‑150 (MQ=255)
                                                                                  tCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATc                                                                  <  2:1153126/150‑1 (MQ=255)
                                                                                    aTTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAc                                                                >  1:687156/1‑150 (MQ=255)
                                                                                     ttATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAcc                                                               >  2:557755/1‑150 (MQ=255)
                                                                                     ttATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAcc                                                               >  2:962733/1‑150 (MQ=255)
                                                                                      tATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCAccc                                                              >  2:1153175/1‑150 (MQ=255)
                                                                                         cACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTAt                                                           >  1:1291122/1‑150 (MQ=255)
                                                                                          aaacccAACTTCAAAGCATGATCGAAAAACTGGATATGAACAAATATGCCGTGGTGATGGTTATCAGCAATGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTAtt                                                          >  2:450862/6‑150 (MQ=255)
                                                                                             cGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGa                                                       <  1:1230883/150‑1 (MQ=255)
                                                                                              ccAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGaa                                                      <  2:1045941/149‑1 (MQ=255)
                                                                                               cAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAAc                                                     >  2:1493819/1‑150 (MQ=255)
                                                                                                    tCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTc                                                <  1:1422623/150‑1 (MQ=255)
                                                                                                      aaaGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTccc                                              >  1:204631/1‑150 (MQ=255)
                                                                                                       aaGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCt                                             >  1:840249/1‑150 (MQ=255)
                                                                                                        aGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTg                                            <  2:1367451/150‑1 (MQ=255)
                                                                                                            tGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGCCCCCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCaaa                                        >  2:237660/1‑150 (MQ=255)
                                                                                                            tGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCaaa                                        >  1:1303627/1‑150 (MQ=255)
                                                                                                             gATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCaaaa                                       >  2:1457251/1‑150 (MQ=255)
                                                                                                                 gAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACggg                                   >  2:98299/1‑150 (MQ=255)
                                                                                                                  aaaaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACgggg                                  <  1:750536/150‑1 (MQ=255)
                                                                                                                  aaaaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACgggg                                  >  1:75904/1‑150 (MQ=255)
                                                                                                                  aaaaaCTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACgggg                                  <  2:544756/150‑1 (MQ=255)
                                                                                                                          gcgcTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTaaaaa                          >  1:28188/1‑150 (MQ=255)
                                                                                                                           cgcTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCTACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTaaaaaa                         >  2:1219777/1‑150 (MQ=255)
                                                                                                                           cgcTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTaaaaaa                         >  2:485120/1‑150 (MQ=255)
                                                                                                                            gcTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAAt                        <  1:173898/150‑1 (MQ=255)
                                                                                                                             cTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATc                       <  2:30365/150‑1 (MQ=255)
                                                                                                                               gAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGCAGGCGCCCCCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTCGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTg                     >  1:694310/1‑150 (MQ=255)
                                                                                                                                  cAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGgag                  <  2:1182374/150‑1 (MQ=255)
                                                                                                                                       aTGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGaaaaa             >  2:1538064/1‑150 (MQ=255)
                                                                                                                                        tGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAaat            <  1:1263937/150‑1 (MQ=255)
                                                                                                                                        tGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCAGGAGAAAaat            >  2:636486/1‑150 (MQ=255)
                                                                                                                                        gccttcagtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAaat            >  1:519262/8‑150 (MQ=255)
                                                                                                                                         ccttcagtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAaata           >  1:1442077/7‑150 (MQ=255)
                                                                                                                                          cttcagtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTgg                    <  2:1279610/135‑1 (MQ=255)
                                                                                                                                          cttcagtgCTGATAACCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTgg                    >  1:1279610/6‑140 (MQ=255)
                                                                                                                                                tgCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTa    <  1:1457251/150‑1 (MQ=255)
                                                                                                                                                  cTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAg  >  2:950740/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCTGGGTAATCCGTACGGCAACGACCCGCGCATCATCTCCGGTGAATCCGGCGCTGTCGGTTTGGGCGTTCTCGCAGCGGTTCATTATCACCCGCAACGTCAAAGCCTGATGGAAAAACTGGCGCTGAACAAAGATGCCGTGGTGCTGGTTATCAGCACTGAAGGCGACACCGACGTGAAGCACTACCGCGAAGTTGTCTGGGAAGGCAAACACGCTGTAGCACCTTAATCACCCTATTGGAACTGCTCCCTGCAAAACGGGGAGTAAAAAATCTGGAGAAAAATAATGGCTAAG  >  NC_000913/3008477‑3008772

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