Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,335,787 A→G G236G (GGT→GGC)  murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,335,7870AG98.7% 279.4 / ‑3.9 79G236G (GGT→GGC) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (45/32);  total (46/33)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

ATGCATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCA  >  NC_000913/3335640‑3335934
                                                                                                                                                   |                                                                                                                                                   
aTGCATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCg                                                                                                                                                  >  1:847162/1‑151 (MQ=255)
 tGCATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCgg                                                                                                                                                 >  1:149001/1‑151 (MQ=255)
  gCATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGt                                                                                                                                                <  1:865385/151‑1 (MQ=255)
   cATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGtt                                                                                                                                               <  2:570638/151‑1 (MQ=255)
    atatCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGttt                                                                                                                                              <  1:1038140/151‑1 (MQ=255)
    atatCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGttt                                                                                                                                              <  1:681159/151‑1 (MQ=255)
        ccAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGa                                                                                                                                           <  2:715445/150‑1 (MQ=255)
        ccAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGAt                                                                                                                                          >  2:1108383/1‑151 (MQ=255)
         cAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTACCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGAt                                                                                                                                          >  1:1182579/1‑150 (MQ=255)
             cTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGAt                                                                                                                                     >  2:462650/1‑151 (MQ=255)
              tAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATc                                                                                                                                    >  1:924263/1‑151 (MQ=255)
                   cAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCa                                                                                                                               <  2:447890/151‑1 (MQ=255)
                    aGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGc                                                                                                                                >  1:1029032/1‑149 (MQ=255)
                        ttCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGAACCGGCAGAAc                                                                                                                           >  1:689648/1‑150 (MQ=255)
                           gCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGa                                                                                                                       <  2:1101353/151‑1 (MQ=255)
                            ccGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGa                                                                                                                       <  1:713046/150‑1 (MQ=255)
                             cGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGAt                                                                                                                      >  1:662487/1‑150 (MQ=255)
                              gACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAg                                                                                                                    <  2:1216334/151‑1 (MQ=255)
                               aCTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGa                                                                                                                   <  2:1064990/151‑1 (MQ=255)
                                cTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGa                                                                                                                   >  2:603327/1‑150 (MQ=255)
                                  tCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGacac                                                                                                                >  2:796740/1‑151 (MQ=255)
                                    gATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACAcc                                                                                                               >  1:1058077/1‑150 (MQ=255)
                                       gTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACAccgccg                                                                                                           >  1:36426/1‑151 (MQ=255)
                                       gTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACAccgcc                                                                                                            <  2:293318/150‑1 (MQ=255)
                                         ccGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCgccgcc                                                                                                         >  2:1258146/1‑151 (MQ=255)
                                         ccGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCgccgcc                                                                                                         >  2:44741/1‑151 (MQ=255)
                                          cGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCt                                                                                                        >  1:1244564/1‑151 (MQ=255)
                                           gCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCt                                                                                                        >  2:909456/1‑150 (MQ=255)
                                              ccAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCCCCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTaaa                                                                                                     <  1:14227/150‑1 (MQ=255)
                                              ccAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAAc                                                                                                    <  1:351504/151‑1 (MQ=255)
                                              ccAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATACGGCAGAACGCGATAGACACCGCCGCCTAAAc                                                                                                    >  1:1074780/1‑151 (MQ=255)
                                                  cGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTc                                                                                                >  2:131632/1‑151 (MQ=255)
                                                  cGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTc                                                                                                >  1:106374/1‑151 (MQ=255)
                                                     cACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCcac                                                                                             >  2:768182/1‑151 (MQ=255)
                                                       cGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCcaca                                                                                            <  2:1230254/150‑1 (MQ=255)
                                                       cGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCacac                                                                                           <  1:658487/151‑1 (MQ=255)
                                                        gCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACAcc                                                                                          >  2:614652/1‑151 (MQ=255)
                                                         cAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCt                                                                                         >  1:1074679/1‑151 (MQ=255)
                                                         cAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCt                                                                                         <  2:78692/151‑1 (MQ=255)
                                                                  cAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGacg                                                                                >  2:669170/1‑151 (MQ=255)
                                                                    gCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGacga                                                                               >  2:276107/1‑150 (MQ=255)
                                                                      aCGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGacgata                                                                             <  1:841979/150‑1 (MQ=255)
                                                                        ggCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGacgatacg                                                                           >  1:6865/1‑150 (MQ=255)
                                                                         gCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGacgatacgat                                                                         >  2:1264323/1‑151 (MQ=255)
                                                                               agagTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAAGCGCCGCCGCCACCAGGCAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATAGGTGc                                                                   >  1:1286455/1‑151 (MQ=255)
                                                                                  gTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGccca                                                                <  1:760681/151‑2 (MQ=255)
                                                                                      cTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGAcc                                                            >  1:677964/1‑151 (MQ=255)
                                                                                        ggCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCg                                                           <  2:1102562/150‑1 (MQ=255)
                                                                                        ggCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCg                                                           <  2:143244/150‑1 (MQ=255)
                                                                                           tGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTa                                                        >  2:295461/1‑150 (MQ=255)
                                                                                                 ttACGGGTGATTATTTTGCCGCGAGAAATCGCCGCCGCCCTGAGGAATGGGCCGGTTTCGGTACGACCCGGCAAAACGCGATAGACAGCGACCCCTAAAAGTTCCACACATTCGATGACGATACGATCGGTGACCTGAACGCAAAtttt                                                   >  1:1083609/1‑149 (MQ=255)
                                                                                                 ttACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAACGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTc                                                  >  1:1035933/1‑150 (MQ=255)
                                                                                                   aCGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGc                                                >  2:1104731/1‑150 (MQ=255)
                                                                                                   aCGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGc                                                <  2:327301/150‑1 (MQ=255)
                                                                                                   aCGGCAGATAATTTTGCCGAGAGAAATCGACGCCGCCACCAGGAAAGTCCCGGTTTCGATAAGATCCGGCAGAACGCGACAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGc                                                >  1:1268608/1‑150 (MQ=255)
                                                                                                      gCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCAccc                                            <  2:687438/151‑1 (MQ=255)
                                                                                                          aTAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCg                                        >  2:965471/1‑151 (MQ=255)
                                                                                                            aaTTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGt                                       <  2:850731/150‑1 (MQ=255)
                                                                                                            aaTTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTa                                      <  2:1282965/151‑1 (MQ=255)
                                                                                                                ttGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAATCGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATc                                   <  2:338225/150‑1 (MQ=255)
                                                                                                                   ccGCGAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAgg                                >  2:187314/1‑150 (MQ=255)
                                                                                                                    cgcgAGAAATCGCCGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGa                               >  2:564421/1‑150 (MQ=255)
                                                                                                                       gagaAATCGACGCCGCCACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCACCGCCTAAACGTTCCACACCTTCGCTGACGATACGCTCGGTGCCCTGACCGCTAATTTTCGCACACAGCGTGATCAGGAAGtt                           >  1:322196/1‑151 (MQ=255)
                                                                                                                                ccgccgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATc                  >  1:766902/1‑151 (MQ=255)
                                                                                                                                 cgccgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATc                  >  2:1026108/1‑150 (MQ=255)
                                                                                                                                   ccgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATcgac               <  2:522153/151‑1 (MQ=255)
                                                                                                                                    cgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATcgacg              >  1:265529/1‑151 (MQ=255)
                                                                                                                                    cgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATcgacg              <  2:1250154/151‑1 (MQ=255)
                                                                                                                                    cgccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATcgacg              <  1:1068424/151‑1 (MQ=255)
                                                                                                                                     gccACCAGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATcgacga             <  2:720848/151‑1 (MQ=255)
                                                                                                                                      ccaccaGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGAt            <  1:1249199/151‑1 (MQ=255)
                                                                                                                                       caccaGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGAtt           >  1:1095326/1‑151 (MQ=255)
                                                                                                                                        accaGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGAttt          >  1:461452/1‑151 (MQ=255)
                                                                                                                                        accaGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGAttt          >  2:920325/1‑151 (MQ=255)
                                                                                                                                        accaGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGAttt          <  1:808080/151‑1 (MQ=255)
                                                                                                                                           aGGAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTcc        <  1:472669/150‑1 (MQ=255)
                                                                                                                                            ggAAAGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGGGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCgg      >  1:861667/1‑151 (MQ=255)
                                                                                                                                                aGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTc   >  2:524110/1‑150 (MQ=255)
                                                                                                                                                aGTGCCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCa  <  2:6862/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATGCATATCCAGGCTAATCCAGTCTTCGCCGACTTCGATGTCCGCTCCAGCGTCACGCAGTTTCGCCAGCACGGCGTCGAGAGTATCTGGCTGCGCGTTACGGCAGATAATTTTGCCGCGAGAAATCGCCGCCGCCACCAGGAAAGTACCGGTTTCGATACGATCCGGCAGAACGCGATAGACACCGCCGCCTAAACGTTCCACACCTTCGATGACGATACGATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCA  >  NC_000913/3335640‑3335934

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