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

TTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCC  >  NC_000913/3911835‑3912316
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
tttATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAATTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCAAGTTATCGCTACTTCCAAAACCCGCATCCTGATCGGCGAAGACATACCACTGAACGCCACGCGCGCGACCTTCTTCACTGTTGAATTTCAGtttt                                                                                                                                                                                                                                         >  2:451909/1‑251 (MQ=255)
     ttcccAACCTTAACCGTTTTTTCCAGGTTACGCGGCGGGTCACCGTCGGGGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGAGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCa                                                                                                                                                                                                                                     <  2:458285/246‑1 (MQ=255)
          aaCCGTACCCGATTGTGCCAGGGTACGCGGCGGGTCAACGTCGGGGCCTTTTATCGGCTCGACATGGTAATCCAGCAGTTGTAGCGGAACGGTTTAGAAAAGCGGTGCAATCCCCTCTTCCACATGCGGCAGCACGATG‑AGGTGCATGTTATCGCTACTTACATATCCCGCATCCTGATCGGCGAAGACATACAACGTTCCGCCACGCGCGCGACCGGCAGCAATGTTGAAGTTCAGTCTGTCCAGCAAtt                                                                                                                                                                                                                               <  2:277630/251‑1 (MQ=255)
             cGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgt                                                                                                                                                                                                                            >  1:152897/1‑251 (MQ=255)
              gTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCTGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCgtt                                                                                                                                                                                                                           >  1:45627/1‑251 (MQ=255)
                    gATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCgg                                                                                                                                                                                                                     >  1:462389/1‑251 (MQ=255)
                           ccAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGcaac                                                                                                                                                                                                               >  1:376092/1‑250 (MQ=255)
                                           gTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGGTG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATAAAACTGACCGCCACGCCCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCATTTCGTTTTTCGGTGCAACAACAATAACCGGCatat                                                                                                                                                                                              >  2:370295/1‑251 (MQ=255)
                                                cGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGc                                                                                                                                                                                          >  1:197691/1‑250 (MQ=255)
                                                    ggTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGCACGGTGTAGAAGATCGGTGCAATCACATCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCa                                                                                                                                                                                      >  1:266006/1‑250 (MQ=255)
                                                            gtaacagagcaataTGGAAAGCCAGTAGTTGCAGCTGAACTGTGTAGAAGATCGGTGCAACCCCCACTTCCACATGCTGCATCGAG‑TGTATGTGCATGTTATCGCTACTTCCAAAACCCGCTACCTGATTGGCGAAGTCAAACAACTTTCCGTCACGCGCTCGAACTTCTTCAATGCTGAATTTCAGCTCTTCCAGCAATTCGGTGTCCGGTGCAACAACAATAACCGGCATATCGGCAGTAATTAGCGcc                                                                                                                                                                             <  2:69814/238‑1 (MQ=255)
                                                                           ggTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTACAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGACATCAATTAGAGCCAGCGGACCGTGTTTc                                                                                                                                                              >  1:261120/1‑251 (MQ=255)
                                                                              aaGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGt                                                                                                                                                           >  1:230457/1‑251 (MQ=255)
                                                                                              cGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATCTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAgg                                                                                                                                           >  1:48314/1‑251 (MQ=255)
                                                                                                gAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGCATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGAGTTTCAGTTCGCCAGCAGCGTAGGCt                                                                                                                                         >  1:61623/1‑251 (MQ=255)
                                                                                                 aaCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCtt                                                                                                                                        >  1:408288/1‑251 (MQ=255)
                                                                                                  aCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTc                                                                                                                                       >  1:462316/1‑251 (MQ=255)
                                                                                                   cGGTGTAGAAGATCGTTGCAATCACCTCTGCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCa                                                                                                                                      <  1:147218/251‑1 (MQ=255)
                                                                                                         agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTg                                                                                                                                 >  1:234087/1‑250 (MQ=255)
                                                                                                             gaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTGCGCCAGCAGCGTAGGCTTCAGCGTGAATGt                                                                                                                            <  1:214719/251‑1 (MQ=255)
                                                                                                             gaTCGGTGCAATCACCTCTTCCACAGGCGGCATCACGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGt                                                                                                                            <  2:376092/251‑1 (MQ=255)
                                                                                                                          cccTCTTCCACATGCGGCATCTCGATG‑AGGGGCATGTTATCGCTACGTACAAAACCCGCAGCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGTACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCttt                                                                                                               <  2:29434/250‑1 (MQ=255)
                                                                                                                            ctctTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTc                                                                                                              >  1:31435/1‑250 (MQ=255)
                                                                                                                                                 gatg‑atGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAgcg                                                                                        >  2:87102/1‑251 (MQ=255)
                                                                                                                                                   tg‑atgTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGAATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGACTTCCAGcgcg                                                                                      >  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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
TTTATTACTCAACCGTAACCGATTTTGCCAGGTTACGCGGCTGGTCAACGTCGGTGCCTTTGATCAGCGCGACATGGTAAGCCAGCAGCTGCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATG‑ATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCAACAACAATAACCGGCATATCGGCATCAATTAGCGCCAGCGGACCGTGTTTCAGTTCGCCAGCAGCGTAGGCTTCAGCGTGAATGTAAGAGATCTCTTTCAACTTCAATGCGCCTTCCAGCGCGATTGGGTACTGATCGCCACGGCCCAGGAACAGCGCGTGATGTTTGTCAGAGAAATCTTCTGCCAGCGCTTCAATGCGTTTGTCC  >  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: ‑