Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,219,696 C→T G367G (GGC→GGT aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,219,6960CT98.4% 241.3 / ‑3.6 64G367G (GGC→GGTaceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/39);  minor base A (0/1);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

GAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGA  >  NC_000913/4219554‑4219842
                                                                                                                                              |                                                                                                                                                  
gAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGc                                                                                                                                            <  2:195748/151‑1 (MQ=255)
 aaTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCg                                                                                                                                           >  1:1038640/1‑151 (MQ=255)
 aaTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCg                                                                                                                                           <  1:50726/151‑1 (MQ=255)
  aTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCgg                                                                                                                                          <  1:1016322/151‑1 (MQ=255)
  aTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCgg                                                                                                                                          >  2:433461/1‑151 (MQ=255)
  aTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCgg                                                                                                                                          >  1:1061245/1‑151 (MQ=255)
   tGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCgg                                                                                                                                          >  1:731816/1‑150 (MQ=255)
     gTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGaca                                                                                                                                       <  1:260682/151‑1 (MQ=255)
       gATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACAcc                                                                                                                                     <  1:865155/151‑1 (MQ=255)
          gTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCa                                                                                                                                   <  1:500355/150‑1 (MQ=255)
            tGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAgg                                                                                                                                 <  1:682799/150‑1 (MQ=255)
             ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAg                                                                                                                               <  2:121818/151‑1 (MQ=255)
                gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGttt                                                                                                                            >  2:21404/1‑151 (MQ=255)
                gtTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGAGTCGAATGGCGGACACCCAGGAGttt                                                                                                                            <  1:334562/151‑1 (MQ=255)
                  ttACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTg                                                                                                                           <  2:997880/150‑1 (MQ=255)
                    aCGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGaa                                                                                                                         >  2:273983/1‑150 (MQ=255)
                            gggCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTtgtg                                                                                                                <  2:978049/151‑1 (MQ=255)
                                tttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTgg                                                                                                            >  1:538182/1‑151 (MQ=255)
                                 ttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTgg                                                                                                            <  1:103318/150‑1 (MQ=255)
                                 ttGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTg                                                                                                             <  1:308025/149‑1 (MQ=255)
                                      cGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGc                                                                                                      >  1:848581/1‑151 (MQ=255)
                                       gggTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCGCGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGc                                                                                                      <  2:726034/150‑1 (MQ=255)
                                         gTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGc                                                                                                   <  2:478758/151‑1 (MQ=255)
                                           aTTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCat                                                                                                 >  2:885192/1‑151 (MQ=255)
                                              cAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATAtt                                                                                              <  2:695822/151‑1 (MQ=255)
                                                 aGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTcc                                                                                           >  1:524793/1‑151 (MQ=255)
                                                    catcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCg                                                                                        <  2:416477/151‑1 (MQ=255)
                                                    catcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGACGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCg                                                                                        <  2:1010607/151‑1 (MQ=255)
                                                      tcaAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGc                                                                                      >  1:551512/1‑151 (MQ=255)
                                                       caAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGc                                                                                      >  2:771763/1‑150 (MQ=255)
                                                       caAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCa                                                                                     <  2:188699/151‑1 (MQ=255)
                                                       caAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCa                                                                                     <  2:144380/151‑1 (MQ=255)
                                                        aaaGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCAt                                                                                    >  1:157307/1‑151 (MQ=255)
                                                         aaGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCAtt                                                                                   <  2:150208/151‑1 (MQ=255)
                                                          aGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTa                                                                                  <  2:567703/151‑1 (MQ=255)
                                                              aGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATgg                                                                              >  2:886851/1‑151 (MQ=255)
                                                                     cgccgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACt                                                                       >  1:717090/1‑151 (MQ=255)
                                                                     cgccgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACt                                                                       <  1:112060/151‑1 (MQ=255)
                                                                       ccgcAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGc                                                                     <  2:661688/151‑1 (MQ=255)
                                                                           aGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCa                                                                 <  2:663834/151‑1 (MQ=255)
                                                                             aaaGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAgg                                                               <  2:615820/151‑1 (MQ=255)
                                                                              aaGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAgg                                                               >  1:1054603/1‑150 (MQ=255)
                                                                                    tGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAagcagc                                                        >  1:110309/1‑151 (MQ=255)
                                                                                     gTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAagcagc                                                        <  2:255111/150‑1 (MQ=255)
                                                                                     gTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCg                                                       <  1:313431/151‑1 (MQ=255)
                                                                                         gCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaaa                                                   <  1:6222/151‑1 (MQ=255)
                                                                                         gCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATAATTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGa                                                     <  1:355865/149‑1 (MQ=255)
                                                                                          ccGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaaa                                                   >  2:103546/1‑150 (MQ=255)
                                                                                           cGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGaaaaa                                                 <  1:243843/151‑1 (MQ=255)
                                                                                               cACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCa                                             <  2:714249/151‑1 (MQ=255)
                                                                                                 cGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCAcc                                           >  1:612442/1‑151 (MQ=255)
                                                                                                  gTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCg                                          <  1:520766/151‑1 (MQ=255)
                                                                                                   ctCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCg                                          >  1:1009469/2‑150 (MQ=255)
                                                                                                    tCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGa                                         <  2:479468/150‑1 (MQ=255)
                                                                                                     cGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGAtc                                       <  2:722152/151‑1 (MQ=255)
                                                                                                         ccTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCgg                                   <  2:672210/151‑1 (MQ=255)
                                                                                                           tGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCg                                 >  1:866699/1‑151 (MQ=255)
                                                                                                                 cAACTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACaa                            <  1:602353/150‑1 (MQ=255)
                                                                                                                  aaCTGGTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAAt                          >  1:30731/1‑151 (MQ=255)
                                                                                                                      ggTGAAAGAGCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATtgt                       <  2:418301/150‑1 (MQ=255)
                                                                                                                            agagCACGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGc                >  1:987781/1‑151 (MQ=255)
                                                                                                                                 aCGATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCAt             <  2:1061422/149‑1 (MQ=255)
                                                                                                                                   gATCGCGTGGGTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCttt         >  2:855782/1‑151 (MQ=255)
                                                                                                                                           gggTCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGa  <  1:623798/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGA  >  NC_000913/4219554‑4219842

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