Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 865,295 G→A A357V (GCC→GTC)  moeA ← molybdopterin molybdotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913865,2950GA98.2% 188.7 / NA 57A357V (GCC→GTC) moeAmolybdopterin molybdotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/36);  minor base T (0/1);  total (20/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAA  >  NC_000913/865058‑865532
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
gTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTgcag                                                                                                                                                                                                                                  <  2:354246/251‑1 (MQ=255)
  cAAAATCAAAACAGCGCAGAATGATTTGCCGGTTGTAGCGAAGCATCTCCTGATCGCTGAGTTCCTCCATTACAGGCCTCCGAACAACGCGTTAAACGATTCTACTTCCACCCATGCGCCCACTTCCACATTGCCGCGATCGAGTTCCAGCACGATAAAGCAGTTGACGAGGCTAAAGAAGCTAAAGATATTTGAACCCTGATGTCCGATGGGCGTCACTTCCAGATCGCCATCGACGTTACTTTgcagc                                                                                                                                                                                                                                 >  2:34495/1‑250 (MQ=255)
      aTCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCAcgcc                                                                                                                                                                                                                            <  1:52303/251‑1 (MQ=255)
          aaGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGCGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCgcgc                                                                                                                                                                                                                        >  2:364132/1‑251 (MQ=255)
            gCCGCGCAGAATGATTTGCCGGTGGGAGCGCAGCATCTCCTGAGCGCTGAGTGCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTGCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTGCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCTCGCCGCGCTg                                                                                                                                                                                                                      <  2:401159/251‑1 (MQ=255)
                        gATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACg                                                                                                                                                                                                          <  1:233077/251‑1 (MQ=255)
                             gCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCt                                                                                                                                                                                                      <  1:323032/250‑1 (MQ=255)
                             gCCGGTTGTAGCGCAGCATCTCCTGATCGCGGAGTTCCGCCATTACAGGCCCCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATGCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTACGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGCCCTg                                                                                                                                                                                                     <  2:36193/251‑1 (MQ=255)
                                    gTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGtttt                                                                                                                                                                                              >  1:19627/1‑251 (MQ=255)
                                      aGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGttttt                                                                                                                                                                                             <  2:37583/250‑1 (MQ=255)
                                             caTCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCg                                                                                                                                                                                     <  2:169877/251‑1 (MQ=255)
                                                 cccTGCTCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTACAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGAc                                                                                                                                                                                 <  2:160543/250‑1 (MQ=255)
                                                     gATCGCTGAGTTCAGCCATTACAGGCCTCCGAACAACGCGTTAAACGATTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGTTCGCGTTCCAGCACGATCAAGCAGTTGCCGAGGCTCCAGATGCTAACTATATGTAACACCGGAAGTACGGTGGTCATCACTTCCAGATCACCATCGACGTTGCGTTGCAACACGCCGCGCTGGAAATCCAGCCGTCCTGGCGTTTTATCCAGTCGGGAccgtt                                                                                                                                                                             >  2:151662/1‑247 (MQ=255)
                                                        cGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCtgtg                                                                                                                                                                           <  2:155070/250‑1 (MQ=255)
                                                        cGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCtgtg                                                                                                                                                                           >  1:155070/1‑250 (MQ=255)
                                                            gAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGt                                                                                                                                                                        <  1:165057/249‑1 (MQ=255)
                                                            gAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTAc                                                                                                                                                                      <  1:207351/251‑1 (MQ=255)
                                                             aGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATGCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTTGCGTTTTTTTCAGGCGGGACGCTGTGCGTAcg                                                                                                                                                                     <  1:370858/251‑1 (MQ=255)
                                                                tCCGCCATTACAGGCCTCCGACCAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTCcgcg                                                                                                                                                                   <  2:16354/250‑1 (MQ=255)
                                                                    ccagtACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCTCCCATTCGCCCACGTCCACATTGCCGCGATCGCGTTCCGGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTTCTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTGCAGGCGGGACGCTGTGCGTACGCGCttgc                                                                                                                                                              <  2:129713/247‑4 (MQ=255)
                                                                          cAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGTACGCTGTGCGTGCGCGCTGGCGCGCgg                                                                                                                                                        <  1:124772/251‑1 (MQ=255)
                                                                                 ccGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGcc                                                                                                                                                 >  2:344133/1‑251 (MQ=255)
                                                                                   gAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGc                                                                                                                                               >  1:388859/1‑251 (MQ=255)
                                                                                    aacaacGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCt                                                                                                                                              >  1:121849/1‑251 (MQ=255)
                                                                                       aacGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGc                                                                                                                                           >  1:39031/1‑251 (MQ=255)
                                                                                        acGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCg                                                                                                                                          <  2:314595/251‑1 (MQ=255)
                                                                                         cgcgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGAACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCgg                                                                                                                                         <  1:123297/251‑1 (MQ=255)
                                                                                          gcgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCggg                                                                                                                                        >  2:181026/1‑250 (MQ=255)
                                                                                           cgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGgtg                                                                                                                                       <  1:48978/251‑1 (MQ=255)
                                                                                           cgTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGgtg                                                                                                                                       <  1:364511/251‑1 (MQ=255)
                                                                                              tAAACGGTTCTACTTCCACTCATTCGCCCACTTCCACATTGCCGCGATCGCGATCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTACAGTTCGCCATCGAAGTTGCGTTGCAGCACGCCGCACTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGAGGGAAGCTGTTCGTACGCGCTGGCGCGCGGGCAGGCCggtgtgggcggtt                                                                                                                                    >  2:242684/1‑239 (MQ=255)
                                                                                                   ggTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCt                                                                                                                               <  2:121849/251‑1 (MQ=255)
                                                                                                    gTTCTAATTCCACCTATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGCGCAGGCCGCTGGCGGAGTTATCGcct                                                                                                                              >  2:429454/1‑249 (MQ=255)
                                                                                                         aCTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGGTACCGCTTCGttt                                                                                                                         >  1:1019/1‑251 (MQ=255)
                                                                                                            tCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTg                                                                                                                       <  1:420212/250‑1 (MQ=255)
                                                                                                            tCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGc                                                                                                                      >  1:439432/1‑251 (MQ=255)
                                                                                                               ccccATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAg                                                                                                                   <  1:94260/250‑1 (MQ=255)
                                                                                                               aCCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAg                                                                                                                   <  1:27365/251‑1 (MQ=255)
                                                                                                                                     tgcagcgATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCCGCACAgg                                                                                                             >  2:457180/5‑235 (MQ=255)
                                                                                                                                     tGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAgg                                                                                                             <  1:457180/235‑1 (MQ=255)
                                                                                                                                                ggatCCAGCACGATAAAGCAGTTGACGAGGAGAAAGGAGCTAAATATAAGTGAACCCTGATGTCCGGTGGGCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGt                                                                                  <  2:388859/248‑1 (MQ=255)
                                                                                                                                                           gATAAATCCGTTGCCGAGGCTAAAGGAGCTAAATATATGTGACCCCTGATGTCCGGTGGTCTTCCCTTCCAGTTCGCCAGCGTCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGTGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTGCCGCTTAGTTTTGCCATCAAAGGCTGTACCAGTTGATAGAAGGTCAGCTTCGCGGAAACCg                                                                       <  2:80380/251‑1 (MQ=255)
                                                                                                                                                               aaGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTGCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGCCGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCCTAGGTTTGCCAGCAAAGGCTGTACCCGGTGATAGAAGGTACGCGTCGCTGAAACCGGGtt                                                                   >  1:215243/1‑251 (MQ=255)
                                                                                                                                                                         cGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAg                                                                                                                            <  2:335705/184‑1 (MQ=255)
                                                                                                                                                                         cGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAg                                                                                                                            >  1:335705/1‑184 (MQ=255)
                                                                                                                                                                         cGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAgg                                                         <  2:19627/251‑1 (MQ=255)
                                                                                                                                                                            ggCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGGTTTGCCAGCAAAGGCTGTACCAGGTGATAGAAGGTGAGGGGTGCTGAAACCGGGTTGCCCGGCAgggcg                                                      >  2:227967/1‑248 (MQ=255)
                                                                                                                                                                                     aGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTcc                                                                                                                                                                                                       <  2:283873/97‑1 (MQ=255)
                                                                                                                                                                                     aGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTcc                                                                                                                                                                                                       >  1:283873/1‑97 (MQ=255)
                                                                                                                                                                                               atGTGAACCCTGATCTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTGGTGCTGa                                   <  1:356240/251‑1 (MQ=255)
                                                                                                                                                                                                   gAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGtt                                <  1:213654/250‑1 (MQ=255)
                                                                                                                                                                                                           gATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaa                         <  1:366934/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGGTGGTCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGTCCGCAGAACCAGCTATTGCTGAGTTTACcgaacgcgaac                  <  1:295750/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggtCGTCACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGt               <  1:8707/250‑1 (MQ=255)
                                                                                                                                                                                                                          tcACTTCCAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCCGCTGGTAGAAGGTCAGCGCCGCTGAAACCGGGGTGCCCGGCAGGCCGCAGAACCAGCTCTTGCTGAGTTTACCGAACGCGAAAGGTTTAcac         >  2:349683/1‑249 (MQ=255)
                                                                                                                                                                                                                                ccAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGttt    <  2:112455/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGTTCGCCATCGACGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTaa  <  2:152636/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
GTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGGCCTCCGAACAACGCGTTAAACGGTTCTACTTCCACCCATTCGCCCACTTCCACATTGCCGCGATCGCGTTCCAGCACGATAAAGCAGTTGCCGAGGCTAAAGGAGCTAAATATATGTGAACCCTGATGTCCGGTGGTCGTCACTTCCAGTTCGCCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAA  >  NC_000913/865058‑865532

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: