Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 465,005 G→A V181V (GTC→GTT queC ← 7‑cyano‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913465,0050GA95.1% 125.8 / NA 41V181V (GTC→GTTqueC7‑cyano‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/16);  minor base C (2/0);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGGACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACAAA  >  NC_000913/464759‑465245
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
gAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACt                                                                                                                                                                                                                                                >  2:135489/1‑250 (MQ=255)
          tGCGCGAAAAGCTGGAGCAGATGGTAAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGtt                                                                                                                                                                                                                                     <  2:347394/250‑1 (MQ=255)
               gAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGt                                                                                                                                                                                                                               >  1:119996/1‑251 (MQ=255)
               gAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTAGAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGt                                                                                                                                                                                                                               <  2:80924/251‑1 (MQ=255)
                           cAGATGGGTAAGTAAACGTTTTGTGGTGCCGGATGTTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCa                                                                                                                                                                                                                   <  2:119996/251‑1 (MQ=255)
                               tGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCCACGTTTCGTTACGAACTAAATCAAGTTTGCCGTAATAATCTGCCAgcgc                                                                                                                                                                                                               >  1:299893/1‑251 (MQ=255)
                                     attaagctgttggttgtgCCGGATACTCAAGCCTATACCGGCGACACGCGGAAGAATTTCCTCAACCCGGTTTGCTGCTTCATCGCTGCCATCCCCGTCGGTTTATCGGCCAGATAATGATTCAAACCGGTGGCGCGTAAATTGCATGCCGCACAATGACCGCAACCGTCGTCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAGCTAAATCCAGTTTGCCGTAAGACTCTGCCAACGCCCAGGt                                                                                                                                                                                                         <  2:25844/236‑1 (MQ=255)
                                                    gtgCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAAt                                                                                                                                                                                          <  1:483419/251‑1 (MQ=255)
                                                               tctaGCCGCATCCGGCGACACACGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTAGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGAAACCGTCGCCTTTCAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATACACATCAgcg                                                                                                                                                                               >  2:252294/4‑251 (MQ=255)
                                                                cAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTGGCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAgc                                                                                                                                                                                <  1:52164/249‑1 (MQ=255)
                                                                       cATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTc                                                                                                                                                                        >  1:480206/1‑250 (MQ=255)
                                                                        aTCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGGTTCa                                                                                                                                                                       >  2:431791/1‑250 (MQ=255)
                                                                         tCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCCGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCaa                                                                                                                                                                      >  1:74752/1‑250 (MQ=255)
                                                                               cGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGa                                                                                                                                                                 >  1:360275/1‑250 (MQ=255)
                                                                                 acacCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAatt                                                                                                                                                              >  2:252270/1‑250 (MQ=255)
                                                                                                aCCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCACCCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGCACTAAATCCAGt                                                                                                                                                                                                                                      >  1:198516/1‑164 (MQ=255)
                                                                                                aCCTCAACCCGGGTTTCTGCGTCATCGCTGCCATCTCTGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGTCGCGTAAATTACATGCCGCACAATTTCCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGGTTCGGTACGAACTAAATCTAGt                                                                                                                                                                                                                                      <  2:198516/164‑1 (MQ=255)
                                                                                                           gTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAc                                                                                                                                       <  2:139351/248‑1 (MQ=255)
                                                                                                           gTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAc                                                                                                                                       >  1:139351/1‑248 (MQ=255)
                                                                                                            ttttCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGc                                                                                                                                    >  1:258353/1‑250 (MQ=255)
                                                                                                                 tGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGtt                                                                                                                              >  2:540352/1‑251 (MQ=255)
                                                                                                                                  cACCGTCGGTTTATCGGCCAGATAATGGTTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACt                                                                                                             >  2:9719/1‑251 (MQ=255)
                                                                                                                                   cccttCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCATGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCGTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTc                                                                                                            <  1:88035/247‑1 (MQ=255)
                                                                                                                                      gTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATc                                                                                                         >  2:180521/1‑251 (MQ=255)
                                                                                                                                           tttatcacccAGATAATGATTCCAACCGTTGGCGCGTAAATTCCATGCCGCCAAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGATACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCAACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGAAGGCATGGTTTAGTGCTTTAACAAACTCATCGCgg                                                                                                     >  2:62175/9‑250 (MQ=255)
                                                                                                                                                cGGCCAGATCATGATTCCAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATcc                                                                                               >  2:489658/1‑251 (MQ=255)
                                                                                                                                                        tAATGATGCAAACCGGTGGCGCGTAAATTACATGCCGCACAATGACCGCAAACGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCACATTTGCCGTAATAATATGCCAACGCCCAGGTTTCCGCTTTATCAATCCCAATCAGCGGCGTTTCAAAAAGAATATCTTTCGCCATGCACAAACTGACGGCATGGGTTAATGCTTTCACAAACTCCTCGCGGCatacagggaagcc                                                                                       >  2:292121/1‑239 (MQ=255)
                                                                                                                                                            gATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCATGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATAGCTATCGCCATGCCCAAACGGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGgag                                                                                   <  2:230372/251‑1 (MQ=255)
                                                                                                                                                                 aaaCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCAt                                                                                                          >  1:209137/1‑223 (MQ=255)
                                                                                                                                                                 aaaCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGGTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCAt                                                                                                          <  2:209137/223‑1 (MQ=255)
                                                                                                                                                                   aCCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCCGGTCAACGTTTCGTTACGCACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCACTCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCg                                                                            >  1:536002/1‑251 (MQ=255)
                                                                                                                                                                        tGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCg                                                                       >  2:137475/1‑251 (MQ=255)
                                                                                                                                                                        tGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAAGCCGGGTAGCCGGAGAAATCCGTTTc                                                                        <  2:518382/250‑1 (MQ=255)
                                                                                                                                                                                 aaTTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTCTCGCAGACGCCa                                                              <  1:444633/251‑1 (MQ=255)
                                                                                                                                                                                          ccGCACAAAGACCGCAACCGTCGACTTTACAGCCGTTATAGCAGGTCAACGGTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCTCCCAGGTTTCCGCTATATCAATCAACATCAGCGGCGTTTCAAAACGCATATCTTTAGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACACACTCATCGCGGCAATCCAGGTAACCGGAGAAATCCTTTTCGCAGACGCCAATCATTACg                                                     >  2:373360/1‑251 (MQ=255)
                                                                                                                                                                                             cacaATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCt                                                  >  1:443566/1‑251 (MQ=255)
                                                                                                                                                                                                 atgtccGCAACCGGTGCCTTTAAAGGCGTTATAGCAGGTCAACTGTTCGTTACGAACTAAATCCAGTTTGCCGTAATAAGCTGCCAGCGCCCAGGTTTCCGCTTTAGGAATCCACTTCAGCGGCGTTAGAAAACGAAGAGCATTCGCCATGTCCAAACTAACGGCATCGTTGAGTGCTTTCACAACCTCATCGCGGCAATCCGTGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTg                                              <  1:550736/247‑1 (MQ=255)
                                                                                                                                                                                                  ggACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGGCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTAGCAATCCACATCAGCGGCGTTTCAAAACGAATAGCTTTCGCCATGCCCAAACTGACGGCAGGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGAGGCCAGTAAGTACGGCTTCTGc                                             <  2:536002/250‑1 (MQ=255)
                                                                                                                                                                                                                      aaaGCCGTTATAGCAGGTCAACGTTTCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCatat                          >  1:287289/1‑250 (MQ=255)
                                                                                                                                                                                                                                 agccggTCAACGCTTCGTTACGAACTAAATCCAGTTTGACGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACATGATACGCCAATATTgccgcccac              >  2:67132/5‑248 (MQ=255)
                                                                                                                                                                                                                                              ttCGTTACGAACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACaaa  <  1:540352/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
GAAGGGGAATTGCGCGAAAAGCTGGAGCAGATGGTTAAGTAAACGTTTTGTGGTGCCGGATGCTCAAGCCGCATCCGGCGACACCCGGAATAATTACCTCAACCCGGTTTTCTGCTTCATCGCTGCCATCACCGTCGGTTTATCGGCCAGATAATGATTCAAACCGTTGGCGCGTAAATTACATGCCGCACAATGACCGCAACCGTCGCCTTTAAAGCCGTTATAGCAGGTCAACGTTTCGTTACGGACTAAATCCAGTTTGCCGTAATAATCTGCCAGCGCCCAGGTTTCCGCTTTATCAATCCACATCAGCGGCGTTTCAAAACGAATATCTTTCGCCATGCCCAAACTGACGGCATGGTTTAGTGCTTTCACAAACTCATCGCGGCAATCCGGGTAGCCGGAGAAATCCGTTTCGCAGACGCCAGTAATTACGGCTTCTGCTTTTACCTGATACGCATATATTGCCGCCAGCGTCAGGAACAAA  >  NC_000913/464759‑465245

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