Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 808,398 A→G 6.0% L287P (CTC→CCC)  bioA ← 7,8‑diaminopelargonic acid synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913808,3980AG 6.0% 203.0 / 3.1 67L287P (CTC→CCC) bioA7,8‑diaminopelargonic acid synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTC  >  NC_000913/808251‑808546
                                                                                                                                                   |                                                                                                                                                    
aCCTGTTGCTGCCAGTCGCCACATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAg                                                                                                                                                     >  6:215026/1‑149 (MQ=255)
     ttGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCTTAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAgtgt                                                                                                                                                  <  4:296301/147‑1 (MQ=255)
       gctgcCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGc                                                                                                                                               <  2:55596/148‑1 (MQ=255)
       gctgcCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCg                                                                                                                                              >  1:157021/1‑149 (MQ=255)
               tCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGt                                                                                                                                      >  5:160261/1‑149 (MQ=255)
                 gCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCa                                                                                                                                    >  1:247265/1‑149 (MQ=255)
                   cAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCAtt                                                                                                                                  >  3:284433/1‑149 (MQ=255)
                   cAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCAtt                                                                                                                                  <  2:157287/149‑1 (MQ=255)
                   cAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCAtt                                                                                                                                  <  8:15103/149‑1 (MQ=255)
                    aGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGc                                                                                                                                               >  7:223536/1‑135 (MQ=255)
                    aGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGc                                                                                                                                               <  8:223536/135‑1 (MQ=255)
                       ttCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGACCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTgc                                                                                                                              >  6:115469/1‑149 (MQ=255)
                          gagaATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTgccgc                                                                                                                           >  7:65217/1‑149 (MQ=255)
                          gagaATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTgccgc                                                                                                                           >  2:145315/1‑149 (MQ=255)
                             aaTCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCgg                                                                                                                        <  2:247265/149‑1 (MQ=255)
                              aTCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGt                                                                                                                       <  7:244592/149‑1 (MQ=255)
                               tCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGt                                                                                                                       >  8:248208/1‑148 (MQ=255)
                                  ccAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAg                                                                                                                   <  4:66350/149‑1 (MQ=255)
                                             tttGCTGCCGCGCAGGCCAGCGGATTGACCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTAGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAg                                                                                                        >  1:101651/1‑149 (MQ=255)
                                                gctgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGca                                                                                                     >  2:195375/1‑149 (MQ=255)
                                                  tgcCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGcac                                                                                                    >  6:203481/1‑148 (MQ=255)
                                                     cgcgCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACaaaa                                                                                                <  7:248208/149‑1 (MQ=255)
                                                      gcgcAGGCCAGCGGATGGCCCATAACAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGGTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAAt                                                                                               <  6:142271/149‑1 (MQ=255)
                                                               aGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGcac                                                                                                    <  2:11447/135‑1 (MQ=255)
                                                               aGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGcac                                                                                                    >  1:11447/1‑135 (MQ=255)
                                                                    aTTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGAttt                                                                                 <  4:195253/149‑1 (MQ=255)
                                                                     ttGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGa                                                                                                         <  1:190786/124‑1 (MQ=255)
                                                                     ttGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGa                                                                                                         >  2:190786/1‑124 (MQ=255)
                                                                        cccATAAAAGTTGGCCCAAGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCAAAGTGCCGCCGGTTAAGGCATTACCGAGGCACAAAATGTCCGGCGCGATTTCTg                                                                              >  5:359524/1‑148 (MQ=255)
                                                                         ccATAACAGTTGTCCCCTGAATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGAAACCTCGCGCGTTGTTATTTGGGCTGAAATGGTCATTGTGCCTCCGGCTAAGGCTTTACAGAGGAACAAAATGTCCGGCGCGATTTCTGCa                                                                            <  1:354033/149‑1 (MQ=255)
                                                                         ccATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCGCCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGGTAAGGCTTTACCGAGGCACAAACTGTCCGGCGCGATGTCTGCa                                                                            >  2:354033/1‑149 (MQ=255)
                                                                         ccATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCa                                                                            <  5:53480/149‑1 (MQ=255)
                                                                          cATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCAt                                                                           >  5:112512/1‑149 (MQ=255)
                                                                           aTAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATg                                                                          >  1:247132/1‑149 (MQ=255)
                                                                           aTAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATg                                                                          >  6:240940/1‑149 (MQ=255)
                                                                              aaaGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTcc                                                                                           <  3:91867/129‑1 (MQ=255)
                                                                              aaaGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTcc                                                                                           >  4:91867/1‑129 (MQ=255)
                                                                              aaaGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGGCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATTTCCGGCGCGATTTCTGCATGtt                                                                        <  5:215026/148‑1 (MQ=255)
                                                                                  ttGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAg                                                                   <  1:145315/149‑1 (MQ=255)
                                                                                    ggCCCATGCATAAAGCAACCGGCTTCACCGTTTCTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                  >  4:46722/1‑148 (MQ=255)
                                                                                     gCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGc                                                                 <  8:65217/148‑1 (MQ=255)
                                                                                       ccATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCaaa                                                              >  7:192277/1‑149 (MQ=255)
                                                                                             aTAAAGCAACCGGCTTCACCGTTTCTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGttt                                                        <  3:46722/149‑1 (MQ=255)
                                                                                                 aGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCg                                                    <  2:101651/149‑1 (MQ=255)
                                                                                                  gCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCgg                                                   >  7:309586/1‑149 (MQ=255)
                                                                                                    aaCCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTa                                                 >  8:235833/1‑149 (MQ=255)
                                                                                                    aaCCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTa                                                 <  5:240940/149‑1 (MQ=255)
                                                                                                        ggCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGAc                                             <  4:284433/149‑1 (MQ=255)
                                                                                                            tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGGGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                  >  2:274008/1‑124 (MQ=255)
                                                                                                            tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGGGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                  <  1:274007/124‑1 (MQ=255)
                                                                                                            tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGGGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                  >  2:274007/1‑124 (MQ=255)
                                                                                                            tCACCGTTACTGATGGTTTCTGCAACCTCGCGCGGGGTGGGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAgg                                                                  <  1:274008/124‑1 (MQ=255)
                                                                                                              aCCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCAt                                                                           >  2:84364/1‑113 (MQ=255)
                                                                                                              aCCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCAt                                                                           <  1:84364/113‑1 (MQ=255)
                                                                                                              aCCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATc                                       >  3:250479/1‑149 (MQ=255)
                                                                                                                      tGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGc                                >  2:96108/1‑148 (MQ=255)
                                                                                                                      tGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAAAAGTTTCCCGGTACGCCCAAATCCAGTgtc                                >  2:145387/1‑146 (MQ=255)
                                                                                                                      agATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGATTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGc                                >  3:220668/2‑148 (MQ=255)
                                                                                                                           gTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGAtc                            <  3:294707/147‑1 (MQ=255)
                                                                                                                           gTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGAtc                            >  4:294707/1‑147 (MQ=255)
                                                                                                                                  cAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCtcgtcg                      <  2:247132/146‑1 (MQ=255)
                                                                                                                                    aCCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAAt                 >  3:73927/1‑149 (MQ=255)
                                                                                                                                         gcgcgTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCTCATCCAGTGGCGATCTCGTCTGCAATCAGc             >  8:354471/1‑148 (MQ=255)
                                                                                                                                         gcgcgTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGc             >  6:277164/1‑148 (MQ=255)
                                                                                                                                         gcgcgTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGc             <  5:203481/148‑1 (MQ=255)
                                                                                                                                             gtggtgAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGAt        <  1:96108/149‑1 (MQ=255)
                                                                                                                                             gtggtgAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGAt        <  5:115469/149‑1 (MQ=255)
                                                                                                                                                   aGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTc  <  7:78467/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTCTGCAACCTCGCGCGTGGTGAGTGTGGCGGAAAGGGTCATTGTGCCGCCGGTTAAGGCTTTACCGAGGCACAAAATGTCCGGCGCGATTTCTGCATGTTCACAGGCAAACAGTTTCCCGGTACGACCAAATCCAGTGGCGATCTCGTCGGCAATCAGCAAGATACCTTC  >  NC_000913/808251‑808546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑