Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,388,833 G→C 8.6% T636T (ACC→ACG) ‡ PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,388,8330GC 8.6% 208.3 / 10.2 70T636T (ACC→ACG) ‡PP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base C (4/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

GTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388685‑4388979
                                                                                                                                                    |                                                                                                                                                   
gTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATg                                                                                                                                                     <  5:53153/149‑1 (MQ=255)
   ttCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGt                                                                                                                                                   <  6:161388/148‑1 (MQ=255)
      gTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACg                                                                                                                                              >  1:47914/1‑149 (MQ=255)
      gTGCAGCTCGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGACTACCATGGTTCGCCCCCCGACATTCGGCGTCATGGAC‑ACg                                                                                                                                              >  8:224950/1‑149 (MQ=255)
           gCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGCCgaa                                                                                                                                         <  6:375171/149‑3 (MQ=255)
                gTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGAc                                                                                                                                    <  2:182298/149‑1 (MQ=255)
                   cgcCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTg                                                                                                                                  <  5:12265/148‑1 (MQ=255)
                    gccggccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGACTCGCGTCGCAAATCCCGCCATTCACTCATAGGTGCCTCCCATCGTTCTCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTgct                                                                                                                                >  8:10877/1‑147 (MQ=255)
                     ccggccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGtt                                                                                                                               <  1:271088/149‑1 (MQ=255)
                         ccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCa                                                                                                                           <  6:111888/149‑1 (MQ=255)
                         ccgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCa                                                                                                                           <  6:162479/149‑1 (MQ=255)
                          cgACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          >  7:125434/1‑149 (MQ=255)
                           gACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAcc                                                                                                                         >  8:267125/1‑149 (MQ=255)
                             ccGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGa                                                                                                                       >  4:199250/1‑149 (MQ=255)
                              cGGCACGCTGTCGCTGTCAACCCAGAACGGGACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      >  1:35441/1‑149 (MQ=255)
                              cGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      >  2:229868/1‑149 (MQ=255)
                              cGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGAc                                                                                                                      <  3:294879/149‑1 (MQ=255)
                                  aCGCTGTCGCTGTCAACCCAGAACGGCACGGCGCGAAACGACGCCTGGCGTCGCAAATCACGCCATTCACTGATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑TCGCGGGGCTGACTGGTTTTCACCGAGTCGa                                                                                                                  <  7:224950/149‑1 (MQ=255)
                                     ctgtcgctgtcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGacac                                                                                                               <  7:281896/149‑1 (MQ=255)
                                        tcgctgtcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACg                                                                                                              >  1:150436/1‑147 (MQ=255)
                                          gctgtcAACCCAGAACGGAACTCCGCGAAACGACGCCTCGCGTGGCAAATCACGCCATTCACTCTTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcag                                                                                                           <  7:10877/148‑1 (MQ=255)
                                             gtcAACCAAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGAGCCTCCCATCGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGCCTGGTTTTCCACGACTCGACACGcaggc                                                                                                         >  7:143066/1‑147 (MQ=255)
                                              tcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcaggccag                                                                                                      <  2:19277/149‑1 (MQ=255)
                                              tcAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGcaggcca                                                                                                       <  8:156536/148‑1 (MQ=255)
                                                 aCCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          <  6:36614/126‑1 (MQ=255)
                                                 aCCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCAc                                                                                                                          >  5:36614/1‑126 (MQ=255)
                                                   ccAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGGCTGGTTTTCGCCGACTCGACCCGCAGGCCCggcggc                                                                                                 >  8:276063/1‑149 (MQ=255)
                                                    cAGAGCGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGcc                                                                                                                                            <  3:358686/105‑1 (MQ=255)
                                                    cAGAGCGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGACCATC‑ACGcc                                                                                                                                            >  4:358686/1‑105 (MQ=255)
                                                     aGAACGGCCCGCCGCCAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGa                                                                                               <  2:324036/149‑1 (MQ=255)
                                                         cGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTa                                                                                           >  3:18799/1‑149 (MQ=255)
                                                            cACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATc                                                                                         >  8:190140/1‑148 (MQ=255)
                                                               gccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAg                                                                                     <  4:15585/149‑1 (MQ=255)
                                                               gccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAg                                                                                     <  2:47914/149‑1 (MQ=255)
                                                                ccgcGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGCCATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAgg                                                                                    >  5:375065/1‑149 (MQ=255)
                                                                      aaCGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCAccc                                                                              >  4:332386/1‑149 (MQ=255)
                                                                          acgCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCATCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCAc                                                                          <  1:289685/149‑1 (MQ=255)
                                                                            gCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCt                                                                        >  2:386683/1‑149 (MQ=255)
                                                                             ccTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTc                                                                       <  6:19602/149‑1 (MQ=255)
                                                                                      gCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAg                                                              <  2:150436/149‑1 (MQ=255)
                                                                                           tCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgct                                                         >  4:204264/1‑149 (MQ=255)
                                                                                             aCGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgc                                                       <  7:190140/149‑1 (MQ=255)
                                                                                               gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCTTCATTCGGCGTGATGGTC‑AGGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAATGCCAGCCGCACCTCCCCTTCCATCTGCTGCtt                                                     <  2:14330/149‑1 (MQ=255)
                                                                                               gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTCGACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgct                                                      >  1:305142/1‑149 (MQ=255)
                                                                                               gCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAgctgctgcat                                                     >  4:254382/1‑147 (MQ=255)
                                                                                                 cATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGc                                                   >  5:21241/1‑149 (MQ=255)
                                                                                                  aTTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGcc                                                  >  4:367685/1‑149 (MQ=255)
                                                                                                  aTTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGc                                                   >  7:287988/1‑148 (MQ=255)
                                                                                                   ttCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGcca                                                 <  3:204264/149‑1 (MQ=255)
                                                                                                      aCTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGccacc                                               >  6:325617/1‑148 (MQ=255)
                                                                                                          aTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCCCCCACCt                                          >  4:96950/1‑149 (MQ=255)
                                                                                                          aTAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCt                                          <  8:287988/149‑1 (MQ=255)
                                                                                                          aTAGGTGCCTCCCATGGTTAGACGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTTCTTGccacccac                                            <  6:159498/147‑1 (MQ=255)
                                                                                                             ggTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTgc                                        <  4:18799/148‑1 (MQ=255)
                                                                                                                 ccTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCCCCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGa                                   >  2:89819/1‑149 (MQ=255)
                                                                                                                     ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGTCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGAtt                                 <  7:267125/147‑1 (MQ=255)
                                                                                                                     ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGAtt                                 <  3:332386/147‑1 (MQ=255)
                                                                                                                     ccATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGa                               >  8:95579/1‑149 (MQ=255)
                                                                                                                     ccATGGTTCGCCGCCC‑ACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGAc                              >  5:176904/1‑149 (MQ=255)
                                                                                                                      cATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGAc                              >  4:49666/1‑149 (MQ=255)
                                                                                                                          gTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAg                          <  1:229868/149‑1 (MQ=255)
                                                                                                                           ttCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAg                          >  8:137667/1‑148 (MQ=255)
                                                                                                                             cgccgcCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTa                        >  6:250311/1‑148 (MQ=255)
                                                                                                                              gccgccCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTaaa                      <  5:325617/149‑1 (MQ=255)
                                                                                                                                       agcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGc             >  2:131461/5‑149 (MQ=255)
                                                                                                                                        gcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            >  5:78849/4‑149 (MQ=255)
                                                                                                                                        gcccGGCGTGACCATC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            >  3:160214/4‑149 (MQ=255)
                                                                                                                                        aTTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCt            <  3:301552/149‑1 (MQ=255)
                                                                                                                                           cGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCgg         >  6:142497/1‑149 (MQ=255)
                                                                                                                                           cGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCg          <  8:125434/148‑1 (MQ=255)
                                                                                                                                                  aTGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTg  >  6:28400/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTATTCGTGCAGCTGGGTGCGCCGGCCGACCGGCACGCTGTCGCTGTCAACCCAGAACGGCACGCCGCGAAACGACGCCTCGCGTCGCAAATCACGCCATTCACTCATAGGTGCCTCCCATGGTTCGCCGCCCGACATTCGGCGTGATGGTC‑ACGCCGGGCTGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTG  >  NC_002947/4388685‑4388979

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