Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,834 A→G 7.8% I600T (ATT→ACT)  PP_1963 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,8340AG 7.8% 180.1 / 3.9 64I600T (ATT→ACT) PP_1963hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/26);  new base G (2/3);  total (35/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 = 1.00e+00

CGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTC  >  NC_002947/2224689‑2224982
                                                                                                                                                 |                                                                                                                                                    
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                                                                                                      >  1:45987/1‑149 (MQ=255)
                                             ttCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                      >  3:177370/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                         >  8:206107/1‑148 (MQ=255)
                                                                                                                           ttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTc                         <  7:206107/148‑1 (MQ=255)
                                                                                                                             cATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTGGAATCCCTTCCAg                      >  8:60612/1‑149 (MQ=255)
                                                                                                                               ttttCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTg                                                                   <  5:70423/102‑1 (MQ=255)
                                                                                                                               ttttCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTg                                                                   >  6:70423/1‑102 (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                                                     >  4:53003/1‑107 (MQ=255)
                                                                                                                                        tCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTcc                                                     <  3:53003/107‑1 (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)
                                                                                                                                                aaTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAtttt   >  1:165066/1‑149 (MQ=255)
                                                                                                                                                 atatGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTc  <  1:103803/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCACCTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTC  >  NC_002947/2224689‑2224982

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