Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,423,090 G→C 9.2% intergenic (‑143/‑23) PP_5645 ← / → PP_3914 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,423,0900GC 9.2% 217.8 / 12.9 76intergenic (‑143/‑23)PP_5645/PP_3914hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (37/30);  new base C (2/5);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

ACAGCCCTCCTGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTCGGCTCCACCGT  >  NC_002947/4422948‑4423231
                                                                                                                                              |                                                                                                                                             
aCAGCCCTCCTGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATgg                                                                                                                                         <  3:159690/149‑1 (MQ=255)
  aGCCCTCCTGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATggcg                                                                                                                                       >  4:391215/1‑149 (MQ=255)
     cctcctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATggcggc                                                                                                                                     >  7:55120/1‑148 (MQ=255)
      ctcctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcag                                                                                                                                   <  8:35332/149‑3 (MQ=255)
      ctcctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcag                                                                                                                                   <  4:275750/149‑3 (MQ=255)
         ctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAgg                                                                                                                                <  2:294164/149‑1 (MQ=255)
         ctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcagt                                                                                                                                  <  7:101091/147‑4 (MQ=255)
         ctGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATCGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCaa                                                                                                                                  <  8:359394/147‑1 (MQ=255)
          tGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGa                                                                                                                               >  8:81728/1‑149 (MQ=255)
           gCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGAt                                                                                                                              >  3:50370/1‑149 (MQ=255)
             aaCCGACTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTc                                                                                                                            >  6:361020/1‑149 (MQ=255)
                    tGTTGGATGCCTGTAAAAAAATAGCAAAAGGACGTCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCaag                                                                                                                      >  8:96542/1‑146 (MQ=255)
                           tGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATc                                                                                                              <  3:391215/149‑1 (MQ=255)
                            gCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCa                                                                                                             >  3:289529/1‑149 (MQ=255)
                              cTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTACAATTATTAGGGACTCTCACTAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGAAGCTATATGGCGGCCAAGGATTCCCCATGTCTGAATAACg                                                                                                           >  7:213387/1‑149 (MQ=255)
                                 tAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGTCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCct                                                                                                        >  8:156203/1‑149 (MQ=255)
                                 tAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCct                                                                                                        >  1:95293/1‑149 (MQ=255)
                                   aaaaaaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGCCTCTCCCGAATAATGGCCGGACTTGTTGTAAACCCGCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACCCctgc                                                                                                      >  7:363651/1‑149 (MQ=39)
                                    aaaaaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCctgct                                                                                                     >  8:234228/1‑149 (MQ=255)
                                    aaaaaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCctgct                                                                                                     >  3:165152/1‑149 (MQ=255)
                                       aaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCctg                                                                                                       <  4:187986/144‑1 (MQ=255)
                                       aaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCctg                                                                                                       >  3:187986/1‑144 (MQ=255)
                                       aaTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCAAGTCTGAATCACGCCTGCTGAc                                                                                                  >  6:150540/1‑149 (MQ=255)
                                        aTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGAcc                                                                                                 <  2:95293/149‑1 (MQ=255)
                                        aTAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGAcc                                                                                                 >  2:68944/1‑149 (MQ=255)
                                            ccAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGc                                                                                              >  3:284240/1‑148 (MQ=255)
                                             cAAATGACATCATCATTAAAGTACCTCAAATTATTAGTGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCga                                                                                            <  7:357036/149‑1 (MQ=38)
                                             cAAAGGCCATGATGATTCAAGTCCCGTCAATTATTATGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCTGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCga                                                                                            <  6:215878/149‑1 (MQ=38)
                                               aaGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCgaga                                                                                          <  5:361020/149‑1 (MQ=255)
                                                   ccATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTg                                                                                      <  5:83087/149‑1 (MQ=255)
                                                    cATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTgg                                                                                     <  6:230529/149‑1 (MQ=255)
                                                    cATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcagt                                                                                                                                  <  4:424576/104‑4 (MQ=255)
                                                    cATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcagt                                                                                                                                  >  3:424576/1‑101 (MQ=255)
                                                     aTCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGt                                                                                    <  4:284240/149‑1 (MQ=255)
                                                          cATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTc                                                                               <  8:100750/149‑1 (MQ=255)
                                                            ttCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGTCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTcga                                                                             <  7:156203/149‑1 (MQ=255)
                                                             tCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTcgac                                                                            >  6:395993/1‑149 (MQ=255)
                                                               aaGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTcgacga                                                                          >  6:205066/1‑149 (MQ=255)
                                                               aaGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCAGGTCTGAATCCCTCCTGCTGACCGGCGAGCCCTGGTCTGTcgacga                                                                          >  8:400691/1‑149 (MQ=255)
                                                               aaGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCATCTCACCTCTATTCACTGACGCGATATTGCGGCAAAGGATTCCCCATGTCTGAATCACGCCTTCTGACCGTCGAGACATTGTATGTcgacga                                                                          >  8:71764/1‑149 (MQ=255)
                                                                 gTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTcgacgac                                                                         >  2:263203/1‑148 (MQ=255)
                                                                   cccTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcag                                                                                                                                   <  4:201120/88‑3 (MQ=255)
                                                                   cccTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCCATATcgcggcag                                                                                                                                   >  3:201120/1‑86 (MQ=255)
                                                                    ccTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGaa                                                                      >  6:337578/1‑148 (MQ=255)
                                                                     cTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAAc                                                                    >  5:99960/1‑149 (MQ=255)
                                                                        cAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAAcc                                                                   <  8:55120/147‑1 (MQ=255)
                                                                          attattAGTGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCGTGCTGACCGGCGAGACCTTGTCTGTCGACGACGAAACCTTcc                                                               <  7:96542/149‑1 (MQ=255)
                                                                                 gggACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTga                                                        <  5:337578/149‑1 (MQ=255)
                                                                                        tcCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCACTGTGAGAGcc                                                   >  1:405205/1‑147 (MQ=255)
                                                                                          ccGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGACACCTTCCGCAGTGAGAGCCTAgg                                               >  8:239846/1‑149 (MQ=255)
                                                                                          ccGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAgg                                               >  8:83283/1‑149 (MQ=255)
                                                                                          ccGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAgg                                               >  4:417929/1‑149 (MQ=255)
                                                                                           cGAATATTTGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCTATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCGGCTGTGAGAGCCTAGGc                                              <  8:213387/149‑1 (MQ=255)
                                                                                            gAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTc                                                                                                                            <  8:124925/70‑1 (MQ=255)
                                                                                            gAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGc                                              >  3:93249/1‑148 (MQ=255)
                                                                                            gAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGACGCGATATGGCGGCCAAGGATTc                                                                                                                            >  7:124925/1‑70 (MQ=255)
                                                                                              atatTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGa                                            <  4:50370/148‑1 (MQ=255)
                                                                                              atatTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGa                                            <  4:289529/148‑1 (MQ=255)
                                                                                               tatTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGt                                          <  5:115514/149‑1 (MQ=255)
                                                                                                 ttGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTg                                        <  1:263203/149‑1 (MQ=255)
                                                                                                     ccGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTc                                     >  2:215492/1‑148 (MQ=255)
                                                                                                       ggACTTGTTGTAACCCCTCTCCCCTATATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGAt                                  <  5:101252/149‑1 (MQ=255)
                                                                                                        gACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATg                                 <  7:71764/149‑1 (MQ=255)
                                                                                                          cTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGAc                               <  7:81728/149‑1 (MQ=255)
                                                                                                             gttgtAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGAcc                              <  1:68944/147‑1 (MQ=255)
                                                                                                                    cccctctcccctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAAc                     >  4:67179/1‑149 (MQ=255)
                                                                                                                      cctctcccctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCGTGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGt                                          <  2:212257/126‑1 (MQ=255)
                                                                                                                      cctctcccctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCGTGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGt                                          >  1:212257/1‑126 (MQ=255)
                                                                                                                            ccctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTg                <  5:205066/146‑1 (MQ=255)
                                                                                                                            ccctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGAACTCGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTc             >  6:403170/1‑149 (MQ=255)
                                                                                                                             cagcttttCACTGCCGCGATATGGCGTCCAAGGATTCCCCATGTCTGAATCACGCCTGCGTACCGGCGAGACCTGGTCTGGCTACGACGAAACCTTCCGCGGTGAGAGCCTAGGCGATTTGCTCGATGACCACGAAGAACTCGTGGTCg            <  8:363651/143‑1 (MQ=255)
                                                                                                                              ctctATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTCgg           >  2:223243/1‑149 (MQ=255)
                                                                                                                                  aTTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGCAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGACCTCGTGGTCGGCTc        >  3:158191/1‑148 (MQ=255)
                                                                                                                                   ttCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGATCGATGACCACTAAGAACTCGAGGTAGGCTCCa      >  8:114494/1‑149 (MQ=255)
                                                                                                                                      aCTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTCGGCTCCACCg   >  6:400956/1‑149 (MQ=255)
                                                                                                                                       cTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTCGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTCGGCTCCACCGt  <  7:400691/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
ACAGCCCTCCTGCAACCGTCTGTTGGATGCCTGTAAAAAAATAGCCAAAGGCCATCAGCATTCAAGTCCCTCCAATTATTAGGGACTCTCCCGAATATTGGCCGGACTTGTTGTAACCCCTCTCCCCTCTATTCACTGCCGCGATATGGCGGCCAAGGATTCCCCATGTCTGAATCACGCCTGCTGACCGGCGAGACCTGGTCTGTCGACGACGAAACCTTCCGCTGTGAGAGCCTAGGCGAGTTGCTCGATGACCACGAAGAACTCGTGGTCGGCTCCACCGT  >  NC_002947/4422948‑4423231

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: