Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,622,932 A→C 10.9% M269R (ATG→AGG) ‡ PP_4941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,622,9320AC10.9% 177.0 / 11.3 63M269R (ATG→AGG) ‡PP_4941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base C (4/3);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGG  >  NC_002947/5622787‑5623078
                                                                                                                                                 |                                                                                                                                                  
ggCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATgg                                                                                                                                                 <  6:49336/149‑1 (MQ=255)
   tCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcag                                                                                                                                              >  2:129960/1‑149 (MQ=255)
     gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcag                                                                                                                                              <  1:129960/147‑1 (MQ=255)
     gCAGGTCGGCCAGCAGGGTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcaggc                                                                                                                                            <  7:242107/149‑1 (MQ=255)
                     gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAg                                                                                                                             >  5:274644/1‑148 (MQ=255)
                        gtgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAg                                                                                                                          <  5:18756/148‑1 (MQ=255)
                        gtgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGc                                                                                                                         <  4:212691/149‑1 (MQ=255)
                         tgCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGTAGGCAAAAGTAGcc                                                                                                                        <  8:39095/149‑1 (MQ=255)
                           ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                               >  1:28386/1‑124 (MQ=255)
                           ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                               <  2:28386/124‑1 (MQ=255)
                           ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAg                                                                                                                          <  6:154070/145‑1 (MQ=255)
                           ccAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAg                                                                                                                          >  5:154070/1‑145 (MQ=255)
                            cAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATGCAGGCAGca                                                                                                                                        <  7:123232/130‑2 (MQ=255)
                              ggtgggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCt                                                                                                                   <  5:237124/149‑1 (MQ=255)
                                tgggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGt                                                                                                                 <  6:274644/149‑1 (MQ=255)
                                   gtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATg                                                                                                              <  6:327460/149‑1 (MQ=255)
                                      cgcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCaa                                                                                                           <  1:185182/149‑1 (MQ=255)
                                             gcgGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACt                                                                                                    >  3:106108/1‑149 (MQ=255)
                                              cgGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTc                                                                                                   >  7:239440/1‑149 (MQ=255)
                                                   gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTTCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              <  8:239440/149‑1 (MQ=255)
                                                   gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              >  2:259371/1‑149 (MQ=255)
                                                     ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                            >  2:173498/1‑149 (MQ=255)
                                                         tAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCgg                                                                                        >  6:185579/1‑149 (MQ=255)
                                                           acgggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGc                                                                                       <  4:106108/148‑1 (MQ=255)
                                                             gggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATg                                                                                                                                                  >  7:282534/1‑87 (MQ=255)
                                                             gggaacgggTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATg                                                                                                                                                  <  8:282534/87‑1 (MQ=255)
                                                                     gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGc                                                                                                             >  1:232258/1‑116 (MQ=255)
                                                                     gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGc                                                                                                             <  2:232258/116‑1 (MQ=255)
                                                                     gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCa                                                                            >  6:39006/1‑149 (MQ=255)
                                                                     gTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCa                                                                            >  4:194112/1‑149 (MQ=255)
                                                                         aaCTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCacgac                                                                        >  2:154389/1‑149 (MQ=255)
                                                                              gCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATg                                                                    <  5:39006/148‑1 (MQ=255)
                                                                              gCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGc                                                                   >  5:176791/1‑149 (MQ=255)
                                                                                 gAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGt                                                                <  1:259371/149‑1 (MQ=255)
                                                                                     acGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTcc                                                              >  1:17391/1‑147 (MQ=255)
                                                                                     acGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCa                                                            <  6:176791/149‑1 (MQ=255)
                                                                                        cgGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGGCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGc                                                         <  2:216387/149‑1 (MQ=255)
                                                                                               cGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTAGGGTGCTGCCAGCATGGCAGGCAGCTGCATGGAAAAGTAGCCCAGCTGTATTCAACGCGACTCCCGCACGCCGGCCCAACGCGGGACCACGATGCTGTCCCACGAACGCAgg                                                  >  2:302294/1‑149 (MQ=255)
                                                                                                  ccAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgcgc                                               <  3:194112/149‑1 (MQ=255)
                                                                                                       tCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCtt                                          >  4:64727/1‑149 (MQ=255)
                                                                                                       tCGGCTTCACGTGAGGCGGCCGGCTGGACGGGTGCTGCCAGAATGGCAGGAAGCGGCAGCGAAAAATAGCCAAGCAGTATGCAAACCGACTACGGCAAGCCGGCCAAACCCGGCACGCCGAAGCATTCCCAGGCACGCCGGCccccaa                                           >  6:44648/1‑146 (MQ=255)
                                                                                                         ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTc                                         <  3:209865/148‑1 (MQ=255)
                                                                                                              cACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgc                                   >  7:98189/1‑149 (MQ=255)
                                                                                                                cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcgc                                 <  8:98189/149‑1 (MQ=255)
                                                                                                                 gTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGcgcgcg                                <  3:64727/149‑1 (MQ=255)
                                                                                                                  tGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcgg                               >  3:77583/1‑149 (MQ=255)
                                                                                                                       cggccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                          <  4:77583/149‑1 (MQ=255)
                                                                                                                               tGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGc                  <  1:154389/149‑1 (MQ=255)
                                                                                                                                  tCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAg               <  2:51098/149‑1 (MQ=255)
                                                                                                                                   ctgcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                       <  4:247142/119‑1 (MQ=255)
                                                                                                                                   ctgcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAg                                       >  3:247142/6‑124 (MQ=255)
                                                                                                                                    gggTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc             <  1:173498/149‑1 (MQ=255)
                                                                                                                                     gcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc             <  4:166902/145‑1 (MQ=255)
                                                                                                                                     gcagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc             >  3:166902/4‑148 (MQ=255)
                                                                                                                                      cagctgcCTGCCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGcc            >  2:123578/3‑148 (MQ=255)
                                                                                                                                        gctgcCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCAt         >  3:8087/1‑149 (MQ=255)
                                                                                                                                         ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCCGGCACGCAGGCGCCCCGTCAGGCGCGCGGCCTCGGCCTTTGCCCGGCCCATc        >  7:258735/1‑149 (MQ=255)
                                                                                                                                           gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGc                                                     >  6:95834/1‑102 (MQ=255)
                                                                                                                                             ctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgct    <  5:185579/149‑1 (MQ=255)
                                                                                                                                              tgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgctg   >  7:219880/1‑149 (MQ=255)
                                                                                                                                               gcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg  >  5:119276/1‑149 (MQ=255)
                                                                                                                                                cATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCAt         >  7:143910/1‑141 (MQ=255)
                                                                                                                                                cATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCAt         <  8:143910/141‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGATATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGG  >  NC_002947/5622787‑5623078

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