Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,477,550 A→T 6.8% I139N (ATT→AAT)  PP_3088 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,477,5500AT 6.8% 173.2 / 2.6 59I139N (ATT→AAT) PP_3088hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base T (2/2);  total (28/31)
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.40e-01

CTGGAAGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAGTG  >  NC_002947/3477402‑3477692
                                                                                                                                                    |                                                                                                                                              
cTGGAAGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTTTTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCaa                                                                                                                                                <  1:241360/149‑1 (MQ=255)
   gAAGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGGTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGt                                                                                                                                             >  5:176664/1‑149 (MQ=255)
   gAAGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATg                                                                                                                                              <  4:292538/148‑1 (MQ=255)
     aGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCg                                                                                                                                           <  7:82579/149‑1 (MQ=255)
         ctcGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGt                                                                                                                                         >  7:157582/1‑147 (MQ=255)
            ggCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGTAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGACATTGCCGGTggc                                                                                                                                      <  5:296596/147‑3 (MQ=255)
              cGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                  <  7:166488/149‑1 (MQ=255)
              cGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                  <  5:296287/149‑1 (MQ=255)
               gTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGGGCGTgg                                                                                                                                 <  1:188973/149‑1 (MQ=255)
                ttGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGTCAAGCCACCGGCAATGTCGGTGAGCGTGGt                                                                                                                                <  6:139990/149‑1 (MQ=255)
                ttGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGt                                                                                                                                <  4:262215/149‑1 (MQ=255)
                    tCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                  <  1:48993/143‑1 (MQ=255)
                    tCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTg                                                                                                                                  >  2:48993/1‑143 (MQ=255)
                        gCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAg                                                                                                                          <  7:461/147‑1 (MQ=255)
                          cTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAg                                                                                                                      <  2:287688/149‑1 (MQ=255)
                                   gTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTcgc                                                                                                             >  8:31460/1‑149 (MQ=255)
                                      cTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTcgcgcc                                                                                                          >  8:128183/1‑149 (MQ=255)
                                           tGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCcgcgcg                                                                                                     >  6:135963/1‑149 (MQ=255)
                                           tGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCcgcgcg                                                                                                     >  3:14606/1‑149 (MQ=255)
                                             cATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCgcgcgc                                                                                                    <  7:128183/148‑1 (MQ=255)
                                             cATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCt                                                                                                   >  8:207435/1‑149 (MQ=255)
                                              aTCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTc                                                                                                  >  5:43063/1‑149 (MQ=255)
                                              aTCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTc                                                                                                  >  4:273638/1‑149 (MQ=255)
                                                    ttGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTcgcgccgcgc                                                                                                      >  1:283941/1‑139 (MQ=255)
                                                    ttGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTcgcgccgcgc                                                                                                      <  2:283941/139‑1 (MQ=255)
                                                      ggggTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAg                                                                                            <  7:31460/147‑1 (MQ=255)
                                                               cAAGGCGAGGTTCTCTACATGGCGTGGTTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGa                                                                                 <  6:176664/149‑1 (MQ=255)
                                                                 aGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAgg                                                                               <  4:201438/149‑1 (MQ=255)
                                                                    cGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTg                                                                            <  1:151663/149‑1 (MQ=255)
                                                                      aGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTggg                                                                          <  5:135963/149‑1 (MQ=255)
                                                                        gTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCTTGGTCTTCAGCCAGTACAGTCGCGCCGAGCGCTCGTCCAGGCCGTCTTCGAGGGTgggg                                                                         >  3:174621/1‑148 (MQ=255)
                                                                           ctctACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATTATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCATCCAGGACAGTCGCGCCGCGCGCGCGTGCAGGCCGTCTTCGAGGGTGGGGTAct                                                                     >  8:217434/1‑148 (MQ=255)
                                                                               aCATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAgtgt                                                                 >  4:282126/1‑149 (MQ=255)
                                                                               aCATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAgtgt                                                                 >  8:28094/1‑149 (MQ=255)
                                                                                   ggCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTGGAGGGTGGGGTACAgt                                                                   >  5:66113/1‑143 (MQ=255)
                                                                                   ggCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAgt                                                                   <  6:66113/143‑1 (MQ=255)
                                                                                    gCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTccc                                                            <  7:28094/149‑1 (MQ=255)
                                                                                          ttCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACg                                                      >  1:195840/1‑149 (MQ=255)
                                                                                             gTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCaa                                                   >  4:75049/1‑149 (MQ=255)
                                                                                                    cgcAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAg                                            >  7:149293/1‑149 (MQ=255)
                                                                                                      cAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCt                                          >  5:296208/1‑149 (MQ=255)
                                                                                                           cGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCg                                     <  7:207435/149‑1 (MQ=255)
                                                                                                            gCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCgg                                    >  2:69379/1‑149 (MQ=255)
                                                                                                              cAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTc                                  <  3:282126/149‑1 (MQ=255)
                                                                                                                        ccccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGa                           <  2:195840/146‑1 (MQ=255)
                                                                                                                        ccccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGa                           <  6:43063/146‑1 (MQ=255)
                                                                                                                        ccccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGc                        >  4:168443/1‑149 (MQ=255)
                                                                                                                        ccccTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGc                        <  4:14606/149‑1 (MQ=255)
                                                                                                                              gATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAg                   <  3:168443/148‑1 (MQ=255)
                                                                                                                                tAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGATGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGc                               <  2:145820/134‑1 (MQ=255)
                                                                                                                                tAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGc                               >  1:145820/1‑134 (MQ=255)
                                                                                                                                        aaGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGACCGAACAGCTcctg                          <  6:44712/131‑4 (MQ=255)
                                                                                                                                        aaGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGACCGAACAGCTcctg                          >  5:44712/1‑128 (MQ=255)
                                                                                                                                         gctcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTc        >  3:80035/4‑148 (MQ=255)
                                                                                                                                           tcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAgt                                                                   <  7:163953/86‑1 (MQ=255)
                                                                                                                                           tcACCGACATTGCCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAgt                                                                   >  8:163953/2‑87 (MQ=255)
                                                                                                                                           ccACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCg       >  3:65448/1‑147 (MQ=255)
                                                                                                                                           ccACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGaa     <  3:75049/149‑1 (MQ=255)
                                                                                                                                            cACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAg    >  3:133633/1‑149 (MQ=255)
                                                                                                                                              ccGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAgtg  <  6:296208/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CTGGAAGATCTCGGCGTTGATCTTGCCTTCGTTCAGTGCTGTCTGCATCGCCTTGGGGTTTTGCAAGGCGAGGTTCTCTACATGGCGTGATTCGTCATAGCGCAACACGCCAAAATTCTGCCCCTGGATAATGGGCAAGCCACCGGCAATGTCGGTGAGCGTGGTCTTCAGCCAGGACAGTCGCGCCGCGCGCTCGTCCAGGCCGTCTTCGAGGGTGGGGTACAGTGTTTCCCAGAACGCAACCAGGAGCTGTTCGGTCAGCTCGAGGCCAAAGGTGATGCCTTCGAAGTG  >  NC_002947/3477402‑3477692

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