Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,631,143 T→A 7.1% L203Q (CTG→CAG)  lon‑II → DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,631,1430TA 7.1% 185.6 / 3.1 66L203Q (CTG→CAG) lon‑IIDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base A (3/2);  total (34/32)
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.60e-01

GCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAA  >  NC_002947/2631006‑2631290
                                                                                                                                         |                                                                                                                                                    
gCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGaaa                                                                                                                                           <  6:351457/149‑1 (MQ=255)
 ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGattt                                                                                                                                          <  4:90209/149‑4 (MQ=255)
 ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAg                                                                                                                                                  >  4:294470/1‑141 (MQ=255)
 ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAg                                                                                                                                                  <  3:294470/141‑1 (MQ=255)
   aGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgc                                                                                                                                        <  3:116816/149‑1 (MQ=255)
    gggCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAgcg                                                                                                                                       >  5:50355/1‑149 (MQ=255)
                gACACCATGGCCGCGCACATGGCCCTGACGATCCAGCAGAACCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCg                                                                                                                           <  8:345880/149‑1 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCggg                                                                                                                         >  2:82157/1‑146 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCggg                                                                                                                         <  1:82157/146‑1 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                       >  2:277188/1‑148 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                       >  2:26771/1‑148 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                       >  5:317927/1‑148 (MQ=255)
                     cATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTAGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACATCCTGCAGGACGAAAAGCGCAATCGCGGCCgggac                                                                                                                       >  7:332214/1‑146 (MQ=255)
                      ctGGCCGCGCACATGGCCCTGAAGTTCGAGCAGAAGCAGGTAATTCTTGAACTCGTCGACTTGACTACCCGTGTCGAGCCCGTGCTCGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaa                                                                                                                     <  5:290360/148‑1 (MQ=255)
                      aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttcg                                                                                                                                        <  2:243108/130‑6 (MQ=255)
                      aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttcg                                                                                                                                        >  1:243108/1‑125 (MQ=255)
                      aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                       >  7:214990/1‑147 (MQ=255)
                       tGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTCTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaag                                                                                                                    <  1:124732/149‑1 (MQ=255)
                       gactagcgcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaag                                                                                                                    <  7:204309/144‑1 (MQ=255)
                         gCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaaga                                                                                                                   <  4:43957/148‑1 (MQ=255)
                          ccGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaagaag                                                                                                                 <  1:339813/149‑1 (MQ=255)
                             cgcACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAg                                                                                                              <  7:7479/149‑1 (MQ=255)
                                acaTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAACTCGTCGACTTGTCTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAATTCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGTCCGGGTCAAGAAGCAGATg                                                                                                           <  6:96336/149‑1 (MQ=255)
                                  aTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCGGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGa                                                                                                         <  2:297439/149‑1 (MQ=255)
                                     gCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAgcg                                                                                                      >  2:86646/1‑149 (MQ=255)
                                       ccTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAgcg                                                                                                      >  1:52699/1‑147 (MQ=255)
                                        cTGAAGATCAAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGcag                                                                                                   >  2:118494/1‑149 (MQ=255)
                                             gATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGcagccagc                                                                                              >  7:138897/1‑149 (MQ=255)
                                              aTCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGCTCAAGAAGCAGATGGAGCGCAGCCAgcg                                                                                             >  6:84001/1‑149 (MQ=255)
                                                   gCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGt                                                                                        <  1:118494/149‑1 (MQ=255)
                                                      gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGca                                                                                                    <  7:142377/134‑1 (MQ=255)
                                                      gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGca                                                                                                    >  8:142377/1‑134 (MQ=255)
                                                      gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtac                                                                                      >  7:137149/1‑148 (MQ=255)
                                                      gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCGAGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtac                                                                                      >  1:103293/1‑148 (MQ=255)
                                                        aGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtactac                                                                                   <  5:211464/149‑1 (MQ=255)
                                                              aaaTTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAAt                                                                             <  7:218814/149‑1 (MQ=255)
                                                               aaTTCTTGAAATCCTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCTCGTCGATAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAAt                                                                             <  2:130508/148‑1 (MQ=255)
                                                                      gAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCa                                                                        <  2:52699/146‑1 (MQ=255)
                                                                            gtcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAgg                                                                <  6:50355/148‑1 (MQ=255)
                                                                            gtcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAgg                                                                <  3:253602/148‑1 (MQ=255)
                                                                              cgACCTGCCTACCCGTGTCGAGCCCGTGCGGGCTCTGCTGGATGCCGAACTCCACCTCCTCCAGGTCGAAAAGCGCATTCGCGGCCGGGGCAAGAAGCAGCTGGCGCGCAGCCCGCGCGAGTACTACCTGACTGAGCAGATGAAGGCCa                                                             <  8:268865/149‑1 (MQ=255)
                                                                                  ttGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCAtt                                                           >  5:191114/1‑147 (MQ=255)
                                                                                    gACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCag                                                        >  4:96313/1‑148 (MQ=255)
                                                                                    gACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCag                                                        >  1:350921/1‑148 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCATATGGAGCGCAGCCAGCGCGATTACTACCTGAATGAGCAGATGAAg                                                                 <  1:145434/133‑1 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAg                                                                 >  2:145434/1‑133 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                                >  8:227256/1‑149 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                                <  1:26771/149‑1 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                                >  6:320031/1‑149 (MQ=255)
                                                                                          ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCAATTCAGAAAGAAc                                                >  7:270499/1‑149 (MQ=255)
                                                                                            gtgtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACtt                                              >  3:352318/1‑149 (MQ=255)
                                                                                            gtgtCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCAGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAATAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACtt                                              >  7:271225/1‑149 (MQ=255)
                                                                                               tCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAATAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTgg                                            <  1:277188/148‑1 (MQ=255)
                                                                                                   gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                       <  8:214990/149‑1 (MQ=255)
                                                                                                   gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                       >  2:288676/1‑149 (MQ=255)
                                                                                                       gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                       <  4:139135/145‑1 (MQ=255)
                                                                                                       gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                       >  3:139135/1‑145 (MQ=255)
                                                                                                        tGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGAc                                  >  6:140674/1‑149 (MQ=255)
                                                                                                        tGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGAc                                  >  3:265877/1‑149 (MQ=255)
                                                                                                             gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAgg                             >  6:70880/1‑149 (MQ=255)
                                                                                                             gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAg                              <  5:140674/148‑1 (MQ=255)
                                                                                                               cctgcagGAGGTCGATTTCGGCATCCAGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGc                            >  3:113430/2‑148 (MQ=255)
                                                                                                                       aTGCCGAAATCGACCTCCTGCAGGGCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGaa                   <  2:328942/149‑1 (MQ=255)
                                                                                                                           cGAAATCGCCCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc                >  3:246159/1‑148 (MQ=255)
                                                                                                                            gAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGAGATGGTGACGAAGGCCACAATGAAGtag               >  4:46840/1‑146 (MQ=255)
                                                                                                                               aTCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAaga           <  5:70880/149‑1 (MQ=255)
                                                                                                                                    cctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCGGATGAAGGCCATTCAGAAAGAACTGGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTg       >  7:358741/1‑148 (MQ=255)
                                                                                                                                     ctcctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGa      >  3:111285/1‑148 (MQ=255)
                                                                                                                                        ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGATACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaa   <  1:86646/149‑1 (MQ=255)
                                                                                                                                        ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaaa  <  5:320031/149‑1 (MQ=255)
                                                                                                                                        ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaaa  <  2:103293/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
GCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAA  >  NC_002947/2631006‑2631290

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: