Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,210,756 A→T 6.3% intergenic (+33/+36) PP_3692 → / ← PP_5621 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,210,7560AT 6.3% 182.7 / 2.8 63intergenic (+33/+36)PP_3692/PP_5621hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

TTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTG  >  NC_002947/4210610‑4210898
                                                                                                                                                  |                                                                                                                                                
tttGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAAc                                                                                                                                                >  8:299957/1‑149 (MQ=255)
      cATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGt                                                                                                                                          <  7:299957/149‑1 (MQ=255)
          cccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                       <  4:113087/148‑1 (MQ=255)
           ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                     <  4:180473/149‑1 (MQ=255)
            cAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGc                                                                                                                                    >  6:211978/1‑149 (MQ=255)
              gTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTGACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCt                                                                                                                                   >  1:308769/1‑148 (MQ=255)
               tCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTa                                                                                                                                 <  2:193296/149‑1 (MQ=255)
                 ggaggaGCTTCCCAGCAGTGCCTTGCAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGc                                                                                                                               <  3:119643/149‑1 (MQ=255)
                         ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTgg                                                                                                                       >  5:142033/1‑149 (MQ=255)
                           cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTgagc                                                                                                                                   <  3:135935/135‑4 (MQ=255)
                           cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTgagc                                                                                                                                   >  4:135935/1‑132 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTACCGAATTCTCAGCACGCAAGGAGCAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                     <  5:303799/148‑1 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTga                                                                                                                                     <  4:282117/132‑2 (MQ=255)
                            ccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTga                                                                                                                                     >  3:282117/1‑131 (MQ=255)
                              agcagTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgcagca                                                                                                                  >  7:313638/1‑149 (MQ=255)
                                 agTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATAc                                                                                                               <  8:313638/149‑1 (MQ=255)
                                         ggAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCaa                                                                                                        >  7:186677/1‑148 (MQ=255)
                                                aaGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGAGTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAgc                                                                                              <  8:280075/149‑1 (MQ=255)
                                                  ggTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAgcgg                                                                                            <  2:315467/149‑1 (MQ=255)
                                                         gAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTcc                                                                                     >  4:189887/1‑149 (MQ=255)
                                                                aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTa                                                                                <  2:400991/147‑1 (MQ=255)
                                                                aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTa                                                                                <  3:318582/147‑1 (MQ=255)
                                                                  gTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCAt                                                                            <  3:189887/149‑1 (MQ=255)
                                                                         cactgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAAt                                                                      <  5:211978/148‑1 (MQ=255)
                                                                           ctgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATAcc                                                                   >  2:158343/1‑149 (MQ=255)
                                                                           ctgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATAcc                                                                   >  1:175250/1‑149 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCAc                                                                                                           >  4:74406/1‑108 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCAc                                                                                                           <  3:74406/108‑1 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGa                                                                 >  1:189576/1‑148 (MQ=255)
                                                                              caTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGCCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGa                                                                 >  1:358143/1‑148 (MQ=255)
                                                                               aTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                >  1:274048/1‑148 (MQ=255)
                                                                                tAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTAGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCa                                                              <  6:150247/149‑1 (MQ=255)
                                                                                   cTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGc                                                           >  1:123595/1‑149 (MQ=255)
                                                                                        ttAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCg                                                      >  6:233198/1‑149 (MQ=255)
                                                                                         tAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                     >  7:99735/1‑149 (MQ=255)
                                                                                          aaCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                     <  2:358143/148‑1 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                     >  7:56584/1‑69 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                     <  8:56584/69‑1 (MQ=255)
                                                                                           aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGc                                                    >  4:75655/1‑148 (MQ=255)
                                                                                            cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAC‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGa                                                  >  1:143159/1‑149 (MQ=255)
                                                                                             gAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGAt                                                 >  7:42177/1‑149 (MQ=255)
                                                                                               aTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGa                                               <  8:42177/149‑1 (MQ=255)
                                                                                                        aCGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGa                                      >  4:251981/1‑149 (MQ=255)
                                                                                                                aGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTc                              >  2:403435/1‑149 (MQ=255)
                                                                                                                   aGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGcc                           <  1:158343/149‑1 (MQ=255)
                                                                                                                   aGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGc                            <  3:75655/148‑1 (MQ=255)
                                                                                                                     cATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                         >  1:148993/1‑149 (MQ=255)
                                                                                                                     cATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                         <  2:175250/149‑1 (MQ=255)
                                                                                                                         cGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                      <  4:85992/148‑1 (MQ=255)
                                                                                                                         cGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                      <  2:123595/148‑1 (MQ=255)
                                                                                                                         cGGGGAGCGGGCTCAGAACCGCGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCAt                      <  4:80512/148‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                     <  4:269683/131‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                     >  3:269683/1‑131 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGaa                   <  5:233198/149‑1 (MQ=255)
                                                                                                                           gTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGaa                   <  8:186677/149‑1 (MQ=255)
                                                                                                                              aGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAAAGTGTAGCGGCGATTCCGATTACCATCCCAATACCGATAAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAg                  >  1:325107/1‑149 (MQ=255)
                                                                                                                                cGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGAc              >  2:258417/1‑149 (MQ=255)
                                                                                                                                  ggCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGAt                                                 <  4:324468/112‑1 (MQ=255)
                                                                                                                                   gCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCa           <  2:308769/149‑1 (MQ=255)
                                                                                                                                     tCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGc         >  1:355690/1‑149 (MQ=255)
                                                                                                                                     tCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGc         >  8:211449/1‑149 (MQ=255)
                                                                                                                                        gAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGa      >  1:148504/1‑149 (MQ=255)
                                                                                                                                               gTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTg  >  3:396/1‑146 (MQ=255)
                                                                                                                                                  aaCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGt                 <  3:59413/128‑1 (MQ=255)
                                                                                                                                                  aaCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGt                 >  4:59413/1‑128 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTTGAGCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAA‑‑CGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTG  >  NC_002947/4210610‑4210898

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: