Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,912,074 G→A S532F (TCC→TTC)  glmS ← L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,912,0740GA97.7% 135.3 / NA 45S532F (TCC→TTC) glmSL‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (29/14);  minor base C (1/0);  total (30/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCC  >  NC_000913/3911835‑3912316
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
tttATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAATTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCAAGTTATCGCTACTTCCAAAACCCGCATCCTGATCGGCGAAGACATACCACTGAACGCCACGCGCGCGACCTTCTTCACTGTTGAATTTCAGtttt                                                                                                                                                                                                                                         >  2:451909/1‑251 (MQ=255)
     ttcccAACCTTAACCGTTTTTTCCAGGTTACGCGGCGGGTCACCGTCGGGGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGAGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCa                                                                                                                                                                                                                                     <  2:458285/246‑1 (MQ=255)
          aaCCGTACCCGATTGTGCCAGGGTACGCGGCGGGTCAACGTCGGGGCCTTTTATCGGCTCGACATGGTAATCCAGCAGTTGTAGCGGAACGGTTTAGAAAAGCGGTGCAATCCCCTCTTCCACATGCGGCAGCACGATGAGGTGCATGTTATCGCTACTTACATATCCCGCATCCTGATCGGCGAAGACATACAACGTTCCGCCACGCGCGCGACCGGCAGCAATGTTGAAGTTCAGTCTGTCCAGCAAtt                                                                                                                                                                                                                               <  2:277630/251‑1 (MQ=255)
             cGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgt                                                                                                                                                                                                                            >  1:152897/1‑251 (MQ=255)
              gTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCTGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgtt                                                                                                                                                                                                                           >  1:45627/1‑251 (MQ=255)
                    gATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCgg                                                                                                                                                                                                                     >  1:462389/1‑251 (MQ=255)
                           ccAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaac                                                                                                                                                                                                               >  1:376092/1‑250 (MQ=255)
                                           gTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGGTGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATAAAACTGACCGCCACGCCCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCATTTCGTTTTTCGGTGCAACAACAATAACCGGCatat                                                                                                                                                                                              >  2:370295/1‑251 (MQ=255)
                                                cGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGc                                                                                                                                                                                          >  1:197691/1‑250 (MQ=255)
                                                    ggTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGCACGGTGTAGAAGATCGGTGCAATCACATCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCa                                                                                                                                                                                      >  1:266006/1‑250 (MQ=255)
                                                            gtaacagagcaataTGGAAAGCCAGTAGTTGCAGCTGAACTGTGTAGAAGATCGGTGCAACCCCCACTTCCACATGCTGCATCGAGTGTATGTGCATGTTATCGCTACTTCCAAAACCCGCTACCTGATTGGCGAAGTCAAACAACTTTCCGTCACGCGCTCGAACTTCTTCAATGCTGAATTTCAGCTCTTCCAGCAATTCGGTGTCCGGTGCAACAACAATAACCGGCATATCGGCAGTAATTAGCGcc                                                                                                                                                                             <  2:69814/238‑1 (MQ=255)
                                                                           ggTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTACAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGACATCAATTAGAGCCAGCGGACCGTGTTTc                                                                                                                                                              >  1:261120/1‑251 (MQ=255)
                                                                              aaGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGt                                                                                                                                                           >  1:230457/1‑251 (MQ=255)
                                                                                              cGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATCTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAgg                                                                                                                                           >  1:48314/1‑251 (MQ=255)
                                                                                                gAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGCATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGAGTTTCAGTTCGCCAGCAGCGTAGGCt                                                                                                                                         >  1:61623/1‑251 (MQ=255)
                                                                                                 aaCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCtt                                                                                                                                        >  1:408288/1‑251 (MQ=255)
                                                                                                  aCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTc                                                                                                                                       >  1:462316/1‑251 (MQ=255)
                                                                                                   cGGTGTAGAAGATCGTTGCAATCACCTCTGCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCa                                                                                                                                      <  1:147218/251‑1 (MQ=255)
                                                                                                         agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                 >  1:234087/1‑250 (MQ=255)
                                                                                                             gaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTGCGCCAGCAGCGTAGGCTTCAGCGTGAATGt                                                                                                                            <  1:214719/251‑1 (MQ=255)
                                                                                                             gaTCGGTGCAATCACCTCTTCCACAGGCGGCATCACGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGt                                                                                                                            <  2:376092/251‑1 (MQ=255)
                                                                                                                          cccTCTTCCACATGCGGCATCTCGATGAGGGGCATGTTATCGCTACGTACAAAACCCGCAGCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGTACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCttt                                                                                                               <  2:29434/250‑1 (MQ=255)
                                                                                                                            ctctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTc                                                                                                              >  1:31435/1‑250 (MQ=255)
                                                                                                                                                 gatgatGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAgcg                                                                                        >  2:87102/1‑251 (MQ=255)
                                                                                                                                                   tgatgTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGACTTCCAGcgcg                                                                                      >  1:403436/1‑251 (MQ=255)
                                                                                                                                                            aTGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACAGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGt                                                                               >  1:405011/1‑250 (MQ=255)
                                                                                                                                                                tATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTg                                                                           >  1:188080/1‑250 (MQ=255)
                                                                                                                                                                 aacgcactttccAAACCCCGCCTCCGGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTGTAATGTTGAATCTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATAGCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGa                                                                          <  2:48314/239‑1 (MQ=255)
                                                                                                                                                                       aCTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGATTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCCGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCAc                                                                   >  2:106699/1‑251 (MQ=255)
                                                                                                                                                                            cAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGccc                                                              >  1:414736/1‑251 (MQ=255)
                                                                                                                                                                             aaaaCCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGccc                                                              >  1:173182/1‑250 (MQ=255)
                                                                                                                                                                               aaCCCCCATCCTGATCGGCGAAGACATCCAACTGCCCGCCACGCGCGCGAACTTCTTCAATGTTGAATGTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATAACGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAgg                                                           <  2:45627/251‑1 (MQ=255)
                                                                                                                                                                                            aTCGGTGAAGAAATACAACGGCCCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTGGTTCGGCGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACATCGCGTGAt                                              <  2:80618/251‑1 (MQ=255)
                                                                                                                                                                                                 cGAAGACATACAACTGACCGCAACGCGCGAGACCTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGTTGTTTg                                         >  1:377458/1‑251 (MQ=255)
                                                                                                                                                                                                   aaGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGTACCGTGTTTCAGTTCGCCAGCATCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCACCGCGCTTGGGTACTGATCGCCCCGGCCCCGGAACAGCGCGTGCTGTGTGTc                                       >  2:120448/1‑251 (MQ=255)
                                                                                                                                                                                                    aGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTc                                       >  1:430992/1‑250 (MQ=255)
                                                                                                                                                                                                     gacagacaACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGAGCGCCACGGCCCAGGAACa                                                       <  2:44676/233‑1 (MQ=255)
                                                                                                                                                                                                     gACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACa                                                       >  1:44676/1‑233 (MQ=255)
                                                                                                                                                                                                             aCTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAAtct                             <  1:2826/251‑1 (MQ=255)
                                                                                                                                                                                                                  ccGCCACGCGCGTGAACTGCTTCAATGTTGAATTTCAGTTATACCAGCAATTCGTTGGCAGGTGCAACAACAATATCCGGCAGATCGGCATCAAGTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGTGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGc                        <  2:357254/251‑1 (MQ=255)
                                                                                                                                                                                                                     ccccgcgcgcgAACTTCTGCAATTTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAAGGCGCCGTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAg                     <  2:92026/248‑1 (MQ=255)
                                                                                                                                                                                                                      cACGCGCTCGTACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTGCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAgc                    <  2:169309/251‑1 (MQ=255)
                                                                                                                                                                                                                      cACGCGAGAGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGGCAGAGATCTCTTTCAACATCACTGCGCCTTCCAGCGCGATTGGGGAATAATCGCCACGGCCCAGGAACACCGCGTGATGTTTGTCAGAGACATCTTATGCCAgc                    >  2:290527/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAgcg                   >  1:329744/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cAATGTTAAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGTCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGAAGTCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTAAATGCGCCTTCCAGCGCTATGGGGAACTGATCGCCACGGCCCAGGAACACCACGTGTAGTTTGTCAGAGACATCTTCTGCCAGAGCTTCAATGGGTTTGTct  >  2:369565/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
TTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCC  >  NC_000913/3911835‑3912316

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