Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 336,125 +T coding (210/237 nt) PP_0278 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947336,1241.T98.4% 216.6 / ‑3.2 64coding (209/237 nt)PP_0278hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (33/30);  total (34/30)
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.31e-01

ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335980‑336263
                                                                                                                                                  |                                                                                                                                            
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                          <  1:390445/149‑4 (MQ=255)
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                          >  6:81431/1‑146 (MQ=255)
  ccATCCCTCAAGCTCCTACCCGTAGAGCCCGGCGTCCCCGCCGACGCGGCCTACGAGCATCTCTCGATCCAGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttacct                                                                                                                                        <  8:304691/149‑6 (MQ=255)
       ccTCAATCTCCTACGCGTATAGCCCGCCGTCCCCTTCGACGCGGCCTACGAGCATTTCTCGATCCTGATGCTCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAATTTCCTGGGGGCAGCCGATTACCTGTGCg                                                                                                                                   <  8:377389/149‑1 (MQ=255)
       ccTCAATCTCCTACGCGTACCCCCCCGCCTCCCCGTCCACGCGGCCTACGAGCATGTCTCCATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCTATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCg                                                                                                                                   <  8:383076/149‑1 (MQ=255)
        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCg                                                                                                                                   >  7:386945/1‑148 (MQ=255)
         tCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                  >  4:252173/1‑148 (MQ=255)
           aaTCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACAt                                                                                                                               >  6:78825/1‑149 (MQ=255)
               tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACAt                                                                                                                               >  5:71380/1‑145 (MQ=255)
               tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACAt                                                                                                                               <  6:71380/145‑1 (MQ=255)
               tcCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGc                                                                                                                            >  5:120474/1‑148 (MQ=255)
                      cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGgcgc                                                                                                                    <  1:201804/149‑1 (MQ=255)
                          gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                 <  1:240064/148‑1 (MQ=255)
                                  cGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGata                                                                                                        <  8:386945/149‑1 (MQ=255)
                                     ccccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                      >  5:313044/1‑148 (MQ=255)
                                       ccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                      >  7:196719/1‑146 (MQ=255)
                                        cGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                      >  6:336859/1‑145 (MQ=255)
                                        cGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                      >  6:107382/1‑145 (MQ=255)
                                           cGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGGAGCGCGAGGGCGATATGGAGGCCGACCATAAATTCCTGGGGGCAGCCGATAACCTGTGCGACATGGCCAAGGCGCTGATGAACCAGATCggaaca                                                                                                >  1:111379/1‑144 (MQ=255)
                                             aCGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGcc                                                                                             <  2:270158/149‑1 (MQ=255)
                                               gcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCaa                                                                                           <  2:276278/149‑1 (MQ=255)
                                                   ccTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAggg                                                                                        >  8:97979/1‑148 (MQ=255)
                                                          gCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGgcg                                                                                 >  3:407019/1‑148 (MQ=255)
                                                                   gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                         <  3:252173/147‑1 (MQ=255)
                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                      >  2:197923/1‑149 (MQ=255)
                                                                      ccTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGa                                                                    >  1:212905/1‑149 (MQ=255)
                                                                        tGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                  >  3:311561/1‑149 (MQ=255)
                                                                          aTGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                  <  1:197923/147‑1 (MQ=255)
                                                                                tACATCAAGCACCTGGTGCGCGAGGGCGATATGGAAGACCACCAAAAATTCCTGGCGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCAGATGAACGATATCGAGATCGCCAAGCGTAAGCCGCATTAGGCATTGACGCGGTCaat                                                          >  1:64879/1‑147 (MQ=255)
                                                                                tACATCAAGCACCTGGTGCGCGAGGACGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAtt                                                          >  3:132695/1‑149 (MQ=255)
                                                                                  cATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtg                                                        >  1:154265/1‑149 (MQ=255)
                                                                                   aTCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATGAGGCATTGACGCGGTCATTtgt                                                       >  2:16021/1‑149 (MQ=255)
                                                                                    tCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTtgtg                                                      <  7:97979/149‑1 (MQ=255)
                                                                                      aaGCACCTTGTGAGCGAGGGCGATATGGATTACGAACATAAATTCCTGGGGGCAGCAGATTACCTGTGCGACATGGCCAAGGCGCTGTTGAACGATATCGAGATCGGCAAGGGTATGGCGCATTAGGCATTTACGCTGTCAttt                                                         <  6:404714/144‑1 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAGGGCGCTGATGAACGAGATCGAGATCGCCAAGGGTAAGGCGCTTTAGGCATTGACGGGGTCAGTTGTggg                                                    >  8:144603/1‑149 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                           <  4:364078/126‑1 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                           >  3:364078/1‑126 (MQ=255)
                                                                                      aaGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                     <  5:107382/148‑1 (MQ=255)
                                                                                      aaGCACCTGGTCCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGAGCGACATGGCCAAGGCGCTGATGAACGATATAGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAttt                                                         >  5:404714/1‑144 (MQ=255)
                                                                                       aGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTgg                                                     >  3:100498/1‑147 (MQ=255)
                                                                                        gCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAg                                                  >  5:381607/1‑149 (MQ=255)
                                                                                        gCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAg                                                  >  7:265299/1‑149 (MQ=255)
                                                                                                gTGCGCGAGGACGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                           <  4:132695/149‑1 (MQ=255)
                                                                                                   cgcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGc                                         <  6:120474/148‑1 (MQ=255)
                                                                                                    gcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGt                                       <  5:78825/149‑1 (MQ=255)
                                                                                                       aGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTcgc                                    <  6:313044/149‑1 (MQ=255)
                                                                                                           cGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAt                                >  2:398523/1‑149 (MQ=255)
                                                                                                              tatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                             <  6:381607/149‑1 (MQ=255)
                                                                                                                     gacgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                           >  2:182065/1‑144 (MQ=255)
                                                                                                                     gacgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                           <  1:182065/144‑1 (MQ=255)
                                                                                                                     gacgacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCaaa                      <  8:265299/149‑1 (MQ=255)
                                                                                                                        gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGt                                       >  2:373974/1‑129 (MQ=255)
                                                                                                                        gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGt                                       <  1:373974/129‑1 (MQ=255)
                                                                                                                          ccATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                              >  7:182477/1‑88 (MQ=255)
                                                                                                                          ccATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                              <  8:182477/88‑1 (MQ=255)
                                                                                                                           cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                           >  2:316695/1‑122 (MQ=255)
                                                                                                                           cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCtt                                           <  1:316695/122‑1 (MQ=255)
                                                                                                                           cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcc                <  2:338681/149‑1 (MQ=255)
                                                                                                                           cATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGcc                >  1:338681/1‑149 (MQ=255)
                                                                                                                              aaaTTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGCATGGGGTGCAAATCCGccccc             >  5:81779/1‑148 (MQ=255)
                                                                                                                              aaaTTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCAAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCg             <  2:212905/149‑1 (MQ=255)
                                                                                                                               taTTCGTGGGGGCAGCCTATTACCTGTGCGACATGGCCAAGGGGCTGATGAACGATATCGTGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCgg            <  6:240699/148‑1 (MQ=38)
                                                                                                                                 ttCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCgg            >  5:240699/1‑147 (MQ=255)
                                                                                                                                  tCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGc                         <  4:142954/133‑1 (MQ=255)
                                                                                                                                  tCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGc                         >  3:142954/1‑133 (MQ=255)
                                                                                                                                        gggCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCag   <  5:336859/149‑1 (MQ=255)
                                                                                                                                         ggCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCaga  <  4:407019/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGA  >  NC_002947/335980‑336263

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: