Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 807,510 A→G L324P (CTG→CCG)  bioA ← 7,8‑diaminopelargonic acid synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913807,5100AG92.1% 526.9 / 27.6 204L324P (CTG→CCG) bioA7,8‑diaminopelargonic acid synthetase
Reads supporting (aligned to +/- strand):  ref base A (5/11);  new base G (99/88);  total (105/99)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 209 total aligned reads displayed.

CGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTC  >  NC_000913/807421‑807601
                                                                                         |                                                                                           
cGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGc                                                                                          <  2:3263374/93‑1 (MQ=255)
cGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGc                                                                                          >  1:4185013/1‑93 (MQ=255)
  gCATCACGGGCGGGGTCAAGTGGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTg                                                                                        <  2:8838460/93‑1 (MQ=255)
  gCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTg                                                                                        <  2:4125508/93‑1 (MQ=255)
  gCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTg                                                                                        >  1:770691/1‑93 (MQ=255)
  gCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTg                                                                                        >  2:8211277/1‑93 (MQ=255)
   cATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTgg                                                                                       <  1:7456044/93‑1 (MQ=255)
    aTCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGc                                                                                      >  2:4609869/1‑93 (MQ=255)
     tCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCg                                                                                     >  1:3642295/1‑93 (MQ=255)
     tCACGGGCGGGGGCAAGTTGCTCCCGCCGCTGTACTTCAATATCCGCCACCTGTTGCTCCCAGTCGCCAGATTCGAGGCTCGCCTGGCTcgcg                                                                                     >  1:3486042/1‑93 (MQ=255)
       ccgggcgggGGCAAGGTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGtt                                                                                   <  1:1444951/92‑1 (MQ=255)
       aCGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCgac                                                                                   <  2:138218/93‑3 (MQ=255)
          ggcgggGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTtgc                                                                                <  1:234856/93‑1 (MQ=255)
           gcgggGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTtgct                                                                               >  1:8026874/1‑93 (MQ=255)
           gcgggGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTtgct                                                                               <  2:6856937/93‑1 (MQ=255)
               gggCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGCTTGCTGccc                                                                           >  1:837814/1‑92 (MQ=255)
                 gCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCcgcg                                                                         >  1:1119676/1‑93 (MQ=255)
                   aaGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCa                                                                       >  1:3318153/1‑93 (MQ=255)
                   aaGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCa                                                                       >  1:7579458/1‑93 (MQ=255)
                   aaGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCa                                                                       <  2:6864714/93‑1 (MQ=255)
                    aGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAg                                                                      >  1:2673316/1‑93 (MQ=255)
                    aGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCCGGCTGGCGTTTGCTGCCGCGCAg                                                                      >  2:2559624/1‑93 (MQ=255)
                      ttGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGc                                                                    >  2:904990/1‑93 (MQ=255)
                       tGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGcc                                                                   >  2:3795183/1‑93 (MQ=255)
                       tGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGGTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGcc                                                                   <  1:3876300/93‑1 (MQ=255)
                        gCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCa                                                                  <  1:9051119/93‑1 (MQ=255)
                         cTCGCGCGGCTGTACTTCAATATCCGCCGCCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAg                                                                 <  2:7612868/93‑1 (MQ=255)
                         cTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAg                                                                 <  1:6151170/93‑1 (MQ=255)
                          tCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGc                                                                >  1:2783412/1‑93 (MQ=255)
                          tCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGc                                                                <  1:72066/93‑1 (MQ=255)
                            gcgcAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCgg                                                              <  2:9334858/93‑1 (MQ=255)
                             cgcAGCTGTACTTAAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCTTTTGCTGCCGCGCAGGCCAGCGGa                                                             <  1:2390850/93‑1 (MQ=255)
                              gcAGCTGTACTTCAATATCCGCTGGCTGTTCCTGCCAGTCGCCCGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCggtt                                                            <  2:4822440/93‑3 (MQ=255)
                              gcAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGAt                                                            <  1:2179270/93‑1 (MQ=255)
                              gcAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGAt                                                            <  1:4690341/93‑1 (MQ=255)
                              gcAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGAt                                                            <  2:1113194/93‑1 (MQ=255)
                              ccAGCTGTACTTCAATATCCGCCCCCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGAt                                                            <  2:5529349/92‑1 (MQ=255)
                               cAGCTGTGCTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGAtt                                                           <  1:9134831/93‑1 (MQ=255)
                               cAGCTGTACTTCAATATCGGCCACCTGTTGCTGCCAGTCGCCAGATGCGAGAATCGCCAGGCTGGCGCTTGCTGGCGCGCGGCGCAGGGGAtg                                                           >  1:6952033/1‑92 (MQ=255)
                                aGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTg                                                          >  1:5055644/1‑93 (MQ=255)
                                   tGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGccc                                                       >  2:3949057/1‑93 (MQ=255)
                                   tGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGccc                                                       >  1:6559749/1‑93 (MQ=255)
                                    gTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCa                                                      >  2:7288092/1‑93 (MQ=255)
                                    gTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCa                                                      >  1:3074333/1‑93 (MQ=255)
                                    gTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCa                                                      <  1:4651448/93‑1 (MQ=255)
                                    gTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCGGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCa                                                      <  2:5604509/93‑1 (MQ=255)
                                     tACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCAt                                                     <  1:8224775/93‑1 (MQ=255)
                                      aCTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATa                                                    >  2:7896618/1‑93 (MQ=255)
                                      aCTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATa                                                    >  2:8341972/1‑93 (MQ=255)
                                       cTTTAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaa                                                   <  2:5823598/93‑1 (MQ=255)
                                       cTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaa                                                   >  2:2826855/1‑93 (MQ=255)
                                        ttCAATATCCGCCGCCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaaa                                                  <  1:9277826/93‑1 (MQ=255)
                                        ttCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaaa                                                  >  2:908327/1‑93 (MQ=255)
                                        ttCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGCTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaaa                                                  >  2:267668/1‑93 (MQ=255)
                                         tCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATaaaa                                                 <  2:6967146/93‑1 (MQ=255)
                                          cAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAg                                                >  1:526861/1‑93 (MQ=255)
                                          cAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAg                                                <  2:1868283/93‑1 (MQ=255)
                                          cAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAg                                                <  2:5414389/93‑1 (MQ=255)
                                          cAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGCCAGCGGATTGCCCCTtaaaat                                               >  1:6477275/1‑91 (MQ=255)
                                             tatCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTg                                             >  2:9318500/1‑93 (MQ=255)
                                              atCCGCCACCTGTTGCTGGCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGGCCATAAAAGTTgg                                            >  2:6412012/1‑93 (MQ=255)
                                               tCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGc                                           >  2:1569013/1‑93 (MQ=255)
                                               tCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGc                                           <  1:702886/93‑1 (MQ=255)
                                               tCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGc                                           >  2:5109939/1‑93 (MQ=255)
                                                   ccACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCAt                                       <  2:631181/93‑1 (MQ=255)
                                                   ccACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCAt                                       <  2:7789649/93‑1 (MQ=255)
                                                   ccACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCAt                                       <  1:6295137/93‑1 (MQ=255)
                                                    cACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGGTTGCTGCCGCGCAGGCCAGCGGGTTGCCCATAAAAGTTGGCCCaag                                      >  1:2100702/1‑91 (MQ=255)
                                                      ccTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCCTGCa                                    >  1:2043521/1‑93 (MQ=255)
                                                      ccTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCCTGCa                                    >  2:3724800/1‑93 (MQ=255)
                                                      ccTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCa                                    >  2:6137475/1‑93 (MQ=255)
                                                        tGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATa                                  >  1:8767951/1‑93 (MQ=255)
                                                            gctgcCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGc                              <  1:4468073/93‑1 (MQ=255)
                                                              tgcCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCaa                            <  2:494000/93‑1 (MQ=255)
                                                              tgcCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCaa                            <  2:7478312/93‑1 (MQ=255)
                                                               gcCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAAc                           <  1:9086000/93‑1 (MQ=255)
                                                               gcCAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAAc                           <  1:1892493/93‑1 (MQ=255)
                                                                ccAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAAcc                          >  1:6130884/1‑93 (MQ=255)
                                                                 cAGTCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCg                         >  2:5302468/1‑93 (MQ=255)
                                                                 cAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCg                         <  1:9638205/93‑1 (MQ=255)
                                                                    tCGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCt                      >  1:6762211/1‑93 (MQ=255)
                                                                     cGCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCtt                     <  2:8683733/93‑1 (MQ=255)
                                                                      gCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTc                    >  2:1685371/1‑93 (MQ=255)
                                                                      gCCAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTAGCCCATGCATAAAGCAACCCGCTTc                    >  2:4367766/1‑93 (MQ=255)
                                                                       ccAGATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCa                   >  1:1631806/1‑93 (MQ=255)
                                                                        cAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCAc                  <  2:7040014/93‑1 (MQ=255)
                                                                          gATTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCg                >  2:7807636/1‑93 (MQ=255)
                                                                           aTTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGt               >  2:7164075/1‑93 (MQ=255)
                                                                           aTTCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGt               >  2:4017609/1‑93 (MQ=255)
                                                                            ttCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGtt              <  2:330820/93‑1 (MQ=255)
                                                                             tCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTa             >  1:9906074/1‑93 (MQ=255)
                                                                             tCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTa             >  1:6665390/1‑93 (MQ=255)
                                                                             tCGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAACAACCGGCTTCACCGTTa             >  2:7588287/1‑93 (MQ=255)
                                                                              cGAGAATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTAc            <  2:2406325/93‑1 (MQ=255)
                                                                               gagaATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAAACGGCTTCACCGTTACt           >  1:7274434/1‑93 (MQ=255)
                                                                                agaATCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTg          <  2:6625917/93‑1 (MQ=255)
                                                                                   aTCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATg       <  2:4265842/93‑1 (MQ=255)
                                                                                   aTCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATg       <  2:7975163/93‑1 (MQ=255)
                                                                                    tCGCCGGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATgg      >  1:4561200/1‑93 (MQ=255)
                                                                                        cAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTc  <  2:3133438/93‑1 (MQ=255)
                                                                                         |                                                                                           
CGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCAACCGGCTTCACCGTTACTGATGGTTTC  >  NC_000913/807421‑807601

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