Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,840,286 A→G 10.3% D411D (GAT→GAC) ‡ PP_3390 ← porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,840,2860AG10.3% 180.1 / 9.1 68D411D (GAT→GAC) ‡PP_3390porin
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (3/4);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

AGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAATGCC  >  NC_002947/3840138‑3840421
                                                                                                                                                    |                                                                                                                                       
aGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCa                                                                                                                                         >  2:323382/1‑149 (MQ=255)
  cTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATc                                                                                                                                       >  6:117655/1‑149 (MQ=255)
  cTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATc                                                                                                                                       <  5:248877/149‑1 (MQ=255)
   tCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCg                                                                                                                                      >  5:27739/1‑149 (MQ=255)
   tCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCg                                                                                                                                      >  7:300854/1‑149 (MQ=255)
       gCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGcg                                                                                                                                   >  5:368328/1‑147 (MQ=255)
         ccGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCaaaa                                                                                                                                >  4:175939/1‑149 (MQ=255)
              aaGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGa                                                                                                                           <  1:19737/149‑1 (MQ=255)
              aaGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGa                                                                                                                           <  8:300854/149‑1 (MQ=255)
               aGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAg                                                                                                                          <  3:9677/149‑1 (MQ=255)
                gCGGCCGTCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAAGCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCTGAATTCGAACGTGTACCCCTCGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAgc                                                                                                                         <  4:257621/149‑1 (MQ=255)
                      ggcGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAgcatgcgc                                                                                                                   >  1:381047/1‑143 (MQ=255)
                        cGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTc                                                                                                                 <  1:323382/149‑1 (MQ=255)
                           aaaaGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCa                                                                                                              >  7:25561/1‑149 (MQ=255)
                               gCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgc                                                                                                          >  3:6261/1‑149 (MQ=255)
                                ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         >  5:242546/1‑149 (MQ=255)
                                ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         <  1:360213/149‑1 (MQ=255)
                                ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         <  7:151457/149‑1 (MQ=255)
                                ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgc                                                                                                          <  6:242546/148‑1 (MQ=255)
                                 cACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcgc                                                                                                        <  8:249319/149‑1 (MQ=255)
                                    gCAAAACGCGCAGATCCAAAGCGCAAGCCAAACCGGCTGCGCCCAACGAAGATTGTGTCTCAATTCAGAATTCGAACGAGTACCCCAAGATCACCCGCTTCTCATCCCCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCaa                                                                                                     >  5:353835/1‑149 (MQ=255)
                                     cAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAAc                                                                                                    <  5:117655/149‑1 (MQ=255)
                                      aaaaCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACg                                                                                                   >  5:45986/1‑149 (MQ=255)
                                       aaaCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGc                                                                                                  <  3:345525/149‑1 (MQ=255)
                                         aCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAg                                                                                                <  4:6261/149‑1 (MQ=255)
                                           gcgcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCAt                                                                                                             >  6:206679/1‑134 (MQ=255)
                                           gcgcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCAt                                                                                                             <  5:206679/134‑1 (MQ=255)
                                             gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCgga                                                                                            <  6:368328/149‑3 (MQ=21)
                                             gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCgga                                                                                            <  2:381047/149‑3 (MQ=21)
                                             gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGac                                                                                            <  3:26959/149‑1 (MQ=255)
                                             gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGa                                                                                             >  4:302602/1‑148 (MQ=255)
                                                gATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCggatga                                                                                         <  7:316475/149‑6 (MQ=21)
                                                 aTCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACAt                                                                                        <  3:302602/149‑1 (MQ=255)
                                                     aaaGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGCCACATTcc                                                                                     >  1:347236/1‑148 (MQ=255)
                                                             aGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccc                                                                            <  5:308344/149‑1 (MQ=255)
                                                              gCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGcccc                                                                           >  4:382533/1‑149 (MQ=255)
                                                               cGAAACCGGCTGCGCCCAACGAAGATTGTGGCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccccc                                                                          <  8:306939/149‑1 (MQ=255)
                                                                gAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccccc                                                                          >  1:65440/1‑148 (MQ=255)
                                                                    ccGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTg                                                                     >  5:407777/1‑149 (MQ=255)
                                                                            gcCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGATCGCCCCCGACTGCACcaca                                                              >  6:50023/1‑148 (MQ=255)
                                                                                aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGca                                                                   >  5:146430/1‑139 (MQ=255)
                                                                                aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGca                                                                   <  6:146430/139‑1 (MQ=255)
                                                                                aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGt                                                         <  1:272329/149‑1 (MQ=255)
                                                                                 aCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTg                                                        <  6:45986/149‑1 (MQ=255)
                                                                                    aaGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGata                                                     >  7:344623/1‑149 (MQ=255)
                                                                                        ttGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCg                                                 <  6:27739/149‑1 (MQ=255)
                                                                                           tgtCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTa                                               <  2:65440/148‑1 (MQ=255)
                                                                                               tGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTg                                                                     >  8:363246/1‑122 (MQ=255)
                                                                                                gAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGc                                         >  2:70190/1‑149 (MQ=255)
                                                                                                   ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTcc                                                                                     <  8:154620/102‑1 (MQ=255)
                                                                                                   ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTcc                                                                                     >  7:154620/1‑102 (MQ=255)
                                                                                                   ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTcc                                      >  7:123213/1‑149 (MQ=255)
                                                                                                   ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTAGCATTGCGCCAACGAAGCGACACATTACTGAGCGCACCCGAATGAACCACATAGGTGATATCCGTAGTCCGCTcc                                      >  1:14580/1‑149 (MQ=255)
                                                                                                             gAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGa                                                                                  >  6:366645/1‑95 (MQ=255)
                                                                                                             gAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGa                                                                                  <  5:366645/95‑1 (MQ=255)
                                                                                                              aaCGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGccc                           <  6:179000/149‑1 (MQ=255)
                                                                                                                  tgtACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCctct                       <  3:175939/149‑1 (MQ=255)
                                                                                                                     aCCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTACACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCg                    >  2:357900/1‑149 (MQ=255)
                                                                                                                       cccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTc                  >  7:270807/1‑149 (MQ=255)
                                                                                                                       cccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTACACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTc                  <  1:357900/149‑1 (MQ=255)
                                                                                                                        ccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCg                 >  2:403468/1‑149 (MQ=255)
                                                                                                                        ccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCg                 >  1:380871/1‑149 (MQ=255)
                                                                                                                         cACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGc                <  2:380871/149‑1 (MQ=255)
                                                                                                                                 tctggTGCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCTCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGc         <  7:256880/145‑1 (MQ=255)
                                                                                                                                 cccGGTTCTCATCAGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGc         >  5:173543/1‑148 (MQ=255)
                                                                                                                                  aaggTTCTCAACCGCTTCATCGGCAAAATTAGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTACTGAGCGCCCCCGACTGCACCAAATAGGTGAAATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGcc        >  2:236031/3‑148 (MQ=255)
                                                                                                                                   cGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCaa      >  1:160203/1‑149 (MQ=255)
                                                                                                                                     gTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCgac         >  8:256880/1‑142 (MQ=255)
                                                                                                                                        ctcATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAAtgcc  <  6:407777/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
AGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAATGCC  >  NC_002947/3840138‑3840421

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: ‑