Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,832 C→G 7.6% E601Q (GAG→CAG) ‡ PP_1963 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,8320CG 7.6% 191.5 / 5.7 65E601Q (GAG→CAG) ‡PP_1963hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base G (2/3);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATT  >  NC_002947/2224689‑2224979
                                                                                                                                               |                                                                                                                                                   
cGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAatat                                                                                                                                                >  7:329645/1‑149 (MQ=255)
     cGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                           <  7:248738/149‑1 (MQ=255)
       ttCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGc                                                                                                                                            >  5:351147/1‑146 (MQ=255)
         ccgcccTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAGCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGAc                                                                                                                                       >  4:252211/1‑149 (MQ=255)
          cgccgTCTATATCGCGAAGCCTGCTTGATGGAGTCTTCATCATGATTGTCGTTGATGAGAAACTACAAAATTTTACATTCCCAGCGCAGGACAGAGGACCGTGCACCTTCGCTTTCATTTTACAGCTAATACTCAATATGAGATAGaag                                                                                                                                      >  5:377106/1‑147 (MQ=255)
                   tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCAt                                                                                                                                <  1:321625/146‑1 (MQ=255)
                     tCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                            <  4:72658/148‑1 (MQ=255)
                     tCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATAAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAAt                                                                                                                           <  4:148225/149‑1 (MQ=255)
                       gcgAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAg                                                                                                                         <  5:254084/149‑1 (MQ=255)
                         gAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAggg                                                                                                                       <  2:72353/149‑1 (MQ=255)
                                 cTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGa                                                                                                                 <  5:98963/147‑1 (MQ=255)
                                      tGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCg                                                                                                          >  6:251702/1‑149 (MQ=255)
                                      tGGCGTCTTCATCATGATTGTCGTTGATCAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           <  8:70511/148‑1 (MQ=255)
                                           tCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAAt                                                                                                     >  8:332232/1‑149 (MQ=255)
                                             ttCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                   >  3:177370/1‑149 (MQ=255)
                                             ttCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                   >  1:45987/1‑149 (MQ=255)
                                              tcatcaTGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                   >  7:183683/1‑148 (MQ=255)
                                               catcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAAtt                                                                                                    <  5:251702/146‑1 (MQ=255)
                                                atcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGaa                                                                                                >  2:142914/1‑149 (MQ=255)
                                                 tcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCt                                                                                                             <  3:355324/135‑1 (MQ=255)
                                                 tcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCt                                                                                                             >  4:355324/1‑135 (MQ=255)
                                                   atGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCt                                                                                             >  3:316215/1‑149 (MQ=255)
                                                    tGATTGTCGTTTATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCATATAGGGTAGTGAGTCTTCGTAAATTCGAACCTg                                                                                            <  3:252211/149‑1 (MQ=255)
                                                    tGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTg                                                                                            <  6:351147/149‑1 (MQ=255)
                                                        tGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCTAATAGGGTAGTGAGTCTTCGTAAAGTCGAACCTGGGag                                                                                        >  4:65294/1‑149 (MQ=255)
                                                               gatgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                           <  6:360024/91‑1 (MQ=255)
                                                               gatgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                           >  5:360024/1‑91 (MQ=255)
                                                                    gaAACTCCAAATTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGGACCTTCGATTTCATTTTCGAGCTCTTACTCAATATGAGCTAGACGCTCCATCAAATAGTGTATTGAGTCTTCGTAATTTCGAACGTGGGAGCAGCCGATGCAc                                                                            <  6:377106/149‑1 (MQ=255)
                                                                    gaAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCAc                                                                            >  7:50896/1‑149 (MQ=255)
                                                                          ccAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTc                                                                       >  5:280195/1‑148 (MQ=255)
                                                                           cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                     >  8:186001/1‑149 (MQ=255)
                                                                            aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGCGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTAt                                                                    >  5:337765/1‑149 (MQ=255)
                                                                              aaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGc                                                                  <  8:329645/149‑1 (MQ=255)
                                                                                ttttACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTg                                                                >  7:117335/1‑149 (MQ=255)
                                                                                       ttCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           >  8:272649/1‑99 (MQ=255)
                                                                                       ttCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                           <  7:272649/99‑1 (MQ=255)
                                                                                        tCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGcc                                                                                   >  5:177745/1‑122 (MQ=255)
                                                                                        tCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGcc                                                                                   <  6:177745/122‑1 (MQ=255)
                                                                                            aGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAtttcttt                                                    <  7:186001/149‑1 (MQ=255)
                                                                                                cAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCCTTCCAt                                                <  4:177370/149‑1 (MQ=255)
                                                                                                   gTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCATCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAg                                             >  3:224601/1‑149 (MQ=255)
                                                                                                   gTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCAt                                                >  6:159072/1‑146 (MQ=255)
                                                                                                   gTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  5:159072/148‑1 (MQ=255)
                                                                                                    tCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTAAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGATCCTGGGAGAAGCCGATGAACTTGTCTATGCTGTAGCATTTCTTTCCATTAg                                             >  3:143167/1‑148 (MQ=255)
                                                                                                        agGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                        >  4:270889/1‑149 (MQ=255)
                                                                                                            ccGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGAtt                                    <  1:142914/149‑1 (MQ=255)
                                                                                                                  aCCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGTACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              <  3:308166/133‑1 (MQ=255)
                                                                                                                  aCCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGTACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                              >  4:308166/1‑133 (MQ=255)
                                                                                                                           ttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTc                      <  7:206107/148‑1 (MQ=255)
                                                                                                                           ttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTc                      >  8:206107/1‑148 (MQ=255)
                                                                                                                             cATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTGGAATCCCTTCCAg                   >  8:60612/1‑149 (MQ=255)
                                                                                                                               ttttCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTg                                                                >  6:70423/1‑102 (MQ=255)
                                                                                                                               ttttCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTg                                                                <  5:70423/102‑1 (MQ=255)
                                                                                                                                   ctaGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                        >  4:178887/3‑122 (MQ=255)
                                                                                                                                   ctaGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                        <  3:178887/120‑1 (MQ=255)
                                                                                                                                   ccAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTgg              <  8:50896/148‑1 (MQ=255)
                                                                                                                                    cAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt            >  2:103803/1‑149 (MQ=255)
                                                                                                                                    cAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGt             <  2:45987/148‑1 (MQ=255)
                                                                                                                                        tCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTcc                                                  <  3:53003/107‑1 (MQ=255)
                                                                                                                                        tCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTcc                                                  >  4:53003/1‑107 (MQ=255)
                                                                                                                                          aTACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATaat      >  7:71841/1‑147 (MQ=255)
                                                                                                                                          aTACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATAc       <  4:316215/148‑1 (MQ=255)
                                                                                                                                          aTACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACt      >  1:224047/1‑149 (MQ=255)
                                                                                                                                            aCTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGCGAGCAGCCGATGCACTAGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCCATTGCTCGAATCCCTTCCAGTCTCGTAGATACTGa    >  5:50921/1‑149 (MQ=255)
                                                                                                                                              tCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAtt  >  7:63011/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATT  >  NC_002947/2224689‑2224979

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: ‑