Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,388,994 G→C 8.0% L583V (CTG→GTG)  PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,388,9940GC 8.0% 195.5 / 5.1 61L583V (CTG→GTG) PP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base C (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGA  >  NC_002947/4388846‑4389141
                                                                                                                                                    |                                                                                                                                                   
tGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACGTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCgcag                                                                                                                                                     >  6:88009/1‑149 (MQ=255)
  aCTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCgcag                                                                                                                                                     >  6:307902/1‑147 (MQ=255)
      gTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAgg                                                                                                                                                 <  3:238637/147‑1 (MQ=255)
           cACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACGTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGccc                                                                                                                                           <  5:88009/148‑1 (MQ=255)
                 cTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCACCGGGCGGCCTGCTGCgg                                                                                                                                    <  6:352303/149‑1 (MQ=255)
                   cGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCggg                                                                                                                                   >  2:53287/1‑148 (MQ=255)
                    gACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCggg                                                                                                                                   >  3:113145/1‑147 (MQ=255)
                         gCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                            <  6:238718/149‑1 (MQ=255)
                            ggccaggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTc                                                                                                                         <  2:47235/149‑1 (MQ=255)
                               caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCTGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTcgcc                                                                                                                      <  5:307902/149‑1 (MQ=255)
                               caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                            >  6:114822/1‑143 (MQ=255)
                               caggcGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCggg                                                                                                                            <  5:114822/143‑1 (MQ=255)
                                    ggcgAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCcgcacg                                                                                                                 >  5:90901/1‑149 (MQ=255)
                                          ttATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGgcc                                                                                                                                               <  2:349292/113‑1 (MQ=255)
                                           tATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAgc                                                                                                          >  2:289243/1‑149 (MQ=255)
                                             tCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCgg                                                                                                                                    >  8:315737/1‑121 (MQ=255)
                                             tCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCgg                                                                                                                                    <  7:315737/121‑1 (MQ=255)
                                              cGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGcc                                                                                                       >  3:305087/1‑149 (MQ=255)
                                              cGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGcc                                                                                                       >  7:194372/1‑149 (MQ=255)
                                                aaGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTAGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCt                                                                                                      >  5:63125/1‑148 (MQ=255)
                                                   gCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAgg                                                                                                  >  5:57659/1‑149 (MQ=255)
                                                        ccGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                             >  3:169837/1‑149 (MQ=255)
                                                          gCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTAcc                                                                                           <  4:113145/149‑1 (MQ=255)
                                                              cTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCa                                                                                        >  5:99804/1‑148 (MQ=255)
                                                                   cTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCt                                                                                  <  6:90901/149‑1 (MQ=255)
                                                                    ttCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCtt                                                                                 >  4:251455/1‑149 (MQ=255)
                                                                         gctgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcg                                                                            >  8:282126/1‑149 (MQ=255)
                                                                           tgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                             >  3:35232/1‑130 (MQ=255)
                                                                           tgctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGCGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTa                                                                                             <  4:35232/130‑1 (MQ=255)
                                                                             ctgctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCTCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcgcg                                                                          <  1:289243/147‑1 (MQ=255)
                                                                                ctTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgcc                                                                      <  7:282126/148‑1 (MQ=255)
                                                                                     cacccaccTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgccgc                                                                    <  8:269312/145‑1 (MQ=255)
                                                                                     cacccaccTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGcgccgc                                                                    >  7:269312/1‑145 (MQ=255)
                                                                                            cTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGc                                                                                   <  2:36263/123‑1 (MQ=255)
                                                                                            cTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGc                                                                                   >  1:36263/1‑123 (MQ=255)
                                                                                               cgcCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTgcg                                                      <  8:194372/149‑1 (MQ=255)
                                                                                                  cGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGa                                                   >  3:5210/1‑149 (MQ=255)
                                                                                                   gATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                  <  6:57659/149‑1 (MQ=255)
                                                                                                       gACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTg                                               >  4:140684/1‑148 (MQ=255)
                                                                                                        aCGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGc                                              <  3:251455/148‑1 (MQ=255)
                                                                                                        aCGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCa                                             >  5:163045/1‑149 (MQ=255)
                                                                                                         cGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCAc                                            >  7:196062/1‑149 (MQ=255)
                                                                                                          gCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACg                                           <  8:196062/149‑1 (MQ=255)
                                                                                                           cAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACgg                                          <  8:205325/149‑1 (MQ=255)
                                                                                                                aaGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCAtt                                      <  4:169837/148‑1 (MQ=255)
                                                                                                                 aGCCCCCGGCTCGGGCAGTAGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTgg                                    <  6:63125/149‑1 (MQ=255)
                                                                                                                   cccccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGcc                                  >  8:173285/1‑149 (MQ=255)
                                                                                                                    ccccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCt                                 <  6:163045/149‑1 (MQ=255)
                                                                                                                      ccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTg                                >  6:269167/1‑148 (MQ=255)
                                                                                                                      ccGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGc                               <  4:5210/149‑1 (MQ=255)
                                                                                                                          cTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCa                           <  4:305087/149‑1 (MQ=255)
                                                                                                                            cGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCAcc                            >  3:119704/1‑146 (MQ=255)
                                                                                                                            cGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCAcc                            <  4:119704/146‑1 (MQ=255)
                                                                                                                                 aGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTTGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGcgcc                    >  1:209565/1‑149 (MQ=255)
                                                                                                                                  gTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGcc                                                           >  5:53708/1‑109 (MQ=255)
                                                                                                                                  gTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGcc                                                           <  6:53708/109‑1 (MQ=255)
                                                                                                                                      gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTcgcgc                                                                         >  6:310072/1‑91 (MQ=255)
                                                                                                                                      gCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTgc               >  2:290263/1‑149 (MQ=255)
                                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                  >  3:47038/5‑111 (MQ=38)
                                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                  <  4:47038/107‑1 (MQ=38)
                                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                  <  2:253438/107‑1 (MQ=38)
                                                                                                                                         cgctgCCCGCACCGGGCGGCCTGCTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGAc                                                  >  1:253438/5‑111 (MQ=38)
                                                                                                                                          gACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTgccgcc           >  4:293463/1‑149 (MQ=255)
                                                                                                                                                 gcagcaGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACt    <  6:99804/149‑1 (MQ=255)
                                                                                                                                                   agcagGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGa  <  3:293463/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGACTGGTTTTCACCGACTCGACACGCAGGCCAGGCGGCGAGTTATCGAAGGCCACCCGCACCTCCCCTTCCAGCTGCTGCTTGCCACCCACCTGCGCCGATTGACGCAGTAAAGCCCCCGGCTCGGGCAGTTGGCCGGACGCCCGCAGCAGGCCGCCCGGTGCGGGCAGCGGGGTCGCCGCACGCAACAGCGCCTCAGGGCTTACCCCAAAGGCTTGCTCGCGCGCCGCCTGGACGCCCTGCGCGACCTGCACGGCATTGGCCTGCACCAGGGCGCCGCTGCCGCCGGCCACTGA  >  NC_002947/4388846‑4389141

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