Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,834 A→G 7.1% 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.1% 202.7 / 4.9 70I600T (ATT→ACT) PP_1963hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (2/3);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTC  >  NC_002947/2224686‑2224982
                                                                                                                                                     |                                                                                                                                                    
tAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCaa                                                                                                                                                      <  3:94789/149‑1 (MQ=255)
     gggCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGa                                                                                                                                                 <  5:325072/149‑1 (MQ=255)
       gCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGc                                                                                                                                               <  3:261308/149‑1 (MQ=255)
        cGTTCCGCCGTCTATATTGCGAAGCCTTCTTGTTTGTGTCTTCATCATGATTGTCGTTTATGAGAAACTCCAAAATTTTAGCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCt                                                                                                                                              <  7:50290/149‑1 (MQ=255)
         gTTCCGCCGTCTATATCGCGAAGCATGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGCtg                                                                                                                                             <  7:268116/149‑2 (MQ=255)
            ccgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGAc                                                                                                                                          >  2:15505/1‑149 (MQ=255)
             cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAg                                                                                                                                                >  3:22400/1‑142 (MQ=255)
             cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAg                                                                                                                                                <  4:22400/142‑1 (MQ=255)
             cgccgTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACg                                                                                                                                         >  1:146981/1‑149 (MQ=255)
                    tataTCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCAt                                                                                                                                   >  8:299168/1‑148 (MQ=255)
                      tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATATTGAGTATGAGc                                                                                                                                               <  8:296353/134‑1 (MQ=255)
                      tatCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATc                                                                                                                                  <  8:287758/147‑1 (MQ=255)
                        tCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                               >  2:362092/1‑148 (MQ=255)
                           cgAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                           <  1:241836/149‑1 (MQ=255)
                           cgAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAgg                                                                                                                           >  4:299974/1‑149 (MQ=255)
                              aGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTa                                                                                                                        <  1:362092/149‑1 (MQ=255)
                                       gATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCtt                                                                                                               >  2:151659/1‑149 (MQ=255)
                                         tGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                              <  2:125892/148‑1 (MQ=255)
                                           gtttCTTCATCATGATTGTCGTTGAAGATAAACACCAAAATTTTACCTTCCCACCGCAGGTCAGAGGCCCGTCCAA‑CTACGCTTTCATTTTCCACCTCATACTCAAAATGAGCTAAACGCTCCATCAAATTGTGTAGTGAGTCttaga                                                                                                            >  8:96095/4‑145 (MQ=255)
                                           gtgtCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGt                                                                                                            <  1:15505/148‑1 (MQ=255)
                                                 tcatcaTGATTGTCGTTGGTGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                      >  8:260274/1‑148 (MQ=255)
                                                  catcatGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTc                                                                                                      <  5:168120/147‑1 (MQ=255)
                                                        gATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                               >  3:267403/1‑116 (MQ=255)
                                                        gATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCaaa                                                                                                                               <  4:267403/116‑1 (MQ=255)
                                                         aTTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTggg                                                                                             >  8:242717/1‑149 (MQ=255)
                                                             tCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGagca                                                                                         >  5:254631/1‑149 (MQ=255)
                                                                    tgaGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGAt                                                                                   <  7:299168/148‑1 (MQ=255)
                                                                         aaCTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACtt                                                                             <  7:242717/149‑1 (MQ=255)
                                                                         aaCTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACtt                                                                             <  2:107428/149‑1 (MQ=255)
                                                                            tCCAAAATTTTACCTTCCGCGCGCAGGTCTGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCTTCAAATAGGGTAGTGAGTCTTCGTAAATTCGaa                                                                                                   <  8:248098/124‑1 (MQ=255)
                                                                            tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                             <  4:228413/98‑1 (MQ=255)
                                                                            tCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATa                                                                                                                             >  3:228413/1‑98 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGCCGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                              >  5:68688/1‑127 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                              <  5:329428/111‑1 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTc                                                                                                              >  6:329428/1‑111 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                              <  6:68688/127‑1 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                        >  2:12248/1‑149 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                        >  1:290570/1‑149 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                        >  7:129495/1‑149 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                        >  4:103752/1‑149 (MQ=255)
                                                                              cAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTa                                                                        >  4:339966/1‑149 (MQ=255)
                                                                               aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGa                                                                                                                    <  5:91150/104‑1 (MQ=255)
                                                                               aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGa                                                                                                                    >  6:91150/1‑104 (MQ=255)
                                                                               aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                              >  1:203524/1‑126 (MQ=255)
                                                                               aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTgg                                                                                              <  2:203524/126‑1 (MQ=255)
                                                                               aaaaTTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGagcagc                                                                                       >  7:129397/1‑133 (MQ=255)
                                                                               aaaaTTTTACCTTACCAGCGTAGGTCAGAGGACAGTGCAA‑CTTCGTTTTCATTTTGCAGCTCTTACTCAATTTCTGCTAGACGCTCCATCATATTGGGTAGTGAGTCTTCGTAAATTCGAAGCTGGGagcagc                                                                                       <  8:129397/133‑1 (MQ=255)
                                                                                          ttCCCAGCGCAGGTCAGAGGACCGTGCACACTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAt                                                             >  1:211263/1‑149 (MQ=255)
                                                                                            cccAGCGCAGGTCGTAGGACCGTCCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGATCCATCAAATAGGGTAGTGAGTCTTCGTAAATTAGAACCAGGGAGCAGCCGATGCACTTGTCTAAGCTGTAGCAtttc                                                          >  6:7380/1‑149 (MQ=255)
                                                                                                gcgcAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCAtttctttc                                                      >  8:283291/1‑149 (MQ=255)
                                                                                                  gcAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCa                                                    <  7:260274/149‑1 (MQ=255)
                                                                                                     ggTCATAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCTTCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                                 <  7:96095/149‑1 (MQ=255)
                                                                                                     ggTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                                 <  3:339966/149‑1 (MQ=255)
                                                                                                     ggTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                                 <  1:12248/149‑1 (MQ=255)
                                                                                                       tCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAgg                                               <  3:299974/149‑1 (MQ=255)
                                                                                                         agagGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                                 >  3:188538/1‑145 (MQ=255)
                                                                                                         agagGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTa                                                 <  4:188538/145‑1 (MQ=255)
                                                                                                           agGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                           >  8:84101/1‑149 (MQ=255)
                                                                                                           agGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAAcc                                           <  1:151659/149‑1 (MQ=255)
                                                                                                             gACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGa                                         <  7:283291/149‑1 (MQ=255)
                                                                                                               ccGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGAtt                                       <  2:290570/149‑1 (MQ=255)
                                                                                                                      c‑cTTCGCTTTCATTTTCCATCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGaa                                <  7:224661/149‑1 (MQ=255)
                                                                                                                      c‑cTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGaa                                >  3:24166/1‑149 (MQ=255)
                                                                                                                             cTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAATTTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCtt                          <  3:103752/149‑1 (MQ=255)
                                                                                                                              tttCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCCAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTc                         >  7:334002/1‑149 (MQ=255)
                                                                                                                                  aTTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGACCCTGGGAGCAGCCGATGCACTTGTCTATGGTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGt                     >  8:105298/1‑149 (MQ=255)
                                                                                                                                       ctaGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGt                >  5:357829/3‑149 (MQ=255)
                                                                                                                                        taGCTCATATTGAGTATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCTTTTCTTTCCATTAGGACCCGATTGCTCGAATCACTTCCAGTGTGGtg               >  8:85525/2‑148 (MQ=255)
                                                                                                                                        cAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt               <  6:254631/149‑1 (MQ=255)
                                                                                                                                         aGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGtt               >  6:32029/1‑148 (MQ=255)
                                                                                                                                                  tCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGAtt     >  3:25173/1‑149 (MQ=255)
                                                                                                                                                     atatGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTc  <  7:84101/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGCGGGGCGTTCCGCCGTCTATATCGCGAAGCCTGCTTGATGGTGTCTTCATCATGATTGTCGTTGATGAGAAACTCCAAAATTTTACCTTCCCAGCGCAGGTCAGAGGACCGTGCAC‑CTTCGCTTTCATTTTCCAGCTCATACTCAATATGAGCTAGACGCTCCATCAAATAGGGTAGTGAGTCTTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTC  >  NC_002947/2224686‑2224982

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