Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 225,366 C→T 7.9% A99T (GCC→ACC)  PP_0173 ← transcriptional factor‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947225,3660CT 7.9% 178.9 / 3.3 63A99T (GCC→ACC) PP_0173transcriptional factor‑like protein
Reads supporting (aligned to +/- strand):  ref base C (27/31);  new base T (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

CTCAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGGCGC  >  NC_002947/225218‑225506
                                                                                                                                                    |                                                                                                                                            
ctCAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGc                                                                                                                                              >  3:328721/1‑149 (MQ=255)
ctCAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGc                                                                                                                                              >  1:71940/1‑149 (MQ=255)
 tCAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGcc                                                                                                                                             >  1:83229/1‑149 (MQ=255)
     tGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAcc                                                                                                                                         <  1:40262/149‑1 (MQ=255)
     tGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGAc                                                                                                                                          >  4:348678/1‑148 (MQ=255)
       gcgcATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCa                                                                                                                                        <  8:128876/148‑1 (MQ=255)
        cgcATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                      <  7:214984/149‑1 (MQ=255)
           aTCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                      >  7:183617/1‑146 (MQ=255)
           aTCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGt                                                                                                                                      <  8:183617/146‑1 (MQ=255)
              ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGTCGGCCAGTGccgtg                                                                                                                                <  3:301257/149‑4 (MQ=255)
              ccGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAc                                                                                                                                >  7:3859/1‑149 (MQ=255)
               cGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGTAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCaa                                                                                                                                 <  2:159774/147‑1 (MQ=255)
               cGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGTAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCaa                                                                                                                                 >  1:159774/1‑147 (MQ=255)
                gCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAcc                                                                                                                               >  6:28175/1‑148 (MQ=255)
                  aaCTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaa                                                                                                                            >  1:248106/1‑149 (MQ=255)
                          cgcCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGc                                                                                                                    >  3:56777/1‑149 (MQ=255)
                            ccAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGctc                                                                                                                  <  6:60737/149‑1 (MQ=255)
                                    tCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                          <  5:355994/149‑1 (MQ=255)
                                     cGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCACCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCCCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAAAGGAATGCTCTAGACATTg                                                                                                         >  1:341879/1‑149 (MQ=255)
                                       cTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCg                                                                                                       <  6:265261/149‑1 (MQ=255)
                                              gcgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGc                                                                                                >  7:253373/1‑149 (MQ=255)
                                               cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                          <  7:260448/138‑1 (MQ=255)
                                               cgcCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCAtt                                                                                                          >  8:260448/1‑138 (MQ=255)
                                                      aGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAatat                                                                                        >  4:360879/1‑149 (MQ=255)
                                                            aCGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCg                                                                                   >  3:306176/1‑148 (MQ=255)
                                                            aCGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATACTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGa                                                                                  <  2:54750/149‑1 (MQ=255)
                                                                ccTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGc                                                                              >  3:220703/1‑149 (MQ=255)
                                                                 cTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGa                                                                                            <  1:327430/134‑1 (MQ=255)
                                                                 cTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGa                                                                                            >  2:327430/1‑134 (MQ=255)
                                                                 cTCACCCTGCGATGGCTGCAACGCATCCACTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAAACGAATGGAATGCTATTGTCATTGCGTTCCAGCCGGAATATAATCGaaacc                                                                              >  1:301066/1‑146 (MQ=255)
                                                                    aCCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCg                                                                           <  4:220703/148‑1 (MQ=255)
                                                                     cccTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGtt                                                                         <  2:83229/149‑1 (MQ=255)
                                                                          cGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAg                                                                      <  4:56777/147‑1 (MQ=255)
                                                                          cGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAg                                                                      <  5:28175/147‑1 (MQ=255)
                                                                               gCTTCAACGCATCCCGTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACAAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGc                                                               <  2:341879/149‑1 (MQ=255)
                                                                                 tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                     <  6:367937/141‑1 (MQ=255)
                                                                                 tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAgg                                                                     >  5:367937/1‑141 (MQ=255)
                                                                                 tGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTc                                                             >  2:200987/1‑149 (MQ=255)
                                                                                   cAACGAATCCCCTCGCAGATGTGTCAGCAAGTAGAAATTACCGCGCAATTCGCCCAGGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGa                                                           <  2:301066/149‑1 (MQ=255)
                                                                                      cGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaat                                                                                                                           >  5:35451/1‑82 (MQ=255)
                                                                                      cGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCgaat                                                                                                                           <  6:35451/82‑1 (MQ=255)
                                                                                       gCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCtt                                                                                                                >  2:146098/1‑92 (MQ=255)
                                                                                       gCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCtt                                                                                                                <  1:146098/92‑1 (MQ=255)
                                                                                         aTCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGc                                                                                    >  2:224206/1‑118 (MQ=255)
                                                                                         aTCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGc                                                                                    <  1:224206/118‑1 (MQ=255)
                                                                                           ccccTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGt                                                   >  5:316163/1‑149 (MQ=255)
                                                                                             ccTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGttt                                                 <  3:360879/149‑1 (MQ=255)
                                                                                                    gATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAcc                                          <  8:3859/149‑1 (MQ=255)
                                                                                                      tgtgTCAGCAAGTCGAAATTACCGCGCACTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGt                                        <  4:94784/149‑1 (MQ=255)
                                                                                                       gtgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAc                                           <  8:192263/145‑1 (MQ=255)
                                                                                                       gtgtCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGAc                                           >  7:192263/1‑145 (MQ=255)
                                                                                                          tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                     <  5:351816/148‑1 (MQ=255)
                                                                                                          tCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAg                                    >  3:11150/1‑149 (MQ=255)
                                                                                                              cAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCaa                                >  3:198331/1‑149 (MQ=255)
                                                                                                                 gTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGt                             <  3:348678/149‑1 (MQ=255)
                                                                                                                  tCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                     >  2:353283/1‑140 (MQ=255)
                                                                                                                  tCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGaa                                     <  1:353283/140‑1 (MQ=255)
                                                                                                                           ccGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGACTGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGgc                    <  8:253373/148‑1 (MQ=255)
                                                                                                                                 aaTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAc             <  4:11150/149‑1 (MQ=255)
                                                                                                                                  aTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAcc            <  7:140582/149‑1 (MQ=255)
                                                                                                                                  aTTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGAc             >  8:140582/1‑148 (MQ=255)
                                                                                                                                         cACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCgg     >  1:170044/1‑149 (MQ=255)
                                                                                                                                          ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCGGGCGCCAg               >  1:350965/3‑138 (MQ=255)
                                                                                                                                          ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAg               <  2:350965/136‑1 (MQ=255)
                                                                                                                                          ttggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGgc    >  8:252476/3‑149 (MQ=255)
                                                                                                                                            ggCACTGGTCGGCCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGgcgc  >  7:164573/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CTCAATGGCGCATCCCGCAACTGCCGCGCCAGTGCATCGCTGATTCGCGCCACCAGCTCCACGCCCTCACCCTGCGATGGCTGCAACGCATCCCCTCGCAGATGTGTCAGCAAGTCGAAATTACCGCGCAATTCGCCCACGGCACTGGCCGACCAGTGCCAACCGAATGGAATGCTCTTGTCATTGCGTTCCAGCGGGAATATCAGCGAAAGCTCGTTGAGGTCACGCTCGACAGTGCGTTTGCTGACCGTGAAGCCAACGTCACGCAGGCGCCAGACCAGCTCGGCGC  >  NC_002947/225218‑225506

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