Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 58,439 G→C 8.7% intergenic (‑350/+338) PP_5424 ← / ← PP_0050 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294758,4390GC 8.7% 184.5 / 5.2 63intergenic (‑350/+338)PP_5424/PP_0050hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base C (2/3);  total (27/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.93e-01

CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGC  >  NC_002947/58291‑58581
                                                                                                                                                     |                                                                                                                                                
cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                  >  2:17051/1‑149 (MQ=255)
         gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGAGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                          <  1:98913/149‑1 (MQ=255)
         gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                          >  4:291623/1‑148 (MQ=255)
         gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgc                                                                                                                                         <  5:143400/149‑2 (MQ=255)
          ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                        <  8:12397/149‑1 (MQ=255)
          ccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                        <  4:38755/149‑1 (MQ=255)
               cATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAagtagt                                                                                                                                   <  5:275661/149‑1 (MQ=255)
                 tCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGa                                                                                                                                 >  4:365281/1‑149 (MQ=255)
                    tCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcgg                                                                                                                              >  6:19211/1‑149 (MQ=255)
                       tGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcggcgg                                                                                                                          <  5:19211/149‑1 (MQ=255)
                         cATCAGATCATGATTGTCACGGGTTCCGCTTTGGATCTCCTTGCAATTGAGTGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcggcgg                                                                                                                          <  7:178645/147‑1 (MQ=255)
                          aTCAGATCTTGATTGACACGCCTGCCGCTTTGGATAGCCTTGAAATTGAGGGCTCCAACTGCGCATTTGTCGGGACTGCGTACCAGCATGATCCCCCTCAAAATTGGCAGCAGAGGAAGCTGGTTACTCTGAAGTAGTGAcggctcgca                                                                                                                        >  5:39400/1‑144 (MQ=255)
                             aGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTg                                                                                                                    >  6:289021/1‑149 (MQ=255)
                              gATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGc                                                                                                                   <  3:291623/149‑1 (MQ=255)
                                    gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTc                                                                                                              <  6:320228/147‑1 (MQ=255)
                                        gaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAt                                                                                                          <  3:306775/148‑1 (MQ=255)
                                         acacGCGTGCCGCTTTGGATAGCCTTGCAATTGCGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTCCCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                         <  4:231796/148‑1 (MQ=255)
                                           acGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCt                                                                                                       >  8:198448/1‑148 (MQ=255)
                                             gcgcgcCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                    >  4:294049/1‑149 (MQ=255)
                                             gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                    <  4:79221/149‑1 (MQ=255)
                                                tGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGTCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAccc                                                                                                <  8:177472/149‑1 (MQ=255)
                                                  ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAAccc                                                                                                >  7:52826/1‑147 (MQ=255)
                                                  ccGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCaa                                                                                              <  7:198448/149‑1 (MQ=255)
                                                   cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCTAGTTTGGCAGCGGAGGAAGCTGGTTCGTCGGAAGTAGTGTCGGCGGGCATTGCCAGTCGGAGTCTGAAcgctag                                                                                             >  8:297619/1‑144 (MQ=255)
                                                   cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                             >  1:281344/1‑149 (MQ=255)
                                                   cGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAg                                                                                             <  3:246534/149‑1 (MQ=255)
                                                    gCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAgg                                                                                            <  2:281344/149‑1 (MQ=255)
                                                    gCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGACGAAGCTGGTACCTCTGAAGTAGTGACGGCGGGCATTGACAGGCGGATACTGAACCCAAgg                                                                                            >  8:318443/1‑149 (MQ=255)
                                                         ggATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGc                                                                                        <  8:52826/148‑1 (MQ=255)
                                                            tAGACTTGCAATTGAGGGCTGAAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTACTCTGAAGTAGTGACGGCGGGCATTGCCACTCGGATTCTGCACCCAAGGTCgagag                                                                                     >  4:331665/1‑144 (MQ=255)
                                                                cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                <  2:289240/149‑1 (MQ=255)
                                                                 ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                >  2:132125/1‑148 (MQ=255)
                                                                 ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAg                                                                                >  8:123339/1‑148 (MQ=255)
                                                                   gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAg                                                                             >  1:243523/1‑149 (MQ=255)
                                                                      aTTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTTGCAGCAGAGGAACCAGCTTCCTCTg                                                                                                                                          <  8:125414/87‑1 (MQ=255)
                                                                      aTTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTTGCAGCAGAGGAACCAGCTTCCTCTg                                                                                                                                          >  7:125414/1‑87 (MQ=255)
                                                                            ggCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCCAAGGTCGCCAGACTAGGAGCTAAAGGa                                                                     <  1:17051/149‑1 (MQ=255)
                                                                            ggCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGAt                                                                    <  2:243523/149‑1 (MQ=255)
                                                                              cTGAAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCACCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCGAGTCGGAATCTGAACCCAAGGACGCCAGACTAGGAGCTAAAGGATCg                                                                  >  2:42583/1‑149 (MQ=255)
                                                                                    cTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGa                                                            <  2:205914/149‑1 (MQ=255)
                                                                                     tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGac                                                           >  8:276623/1‑149 (MQ=255)
                                                                                        gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGactac                                                        <  3:365281/149‑1 (MQ=255)
                                                                                          aTTTGTCGTGCCTGGGTACCAGAATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACTGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActc                                                      <  2:115239/149‑1 (MQ=255)
                                                                                          aTTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCCGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActc                                                      >  7:238665/1‑149 (MQ=255)
                                                                                               tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                                 >  7:289273/1‑149 (MQ=255)
                                                                                                cGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                                 <  6:326241/148‑1 (MQ=255)
                                                                                                 gggCGTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATcc                                               <  7:318443/149‑1 (MQ=255)
                                                                                                 gggCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATcc                                               <  3:294049/149‑1 (MQ=255)
                                                                                                               aGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGt                                  <  5:289021/149‑1 (MQ=255)
                                                                                                                gCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCGAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGt                                  >  5:123516/1‑148 (MQ=255)
                                                                                                                gCATGAGCCCCCACAAATTTCGCACCAGAGGATGCTGATACCTCTGAAGTAGATACGGCGGGCATTGCCAGACGGATACTGAACCCAAGGTCGCCAGACTATGAGCTAAAGGATAGCTGCTACTACTCCTCATCCCTTACAGaaaatt                                  >  6:210389/1‑146 (MQ=255)
                                                                                                                  aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                               <  1:366216/149‑1 (MQ=255)
                                                                                                                     aGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAcgc                             <  4:331121/148‑1 (MQ=255)
                                                                                                                       cccccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGCCTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCCAAAAGTCACGCgtc                          >  8:131040/1‑149 (MQ=255)
                                                                                                                       cccccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtc                          <  8:289273/149‑1 (MQ=255)
                                                                                                                        cccaGCAATTTTGTCAGCAGAGGAAGCTGGTTCCTCTTAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTTCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtc                          <  3:331665/148‑1 (MQ=255)
                                                                                                                         cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAg                                                                       >  4:109846/1‑102 (MQ=255)
                                                                                                                         cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAg                                                                       <  3:109846/102‑1 (MQ=255)
                                                                                                                                ttttGGCATCAGAGGTATCTGGTTCCTCTGTATTAGTTACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTaaa                 <  7:267072/149‑1 (MQ=255)
                                                                                                                                ttttGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTaaa                 <  7:123339/149‑1 (MQ=255)
                                                                                                                                    ggCATCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACgg              <  5:21026/148‑1 (MQ=255)
                                                                                                                                       ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGcccg          >  3:204340/3‑148 (MQ=255)
                                                                                                                                            agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc       >  2:277519/1‑147 (MQ=255)
                                                                                                                                                 gCTGGTTCCTCTGAAGAAGTGACGGCGGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGAATAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCATTACCGAAAAGACACGCGTCGGTCGTAAACGGCCCTATCCAGGc  >  1:197056/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCCAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGC  >  NC_002947/58291‑58581

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