Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,520,081 T→A 5.4% intergenic (+700/‑180) PP_5580 → / → PP_5582 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,520,0810TA 5.4% 234.1 / 2.3 74intergenic (+700/‑180)PP_5580/PP_5582hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base A (2/2);  total (38/36)
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.62e-01

GCTTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGT  >  NC_002947/3519933‑3520227
                                                                                                                                                    |                                                                                                                                                  
gCTTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAAtgcttgct                                                                                                                                                    >  2:149211/1‑149 (MQ=255)
 cTTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAAtgcttgct                                                                                                                                                    >  3:223822/1‑148 (MQ=255)
     gACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATg                                                                                                                                                <  5:73112/148‑1 (MQ=255)
                    gCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTcc                                                                                                                                >  5:272078/1‑149 (MQ=255)
                    gCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTcc                                                                                                                                >  4:121313/1‑149 (MQ=255)
                        ttAGGGCGTCTTGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTTAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTg                                                                                                                              <  8:344529/147‑1 (MQ=255)
                          aGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTa                                                                                                                           >  6:69860/1‑148 (MQ=255)
                           gggCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTa                                                                                                                           >  4:289688/1‑147 (MQ=255)
                                  tGGGGTCATATTTGAGAATCGAGTCCCAGTTATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAatat                                                                                                                  >  3:303001/1‑149 (MQ=255)
                                  tGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAatat                                                                                                                  >  4:352583/1‑149 (MQ=255)
                                       tCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAtt                                                                                                             <  3:121313/149‑1 (MQ=255)
                                          tatTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAg                                                                                                          >  4:58502/1‑149 (MQ=255)
                                            tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGTTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGc                                                                                                         >  1:146189/1‑148 (MQ=255)
                                            tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCATAagcaagcat                                                                                                                                         >  7:338189/1‑115 (MQ=255)
                                            tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCATAagcaagcat                                                                                                                                         <  8:338189/116‑2 (MQ=255)
                                            tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCATAagcaagcat                                                                                                                                         <  6:92200/116‑2 (MQ=255)
                                            tttGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCATAagcaagcat                                                                                                                                         >  5:92200/1‑115 (MQ=255)
                                              tGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCga                                                                                                       <  1:149211/148‑1 (MQ=255)
                                               gagaATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGTTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCgag                                                                                                      <  2:146189/148‑1 (MQ=255)
                                                 gaATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATg                                                                                                   >  2:12560/1‑149 (MQ=255)
                                                  aaTCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATgg                                                                                                  >  2:21954/1‑149 (MQ=255)
                                                    tCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGc                                                                                                 >  7:86658/1‑148 (MQ=255)
                                                       aGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGaaa                                                                                             <  3:289688/149‑1 (MQ=255)
                                                           ccAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGt                                                                                         >  6:92846/1‑149 (MQ=255)
                                                           ccAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGt                                                                                         >  7:285321/1‑149 (MQ=255)
                                                            cAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTg                                                                                        >  6:195903/1‑149 (MQ=255)
                                                              gATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGt                                                                                      >  1:284249/1‑149 (MQ=255)
                                                                  tatCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTgg                                                                                   <  3:58502/148‑1 (MQ=255)
                                                                  tatCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGc                                                                                  <  6:272078/149‑1 (MQ=255)
                                                                    tCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATATGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCtt                                                                                <  5:189257/149‑1 (MQ=255)
                                                                    tCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCa                                                                                                               >  2:166491/1‑118 (MQ=255)
                                                                    tCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCa                                                                                                               <  1:166491/118‑1 (MQ=255)
                                                                      aTGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTg                                                                               >  4:134398/1‑148 (MQ=255)
                                                                       tGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTGTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTa                                                                             <  4:136164/149‑1 (MQ=255)
                                                                       tGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTa                                                                             <  8:285321/149‑1 (MQ=255)
                                                                        gccgccTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGtt                                                                                     <  6:269063/140‑1 (MQ=255)
                                                                        gccgccTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGtt                                                                                     >  5:269063/1‑140 (MQ=255)
                                                                            ccTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTAttttt                                                                        <  4:303001/149‑1 (MQ=255)
                                                                              tGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAtt                                                                                                             <  3:310120/110‑1 (MQ=255)
                                                                              tGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAtt                                                                                                             >  4:310120/1‑110 (MQ=255)
                                                                               ggaggaCGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCAGTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                      >  1:119948/1‑148 (MQ=255)
                                                                                gaggaCGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                      <  3:352583/147‑1 (MQ=255)
                                                                                gaggaCGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTAc                                                                      <  4:223822/147‑1 (MQ=255)
                                                                                      gATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACaaaaa                                                              >  2:54762/1‑149 (MQ=255)
                                                                                      gATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGAAATAGTGCTCATGGTACAGAAATATTCATTCAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACaaaaa                                                              >  6:53840/1‑149 (MQ=255)
                                                                                        tGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTTTTTTTACGCCAAAAAtt                                                            >  4:82660/1‑149 (MQ=255)
                                                                                          gTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTCTTACGACAAAAATTc                                                           >  2:199423/1‑148 (MQ=255)
                                                                                              atatTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                       <  5:92846/148‑1 (MQ=255)
                                                                                              atatTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAg                                                      >  3:129524/1‑149 (MQ=255)
                                                                                               tatTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGt                                                     <  1:7742/149‑1 (MQ=255)
                                                                                                 ttGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTg                                                    <  5:69860/148‑1 (MQ=255)
                                                                                                 ttGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTg                                                    >  1:158462/1‑148 (MQ=255)
                                                                                                    gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCAcc                                                <  1:21954/149‑1 (MQ=255)
                                                                                                    gTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATCCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCAc                                                 <  8:86658/148‑1 (MQ=255)
                                                                                                        tATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCttt                                             >  2:182183/1‑148 (MQ=255)
                                                                                                           gAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACa                                          >  6:125648/1‑148 (MQ=255)
                                                                                                             atatAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAg                                       >  8:104689/1‑149 (MQ=255)
                                                                                                                taGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                       >  4:82195/1‑130 (MQ=255)
                                                                                                                taGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCaaaa                                                       <  3:82195/130‑1 (MQ=255)
                                                                                                                taGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTc                                     <  5:189494/148‑1 (MQ=255)
                                                                                                                  ggTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTa                                   <  5:195903/148‑1 (MQ=255)
                                                                                                                      tgAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGt                              <  2:284249/149‑1 (MQ=255)
                                                                                                                         agcagcACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGa                           <  1:12560/149‑1 (MQ=255)
                                                                                                                          gcagcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGaa                          >  4:112769/1‑149 (MQ=255)
                                                                                                                            agcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAAt                         <  1:54762/148‑1 (MQ=255)
                                                                                                                             gcaCGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTATAGCTGTCGAAtt                        >  4:346937/1‑148 (MQ=255)
                                                                                                                                 ggATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGc                   <  4:129524/149‑1 (MQ=255)
                                                                                                                                   aTAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGt                 <  5:125648/149‑1 (MQ=255)
                                                                                                                                      gcactattgcttgcaTAAGCAAGCATTAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAt              >  6:271841/8‑149 (MQ=255)
                                                                                                                                        aTCAATGCTTGCTTATGCAAGCAATTGTTCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGTTTGGCTTGTATTTTTACGACAATAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAtt             <  3:60678/148‑1 (MQ=255)
                                                                                                                                        aTCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttg            >  5:149649/1‑149 (MQ=255)
                                                                                                                                         tCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGt                                                                              <  3:352266/82‑1 (MQ=255)
                                                                                                                                         tCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGt                                                                              >  4:352266/1‑82 (MQ=255)
                                                                                                                                          tattgcttgcaTAAGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgtt          >  7:286481/4‑149 (MQ=255)
                                                                                                                                           aaTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgttt         >  6:325661/1‑149 (MQ=255)
                                                                                                                                           aaTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGAttgttt         <  6:149649/149‑1 (MQ=255)
                                                                                                                                            aTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCt                                    >  2:189569/1‑121 (MQ=255)
                                                                                                                                            aTGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCt                                    <  1:189569/121‑1 (MQ=255)
                                                                                                                                               cttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATtgtttgtt      <  3:112769/148‑1 (MQ=255)
                                                                                                                                               cttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATtgtttgtt      <  5:325661/148‑1 (MQ=255)
                                                                                                                                               attgcttATGCAAGCTATAGTGCTCCTTGTACAGTAATATTCATTGAGCTATATGTCGAAAGCGTGGTTGGCTTGTATTTTTACTACAAAAATTCTATAGTGCTCCTTTACAAAGTGTAGCTGGCGTATTTATGCGTGATtgtttgtt      <  5:53840/147‑1 (MQ=255)
                                                                                                                                                ttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTc                                     <  2:155414/116‑1 (MQ=255)
                                                                                                                                                ttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTc                                     >  1:155414/1‑116 (MQ=255)
                                                                                                                                                ttgcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTtgttgt    >  1:228691/1‑149 (MQ=255)
                                                                                                                                                  gcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTgtgt  >  4:218005/1‑149 (MQ=255)
                                                                                                                                                  gcttATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTgtgt  >  1:37274/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTTTGACGATGAAGGAAAGGCTATTAGGGCGTCTGGGGTCATATTTGAGAATCGAGTCCCAGATATATCATGCCGCCTGGAGGACGATGGTGAATATTGGTGCTATGAATATAGGTGTGAAGCAGCACGGATAGTATCAATGCTTGCTTATGCAAGCAATAGTGCTCCTGGTACAGAAATATTCATTGAGCGAGATGGCGAAAGCGTGGTTGGCTTGTATTTTTACGACAAAAATTCAAAAGTGCACCTTTACAAAGTCTAGCTGTCGAATTTATGCGTGATTGTTTGTTGTGT  >  NC_002947/3519933‑3520227

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: