Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,147,279 T→A pseudogene (413/420 nt) yghX ← pseudogene, dienlactone hydrolase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,147,2790TA98.7% 294.7 / ‑4.3 78pseudogene (413/420 nt)yghXpseudogene, dienlactone hydrolase family
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (44/32);  minor base G (1/0);  total (45/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGATCAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGATTC  >  NC_000913/3147137‑3147425
                                                                                                                                              |                                                                                                                                                  
ccAGTTCCGAATACGCAAAAGCCGCCGCTTGCGATGCGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACct                                                                                                                                              <  1:1194182/149‑1 (MQ=255)
  aGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACctgctg                                                                                                                                          >  2:996235/1‑151 (MQ=255)
   gTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGt                                                                                                                                         >  2:318180/1‑151 (MQ=255)
      ccGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTg                                                                                                                                       <  2:564444/150‑1 (MQ=255)
      ccGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGc                                                                                                                                      >  1:1077244/1‑151 (MQ=255)
        ggATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGc                                                                                                                                      >  2:506279/1‑149 (MQ=255)
        ggATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAg                                                                                                                                    <  2:378462/151‑1 (MQ=255)
             cgcgACAGCCGCCGCGTTCGATCCGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCAc                                                                                                                               <  2:470180/151‑1 (MQ=255)
                      cgccgcGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTat                                                                                                                      >  1:285597/1‑151 (MQ=255)
                          gcgTTCGATACTCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCAGTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTatcat                                                                                                                   <  2:1290155/150‑1 (MQ=255)
                          gcgTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTatcatc                                                                                                                  >  1:693759/1‑151 (MQ=255)
                          gcgTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTatcatc                                                                                                                  >  2:10174/1‑151 (MQ=255)
                             ttCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCAtt                                                                                                               >  1:679665/1‑151 (MQ=255)
                             ttCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCAtt                                                                                                               >  2:403277/1‑151 (MQ=255)
                             ttCAATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTTGAACAACCTGCTGTTGCAGCTCACGACCTTTAtcatca                                                                                                                 <  1:306798/149‑1 (MQ=255)
                                    cGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCgg                                                                                                         <  1:1164876/150‑1 (MQ=255)
                                        aCCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATaa                                                                                                     <  1:1228041/150‑1 (MQ=255)
                                         ccGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAAc                                                                                                    <  2:146966/150‑1 (MQ=255)
                                          cGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCt                                                                                                  >  2:854500/1‑151 (MQ=255)
                                           gCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTc                                                                                                 <  1:1049652/151‑1 (MQ=255)
                                             cATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCa                                                                                               <  2:252057/151‑1 (MQ=255)
                                              aTAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCaa                                                                                              <  2:513787/151‑1 (MQ=255)
                                               tAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAAc                                                                                             >  1:457976/1‑151 (MQ=255)
                                                aGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCaaag                                                                                            <  2:934087/151‑3 (MQ=255)
                                                 gCAAAATCCGGTAATAGCCACTGTGACTGGCGCTTGCGGATCGCGTTGCATAAACTCAATTGCGGCAAAGAATTCATTTATCAGTTTGTTGGAAAAAGCTGCTGTTGCCGCTAAAGATCTTTATCATCAGTTTCCGGATAAACTTCAAc                                                                                             >  1:42249/1‑149 (MQ=255)
                                                 gCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACgg                                                                                           >  2:1168893/1‑151 (MQ=255)
                                                   aaaaTCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGa                                                                                          >  2:281867/1‑150 (MQ=255)
                                                          ggTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGc                                                                                  >  2:1081874/1‑151 (MQ=255)
                                                            tAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCg                                                                                >  1:983570/1‑151 (MQ=255)
                                                            tAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGGTGGAACAACCTGCTGTTGCAGCTCACGACCTGTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCg                                                                                >  1:505102/1‑151 (MQ=255)
                                                              aTACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTc                                                                              <  2:528415/151‑1 (MQ=255)
                                                                aCCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTc                                                                              >  2:719356/1‑149 (MQ=255)
                                                                 cccACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAgg                                                                           <  1:922720/151‑1 (MQ=255)
                                                                     gtttGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGcc                                                                       >  1:701555/2‑151 (MQ=255)
                                                                     cTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGcc                                                                       >  2:1242880/1‑151 (MQ=255)
                                                                          ccTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGc                                                                  >  2:673231/1‑151 (MQ=255)
                                                                            tGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCg                                                                 >  1:79392/1‑150 (MQ=255)
                                                                              tCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGATATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGat                                                               >  1:427840/1‑150 (MQ=255)
                                                                                 cTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAc                                                           >  1:61298/1‑151 (MQ=255)
                                                                                   tGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAccc                                                         >  2:1249027/1‑151 (MQ=255)
                                                                                   tGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAcc                                                          <  1:798750/150‑1 (MQ=255)
                                                                                   tGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAcc                                                          <  2:1136292/150‑1 (MQ=255)
                                                                                    gCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAcccc                                                        >  1:1291166/1‑151 (MQ=255)
                                                                                    gCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATAcccc                                                        <  1:277853/151‑1 (MQ=255)
                                                                                     cGGATAGCGTTGCATAAACGCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCg                                                       <  1:664896/151‑1 (MQ=255)
                                                                                      ggATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGc                                                      >  2:81922/1‑151 (MQ=255)
                                                                                        aTAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCt                                                    >  2:469038/1‑151 (MQ=255)
                                                                                          aGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCtt                                                   <  2:417238/150‑1 (MQ=255)
                                                                                          aGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTc                                                  >  2:25353/1‑151 (MQ=255)
                                                                                            cGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGc                                                >  1:256052/1‑151 (MQ=255)
                                                                                                    aaaCTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGATCAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAGGCCGTCAGGTTCCAGGGCGATATACCCCGCCTTCGCCACTCGcc                                        <  1:969287/151‑1 (MQ=255)
                                                                                                    aaaCTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGcc                                        >  1:230790/1‑151 (MQ=255)
                                                                                                      aCTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGt                                      <  1:588735/151‑1 (MQ=255)
                                                                                                       cTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTg                                     >  2:772481/1‑151 (MQ=255)
                                                                                                          aaTTGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGcc                                   >  2:1132468/1‑150 (MQ=255)
                                                                                                            ttGCGGCAAAGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCcac                                 >  2:1222052/1‑150 (MQ=255)
                                                                                                                   aaaGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCg                          >  1:662001/1‑150 (MQ=255)
                                                                                                                    aaGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGat                        <  2:51662/151‑1 (MQ=255)
                                                                                                                    aaGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGat                        >  2:1115477/1‑151 (MQ=255)
                                                                                                                     aGAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGat                        >  2:1269694/1‑150 (MQ=255)
                                                                                                                      gAAATCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGatat                      >  1:364550/1‑151 (MQ=255)
                                                                                                                       aaaTCATTCATCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGAtata                     <  1:422761/151‑1 (MQ=255)
                                                                                                                         aTCATTCATCAGTTTGTTGGAGCAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACg                   >  2:1220333/1‑151 (MQ=255)
                                                                                                                           cattcatCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGa                 <  1:489010/151‑1 (MQ=255)
                                                                                                                             ttcatCAGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGAtt               >  1:345653/1‑151 (MQ=255)
                                                                                                                               catcaGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTc              <  1:898286/150‑1 (MQ=255)
                                                                                                                                atcaGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTc              >  1:353226/1‑149 (MQ=255)
                                                                                                                                  caGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCGTCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATGCGGATTCAGt           <  1:117604/150‑1 (MQ=255)
                                                                                                                                  caGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGt           <  2:983609/150‑1 (MQ=255)
                                                                                                                                   aGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTcc         >  2:394725/1‑151 (MQ=255)
                                                                                                                                   aGTTTGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTcc         <  2:56721/151‑1 (MQ=255)
                                                                                                                                     tttGTTGGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCAc       >  2:780073/1‑151 (MQ=255)
                                                                                                                                       tgttgGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACg      >  2:106792/1‑150 (MQ=255)
                                                                                                                                       tgttgGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACg      <  1:826880/150‑1 (MQ=255)
                                                                                                                                        gttgGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGAt    <  1:146845/151‑1 (MQ=255)
                                                                                                                                         ttgGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGAtt   >  1:1004780/1‑151 (MQ=255)
                                                                                                                                          tgGAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGAtt   <  1:226504/150‑1 (MQ=255)
                                                                                                                                           ggAACAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGATtc  <  2:577576/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCAGTTCCGGATACGCGACAGCCGCCGCGTTCGATACGCCACCGCCATAGCAAAATCCGGTAATACCCACTTTGCCTGTCGCTTGCGGATAGCGTTGCATAAACTCAATTGCGGCAAAGAAATCATTCATCAGTTTGTTGGATCAACCTGCTGTTGCAGCTCACGACCTTTATCATCATTTCCCGGATAACCTCCAACGGAACTTAAGCCGTCAGGTGCCAGGGCGATATACCCCGCCTTCGCCACTCGCCGTGCCACATCTTCGATATACGGATTCAGTCCACGATTC  >  NC_000913/3147137‑3147425

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