Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,072,800 C→T 6.6% D298N (GAT→AAT)  epd ← D‑erythrose‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,072,8000CT 6.6% 281.3 / 6.6 91D298N (GAT→AAT) epdD‑erythrose‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base T (3/3);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

GCAGACGCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCACCGT  >  NC_000913/3072654‑3072944
                                                                                                                                                  |                                                                                                                                                
gCAGACGCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGAc                                                                                                                                               >  1:44532/1‑150 (MQ=255)
gCAGACGCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCAACAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGAc                                                                                                                                               >  2:56761/1‑150 (MQ=255)
    aCGCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATg                                                                                                                                           <  1:882666/150‑1 (MQ=255)
         gctTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCGAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAcc                                                                                                                                      <  2:85980/150‑2 (MQ=255)
          ctTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCACCGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACtg                                                                                                                                     <  1:1446327/150‑1 (MQ=255)
           ttGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAccca                                                                                                                                    <  1:137471/150‑4 (MQ=255)
            tGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATTGTCGATGGCAcccaa                                                                                                                                   <  1:1209041/150‑5 (MQ=255)
                 cTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGAt                                                                                                                              <  1:577911/150‑1 (MQ=255)
                  ttACCTGAAAGCAACAGTAGCGATAGCTAACGTCGTGTCGAGGATTCGGTTAGCAAAGCCCGATTCGGTATCGCACGAGACCATCGTTTTGATCAGGTGTGCGCCTCTGACCGGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATc                                                                                                                             <  1:235723/150‑1 (MQ=255)
                   tACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCg                                                                                                                            >  1:813582/1‑150 (MQ=255)
                   tACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCg                                                                                                                            >  2:1433307/1‑150 (MQ=255)
                   tACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCg                                                                                                                            >  1:426328/1‑150 (MQ=255)
                     ccTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTg                                                                                                                          <  1:63294/150‑1 (MQ=255)
                     ccTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTg                                                                                                                          >  2:121279/1‑150 (MQ=255)
                      cTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTgg                                                                                                                         <  1:865136/150‑1 (MQ=255)
                        gAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGtt                                                                                                                       >  2:1406582/1‑150 (MQ=255)
                        gAAAGCAACAGTAGCCATAGCTAACGTCCTGTCGAGCATTCGCTTAGCAAAGCCCCATTCGTTATAGCACCAGCCCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATACACAATGGCACTGTGCGGATCGTGGtt                                                                                                                       >  2:573521/1‑150 (MQ=255)
                         aaaGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCg                                                                                                                                 <  2:1584621/139‑1 (MQ=255)
                         aaaGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCg                                                                                                                                 >  1:1584621/1‑139 (MQ=255)
                         aaaGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTa                                                                                                                      <  1:1214128/150‑1 (MQ=255)
                              aaCAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATc                                                                                                                 <  1:1451526/150‑1 (MQ=255)
                               aCAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCt                                                                                                                <  2:465225/150‑1 (MQ=255)
                                cAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTa                                                                                                               <  2:1420105/150‑1 (MQ=255)
                                 aGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTAc                                                                                                              <  2:974604/150‑1 (MQ=255)
                                  gTAGCCATAGCTAAAGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACCTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACa                                                                                                             >  2:150033/1‑150 (MQ=255)
                                     gCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACagag                                                                                                          >  1:1006635/1‑150 (MQ=255)
                                     gCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACagag                                                                                                          <  2:84299/150‑1 (MQ=255)
                                        aTAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGAcc                                                                                                       <  1:775542/150‑1 (MQ=255)
                                            cTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAAca                                                                                                   >  2:375527/1‑149 (MQ=255)
                                            cTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTGTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACg                                                                                                   <  1:812391/150‑1 (MQ=255)
                                            cTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTAGTTATCGCACCAGACCAACGTTTTGATAAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACg                                                                                                   >  1:153962/1‑150 (MQ=255)
                                                 gtcgtGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAAt                                                                                              >  2:697537/1‑150 (MQ=255)
                                                     tgtCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCg                                                                                          >  1:368217/1‑150 (MQ=255)
                                                      gtCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGt                                                                                         <  2:198655/150‑1 (MQ=255)
                                                      gtCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGt                                                                                         >  2:880985/1‑150 (MQ=255)
                                                       tCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGta                                                                                        <  2:206457/150‑1 (MQ=255)
                                                         gAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGtata                                                                                      >  1:916110/1‑150 (MQ=255)
                                                         gAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGtata                                                                                      >  1:656726/1‑150 (MQ=255)
                                                          aGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAg                                                                                     <  2:35340/150‑1 (MQ=255)
                                                          aGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAg                                                                                     <  1:1428539/150‑1 (MQ=255)
                                                           gCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGt                                                                                    >  1:79665/1‑150 (MQ=255)
                                                             aTTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATTGCACTGTGCGGATCGTTGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCa                                                                                  <  2:44532/150‑1 (MQ=255)
                                                               tCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAAc                                                                                <  2:426328/150‑1 (MQ=255)
                                                                cGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACt                                                                               >  1:923504/1‑150 (MQ=255)
                                                                     aGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATAcc                                                                          >  2:1037770/1‑150 (MQ=255)
                                                                      gCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGGTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCa                                                                         >  2:945611/1‑150 (MQ=255)
                                                                        aaaGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATg                                                                       >  1:828872/1‑150 (MQ=255)
                                                                           gCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGaaa                                                                    >  2:1316185/1‑150 (MQ=255)
                                                                               cATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCa                                                                <  1:721511/150‑1 (MQ=255)
                                                                                   cGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCtt                                                            <  2:1004431/150‑1 (MQ=255)
                                                                                    gTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAGCTATACCATGAAATGCACCTtg                                                           <  2:320956/150‑1 (MQ=255)
                                                                                    gTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTtg                                                           <  2:368217/150‑1 (MQ=255)
                                                                                    gTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTtg                                                           <  2:1622158/150‑1 (MQ=255)
                                                                                    gTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCAGGGTTTGGGTGCCATCGACAATGGCACTCTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTtg                                                           >  2:1488454/1‑150 (MQ=255)
                                                                                       aTCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGtgc                                                        <  2:996940/150‑1 (MQ=255)
                                                                                           cACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTgctgct                                                    >  1:1100248/1‑150 (MQ=255)
                                                                                              cAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCtttt                                                 <  1:1406582/150‑1 (MQ=255)
                                                                                                 aCCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgc                                              <  1:415371/150‑1 (MQ=255)
                                                                                                  ccAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGTCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgca                                             <  2:809060/150‑1 (MQ=255)
                                                                                                  ccAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgca                                             <  2:583501/150‑1 (MQ=255)
                                                                                                   cAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgcag                                            >  2:24846/1‑150 (MQ=255)
                                                                                                     aCGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACAGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTgcagca                                          <  1:375527/150‑1 (MQ=255)
                                                                                                         tttGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAg                                      <  1:128904/150‑1 (MQ=255)
                                                                                                          ttGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAgg                                     >  1:1487354/1‑150 (MQ=255)
                                                                                                           tGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGt                                    <  1:1002774/150‑1 (MQ=255)
                                                                                                             aTCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTg                                  >  1:692225/1‑150 (MQ=255)
                                                                                                               cAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGAc                                <  1:386889/150‑1 (MQ=255)
                                                                                                                aGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACt                               <  2:327356/150‑1 (MQ=255)
                                                                                                                aGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACt                               >  2:877339/1‑150 (MQ=255)
                                                                                                                 ggTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACtt                              >  1:1070470/1‑150 (MQ=255)
                                                                                                                     tgCGTCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCAtt                          >  2:257343/1‑150 (MQ=255)
                                                                                                                     tgCGTCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCAtt                          <  1:257343/150‑1 (MQ=255)
                                                                                                                      gcgcCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTg                         <  2:337388/150‑1 (MQ=255)
                                                                                                                       cgcCACTGCCCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTgg                        >  1:1080533/1‑150 (MQ=255)
                                                                                                                        gcCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGc                       >  1:742214/1‑150 (MQ=255)
                                                                                                                         ccACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCt                      >  1:531811/1‑150 (MQ=255)
                                                                                                                          cACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCtt                     <  2:153962/150‑1 (MQ=255)
                                                                                                                              gACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCt                                                             >  1:705093/1‑106 (MQ=255)
                                                                                                                              gACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCt                                                             <  2:705093/106‑1 (MQ=255)
                                                                                                                              gACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTAc                 <  1:246354/150‑1 (MQ=255)
                                                                                                                                cccGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTTTAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCtttt                   <  2:1376447/146‑1 (MQ=255)
                                                                                                                                cccGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTGTAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCtttt                   >  1:1376447/1‑146 (MQ=255)
                                                                                                                                cccGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAg               <  2:813582/150‑1 (MQ=255)
                                                                                                                                cccGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAg               >  1:1054268/1‑150 (MQ=255)
                                                                                                                                 ccGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAgg              <  1:697537/150‑1 (MQ=255)
                                                                                                                                      gcacagTGCCATTGTCGATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTttct         >  1:1495737/6‑150 (MQ=255)
                                                                                                                                         cagTGCCATTGTCGATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCa      >  1:1360142/3‑150 (MQ=255)
                                                                                                                                         cagTGCCATTGTCGATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCa      >  1:574163/3‑150 (MQ=255)
                                                                                                                                           gTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCAcc    >  2:831503/1‑150 (MQ=255)
                                                                                                                                             gCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCACCGt  >  1:929500/1‑150 (MQ=255)
                                                                                                                                             gCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCACCGt  <  2:531811/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCAGACGCTGCTTGCGTCTTACCTGAAAGCAACAGTAGCCATAGCTAACGTCGTGTCGAGCATTCGGTTAGCAAAGCCCCATTCGTTATCGCACCAGACCAACGTTTTGATCAGGTGTGCGCCACTGACCCGGGTTTGGGTGCCATCGACAATGGCACTGTGCGGATCGTGGTTAAAATCTACAGAGACCAACGGCAATTCCGTATAGTCAACTATACCATGAAATGCACCTTGTGCTGCTTTTTGCAGCAACAGGTTGACTTCATTGGCTTTTACAGGTTTCTTCACCGT  >  NC_000913/3072654‑3072944

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