Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 860,698 G→A 7.7% intergenic (‑27/+77) PP_0741 ← / ← PP_0742 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947860,6980GA 7.7% 188.5 / 4.5 65intergenic (‑27/+77)PP_0741/PP_0742hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (3/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

TGTTCCAGCTGTTTACGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTCGACTCGGCC  >  NC_002947/860550‑860843
                                                                                                                                                    |                                                                                                                                                 
tGTTCCAGCTGTTTACGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCCGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACg                                                                                                                                                   >  2:363800/1‑149 (MQ=255)
         tGTTTACGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGa                                                                                                                                           >  5:308809/1‑148 (MQ=255)
              aCGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACtt                                                                                                                                      <  8:314942/148‑1 (MQ=255)
              aCGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACtt                                                                                                                                      <  8:234718/148‑1 (MQ=255)
                       gTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGttt                                                                                                                            <  5:372449/149‑1 (MQ=255)
                          ggAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTAtt                                                                                                                         >  6:184851/1‑149 (MQ=255)
                          ggAGGCCTTGTGGCCGACGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTAtt                                                                                                                         >  2:168565/1‑149 (MQ=255)
                             ggCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCa                                                                                                                      <  2:288121/149‑1 (MQ=255)
                              gCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGTAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAAGTTGTACAGGTTTATTTCa                                                                                                                      <  4:110305/148‑1 (MQ=255)
                              gCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACAAGTTTCGAGTCACGTACAGGTTTAttttcg                                                                                                                     <  2:269975/149‑4 (MQ=37)
                                cTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGca                                                                                                                   >  6:214915/1‑149 (MQ=255)
                                 ttGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGca                                                                                                                   >  6:11378/1‑148 (MQ=255)
                                      gccgccGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAgcgc                                                                                                             >  2:210782/1‑149 (MQ=255)
                                            gTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGt                                                                                                       <  1:82237/149‑1 (MQ=255)
                                            gTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGt                                                                                                       >  6:333883/1‑149 (MQ=255)
                                               gTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTa                                                                                                     <  5:214915/148‑1 (MQ=255)
                                                tAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTAt                                                                                                    <  8:94089/148‑1 (MQ=255)
                                                tAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATg                                                                                                   <  8:293750/149‑1 (MQ=255)
                                                  gcgcACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGc                                                                                                  <  6:143323/148‑1 (MQ=255)
                                                    gcACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCaaa                                                                                               >  1:165372/1‑149 (MQ=255)
                                                     cACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAAc                                                                                              >  7:143068/1‑149 (MQ=255)
                                                       cGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATACACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTa                                                                                            <  2:231600/149‑1 (MQ=255)
                                                            tgCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGTCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTa                                                                                         <  7:239584/147‑1 (MQ=255)
                                                            tgCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCCGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTa                                                                                         <  1:363800/147‑1 (MQ=255)
                                                              cccGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTa                                                                                            >  3:8083/1‑142 (MQ=255)
                                                              cccGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTa                                                                                            <  4:8083/142‑1 (MQ=255)
                                                              cccGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACAAGTTTCGAGTCACGTACAGGTTTAttttc                                                                                                                      >  2:265573/1‑114 (MQ=25)
                                                              cccGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCCAAAATAAACCTGTACAAGTTTCGAGTCACGTACAGGTTTAttttc                                                                                                                      <  1:265573/116‑3 (MQ=21)
                                                               ccGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTAttttt                                                                                    <  5:333883/149‑1 (MQ=255)
                                                                 ggTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTg                                                                                   >  5:164470/1‑148 (MQ=255)
                                                                 ggTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTg                                                                                   >  1:61406/1‑148 (MQ=255)
                                                                 ggTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGt                                                                                  <  7:54648/149‑1 (MQ=255)
                                                                     aggcaggcGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGtata                                                                               <  5:184851/148‑1 (MQ=255)
                                                                      ggcaggcGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGGACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGt                                                                             <  1:168565/149‑1 (MQ=255)
                                                                       gcaggcGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGt                                                                             >  4:88130/1‑148 (MQ=255)
                                                                         aggcGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGtat                                                                                >  3:176984/1‑143 (MQ=255)
                                                                         aggcGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGtat                                                                                <  4:176984/143‑1 (MQ=255)
                                                                            cccTGATGGCGATCCGCGGCTGGCTGGTGGCGGAAGTGTCGGGCATGGCTGCCGCCGAAAATACCCCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCa                                                                       <  8:255977/147‑1 (MQ=255)
                                                                             gcTGATGTCGATACGCGGCTTGCTGTGGGCGGAATGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTGGATACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCAt                                                                      <  5:10000/149‑1 (MQ=255)
                                                                                     cGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAg                                                              <  6:308809/149‑1 (MQ=255)
                                                                                     cGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAg                                                              <  1:344556/149‑1 (MQ=255)
                                                                                           gcgGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCa                                                                    >  7:260098/1‑137 (MQ=255)
                                                                                           gcgGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCa                                                                    <  8:260098/137‑1 (MQ=255)
                                                                                               cTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGc                                                    <  2:165372/149‑1 (MQ=255)
                                                                                                ttGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGc                                                    >  8:233131/1‑148 (MQ=255)
                                                                                                         gCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCcgcg                                          >  4:122112/1‑149 (MQ=255)
                                                                                                         gCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCcgcg                                          >  3:329793/1‑149 (MQ=255)
                                                                                                         gCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCcgcg                                          >  5:203618/1‑149 (MQ=255)
                                                                                                                ggTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTc                                   >  5:208972/1‑149 (MQ=255)
                                                                                                                ggTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTc                                   <  3:88130/149‑1 (MQ=255)
                                                                                                                    gTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTg                                <  1:198473/148‑1 (MQ=255)
                                                                                                                       cATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAAGTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCa                            >  2:72524/1‑149 (MQ=255)
                                                                                                                       cATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCa                            <  6:203618/149‑1 (MQ=255)
                                                                                                                       cATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCa                            <  4:329793/149‑1 (MQ=255)
                                                                                                                          ggCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCAc                                                                   >  4:325143/1‑107 (MQ=255)
                                                                                                                          ggCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCAc                                                                   <  3:325143/107‑1 (MQ=255)
                                                                                                                           gCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGa                        <  5:11378/149‑1 (MQ=255)
                                                                                                                                cTCCGAAAATAAACCTGTACGTTACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCAt                   <  3:122112/149‑1 (MQ=255)
                                                                                                                                cTCCGAAAATAAACCTCTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTACTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCAt                   >  5:114574/1‑149 (MQ=255)
                                                                                                                                 tcagAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCATGGCCACTGAGCCATg                  >  6:286547/4‑149 (MQ=255)
                                                                                                                                 tCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATg                  >  3:247501/1‑149 (MQ=255)
                                                                                                                                    tgaaaTAAACCTGTACAAGTTTCGAGTCACGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATgcgc               >  3:146430/3‑149 (MQ=37)
                                                                                                                                        aTAAACCTGTACAAGTTTCGAGTCACGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTc           >  5:333330/1‑149 (MQ=38)
                                                                                                                                             ccTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTCGACTc      >  2:103287/1‑149 (MQ=255)
                                                                                                                                              cTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTCGACTCg     <  7:233131/149‑1 (MQ=255)
                                                                                                                                                 tACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTCGACTCGGcc  >  4:316263/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGTTCCAGCTGTTTACGGCACAGGTCGGAGGCCTTGTGGCCGCCGTTGTAGCGCACGCTGTGCCCGGTCAGGCAGGCGCTGATGGCGATACGCGGCTTGCTGTGGGCGGAAGGGTCGTGCATGGCTGCCTCCGAAAATAAACCTGTACGTGACTCGAAACTTGTACAGGTTTATTTCAGCACAGCGCTGAAGTTATGCAAACTAGTATTTTTGTATAGTTCGGCATCACTCAAGGTGTTCGAGCAAGCGCCGCGCTGCTTCCTGGCCACTGAGCCATGCGCCCTCGACTCGGCC  >  NC_002947/860550‑860843

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: