Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,497 T→A 6.6% intergenic (+21/‑131) PP_3010 → / → PP_3011 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,4970TA 6.6% 177.9 / 3.1 61intergenic (+21/‑131)PP_3010/PP_3011hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (2/2);  total (31/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 = 6.88e-01

ACAAGGCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATG  >  NC_002947/3402349‑3402640
                                                                                                                                                    |                                                                                                                                               
aCAAGGCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAAt                                                                                                                                                 <  8:146222/149‑1 (MQ=255)
   aGGCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGc                                                                                                                                              <  5:150475/149‑1 (MQ=255)
   aGGCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGc                                                                                                                                              <  4:261534/149‑1 (MQ=255)
     gCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGccc                                                                                                                                            <  3:31108/149‑1 (MQ=255)
     gCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGccc                                                                                                                                            >  6:89565/1‑149 (MQ=255)
     gCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGccc                                                                                                                                            <  5:109485/149‑1 (MQ=255)
         agGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCccgccg                                                                                                                                        >  4:190495/1‑149 (MQ=255)
           gggCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCccgctctc                                                                                                                                      >  4:125825/1‑145 (MQ=255)
              cTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGc                                                                                                                                              >  1:69278/1‑138 (MQ=255)
              cTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGc                                                                                                                                              <  2:69278/138‑1 (MQ=255)
                         gCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCtat                                                                                                                         <  3:125825/148‑3 (MQ=255)
                             ggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtg                                                                                                                    >  1:230415/1‑149 (MQ=255)
                              gCACGGGAAGCGAGCGCCACTTCAGGGCCTATCAACGTCATGCCTTGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTTATCATCCAACCCGCACACGAATAGCCCGCCGTCGTGAGCGGGCTTTTtgtgt                                                                                                                   <  3:187453/149‑1 (MQ=255)
                              gCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAAACATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgt                                                                                                                   >  1:5734/1‑149 (MQ=255)
                                    gAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCAt                                                                                                             >  3:74170/1‑149 (MQ=255)
                                        cgagcgCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTc                                                                                                                                      <  2:157785/120‑1 (MQ=255)
                                        cgagcgCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTc                                                                                                                                      >  1:157785/1‑120 (MQ=255)
                                         gagcgCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCt                                                                                                        <  2:230415/149‑1 (MQ=255)
                                              ccACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtg                                                                                                                      >  3:250114/1‑130 (MQ=255)
                                              ccACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtg                                                                                                                      <  4:250114/130‑1 (MQ=255)
                                               cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCgg                                                                                                  >  7:268879/1‑149 (MQ=255)
                                                    cACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAAc                                                                                              <  2:5734/148‑1 (MQ=255)
                                                        gCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTc                                                                                          >  1:54201/1‑148 (MQ=255)
                                                               aaCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAg                                                                                  >  7:168697/1‑149 (MQ=255)
                                                               aaCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAg                                                                                  >  2:177622/1‑149 (MQ=255)
                                                               aaCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAg                                                                                  >  7:283108/1‑149 (MQ=255)
                                                                    gATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGAgg                                                                             >  6:52007/1‑149 (MQ=255)
                                                                     aTGCCTGGTACGACCATGGCAGGTTGCTTGATTCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCggg                                                                            <  7:99617/149‑1 (MQ=255)
                                                                     aTGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCggg                                                                            <  1:212158/149‑1 (MQ=255)
                                                                      tGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGt                                                                           <  1:177622/149‑1 (MQ=255)
                                                                        ccTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTg                                                                          >  3:185428/1‑148 (MQ=255)
                                                                          tGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATTGCCTGCCCGGCAACGTTCAAAAAAAGGGc                                                                               <  6:26003/141‑1 (MQ=255)
                                                                          tGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATTGCCTGCCCGGCAACGTTCAAAAAAAGGGc                                                                               >  5:26003/1‑141 (MQ=255)
                                                                          tGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAg                                                                       >  7:98432/1‑149 (MQ=255)
                                                                             tACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                    <  5:89565/149‑1 (MQ=255)
                                                                                gACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTa                                                                 >  5:77428/1‑149 (MQ=255)
                                                                                   cATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATgg                                                              <  8:268879/149‑1 (MQ=255)
                                                                                        cAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGa                                                         <  3:190495/149‑1 (MQ=255)
                                                                                         aGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGag                                                        <  4:74170/149‑1 (MQ=255)
                                                                                                        gcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGAGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGc                                          <  4:185428/148‑1 (MQ=255)
                                                                                                         cACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGc                                          <  5:52007/147‑1 (MQ=255)
                                                                                                           cTGATGTGTTCAATAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAATAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGGTGCTTTGCAttt                                      <  2:206675/149‑1 (MQ=255)
                                                                                                              aTGTGTTCAAAAGCTTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGc                                    <  8:283108/148‑1 (MQ=255)
                                                                                                                 tgtTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGc                                >  8:258548/1‑149 (MQ=255)
                                                                                                                  gtTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcc                               >  1:233191/1‑149 (MQ=255)
                                                                                                                   ttaaaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCTGGTGGAGTCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                              >  1:138840/4‑149 (MQ=255)
                                                                                                                   ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                              <  6:77428/149‑1 (MQ=255)
                                                                                                                   ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCAGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGccc                              >  5:174100/1‑149 (MQ=255)
                                                                                                                        aaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgccc                         >  2:265027/1‑149 (MQ=255)
                                                                                                                          gCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgcccg                        >  1:111030/1‑148 (MQ=255)
                                                                                                                          gCGTAATAATACAACCCCCACAACAATAGCCAGCCGTCGAGAGCGGTCTTTTTGTGTGCCCATCGCCTGACCGCCAACGTTAAAAAAAAGGGCGGGTGCAGCCCGTATGGCCTGAGCAGCGATGCTTTGCATTTGCCGGCaccgccag                        >  2:111033/1‑146 (MQ=255)
                                                                                                                           cGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTa                                                                 <  2:256437/106‑1 (MQ=255)
                                                                                                                              aaTCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAata                   >  4:219136/1‑149 (MQ=255)
                                                                                                                                      aaCCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGt           <  8:168697/149‑1 (MQ=255)
                                                                                                                                       atgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCtt                                             <  3:261025/109‑1 (MQ=21)
                                                                                                                                       atgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCtt                                             >  4:261025/6‑114 (MQ=21)
                                                                                                                                       atgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCtt                                             <  5:181353/109‑1 (MQ=21)
                                                                                                                                       atgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTAATGGTGTGCCCATCGCCTGCCCGGCAACGTACAAAAAAAGGGCGGGTGGAGCCCGAATGGCCTGAGCAGCGTTGCtt                                             >  6:181353/6‑114 (MQ=18)
                                                                                                                                          cccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCATCGTTCAAAAAAAGGGCGTGTGGAGCCGGTTTGGCCTGAGCATCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt       <  1:111033/149‑1 (MQ=255)
                                                                                                                                          cccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt       <  1:265027/149‑1 (MQ=255)
                                                                                                                                           ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGGAACGTTGAAAAACAGGGCGGCTGGAGCCGGTATGGCCTGAGCAGAGTTGCTTTGCATTTGGCGGCCACGCCCGTAGTCTTCGTGGttggtg      >  8:121050/1‑144 (MQ=255)
                                                                                                                                                 cAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatg  >  2:218172/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACAAGGCAGAGGGGCTTGTAGAGTTGCGGGGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATG  >  NC_002947/3402349‑3402640

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: