Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,940 G→C 5.3% D90E (GAC→GAG PP_5039 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,9400GC 5.3% 366.3 / 7.6 112D90E (GAC→GAGPP_5039hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (52/54);  new base C (4/2);  total (56/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 115 total aligned reads displayed.

AGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741791‑5742085
                                                                                                                                                     |                                                                                                                                                 
agcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCgg                                                                                                                                                   <  1:787688/150‑1 (MQ=255)
  caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTc                                                                                                                                                 <  1:221510/150‑1 (MQ=255)
      cGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGccc                                                                                                                                             <  2:973907/150‑1 (MQ=255)
       gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                            >  1:110501/1‑150 (MQ=255)
              cTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgc                                                                                                                                     >  1:1280000/1‑150 (MQ=255)
                    gCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcgg                                                                                                                               >  2:2403820/1‑150 (MQ=255)
                     cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTTGACACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                              <  2:1483747/150‑7 (MQ=255)
                     cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcggc                                                                                                                              <  1:140245/150‑1 (MQ=255)
                         ggTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCgg                                                                                                                                                   >  1:724954/1‑125 (MQ=255)
                         ggTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCgg                                                                                                                                                   <  2:724954/125‑1 (MQ=255)
                             cTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcgcg                                                                                                                      <  2:1280000/150‑1 (MQ=255)
                                    tCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATc                                                                                                               >  1:664627/1‑150 (MQ=255)
                                       aGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcggc                                                                                                                              <  1:169591/132‑1 (MQ=255)
                                       aGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcggc                                                                                                                              >  2:169591/1‑132 (MQ=255)
                                          aCTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                         >  2:681993/1‑150 (MQ=255)
                                          aCTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGt                                                                                                         >  2:681491/1‑150 (MQ=255)
                                           cTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTg                                                                                                        <  1:2143376/150‑1 (MQ=255)
                                           cTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTg                                                                                                        >  2:952018/1‑150 (MQ=255)
                                            ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTgg                                                                                                       <  1:952018/150‑1 (MQ=255)
                                                gCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAg                                                                                                    >  2:2036979/1‑149 (MQ=255)
                                                gCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAg                                                                                                    <  1:2036979/149‑1 (MQ=255)
                                                gCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGc                                                                                                   <  1:2305136/150‑1 (MQ=255)
                                                 cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCt                                                                                                  >  1:2085445/1‑150 (MQ=255)
                                                 cATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGACGTGTCGCTAATCTGCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCt                                                                                                  >  2:64523/1‑150 (MQ=255)
                                                  aTGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTg                                                                                                 <  1:979926/150‑1 (MQ=255)
                                                  aTGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCGGGTCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTg                                                                                                 <  2:2011586/150‑1 (MQ=255)
                                                    gCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGa                                                                                                                 >  1:1808581/1‑132 (MQ=255)
                                                    gCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGa                                                                                                                 <  2:1808581/132‑1 (MQ=255)
                                                       ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            <  2:1901769/150‑1 (MQ=255)
                                                       ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            <  1:681993/150‑1 (MQ=255)
                                                       ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            <  1:681491/150‑1 (MQ=255)
                                                       ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            <  2:110501/150‑1 (MQ=255)
                                                       ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTACAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                            >  2:671473/1‑150 (MQ=255)
                                                         aGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGc                                                                                          <  2:1078342/150‑1 (MQ=255)
                                                                gcagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGACCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                   >  1:1213104/1‑150 (MQ=255)
                                                                   gcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagca                                                                                >  1:356795/1‑150 (MQ=255)
                                                                   gcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcagca                                                                                >  2:339914/1‑150 (MQ=255)
                                                                     aCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                              >  2:431904/1‑150 (MQ=255)
                                                                         acacGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCcca                                                                          >  1:1822186/1‑150 (MQ=255)
                                                                           acGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagc                                                                        <  1:1656539/150‑1 (MQ=255)
                                                                                     tCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCg                                                              >  2:439696/1‑150 (MQ=255)
                                                                                      cGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgg                                                             <  1:1582313/150‑1 (MQ=255)
                                                                                       ggCAGCATCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCACGAAGCCGCGCCCCGTGTCGCGGTAGATCATCCGTGCCATCTCGTCCAGCACACTCAGCAGGCCCAGCCAGCCCttctgg                                                            <  2:419689/150‑6 (MQ=255)
                                                                                       ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                            <  2:1203155/150‑1 (MQ=255)
                                                                                       ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                            <  1:26093/150‑1 (MQ=255)
                                                                                        gcagccGCTGCTAGTGCCAGCCCGTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGACCTCGCCGGCGCGAGGCGTGTCGCGCAACCTCGGGCCCCTCAGCCTGTCGATCAAAATCAGCATGCCCAGCCAGCCCTGCgggg                                                           >  2:1729865/1‑150 (MQ=255)
                                                                                        gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATc                                                                                                               <  1:1138737/98‑1 (MQ=255)
                                                                                        gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                           >  1:599414/1‑150 (MQ=255)
                                                                                        gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                           >  1:587157/1‑150 (MQ=255)
                                                                                         cagcagCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAg                                                                                                    <  1:996064/108‑1 (MQ=255)
                                                                                          agcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGAGCGGTAGAGCGCTGGAGATCATGCGTGGCAGCTGGTCGAGCAAAATAACCAGGCCAACCAAACCCTGCGGGGtt                                                         >  1:760728/1‑150 (MQ=255)
                                                                                           gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        >  2:764760/1‑150 (MQ=255)
                                                                                           gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        >  1:1682535/1‑150 (MQ=255)
                                                                                                  gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGTGGTAGATCATGCGTGGCAGCTGGTCGAGGGTGGGCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCa                                                 >  1:1331473/1‑150 (MQ=255)
                                                                                                  gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGt                                                          <  2:2219152/141‑1 (MQ=255)
                                                                                                  gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCa                                                 <  1:339914/150‑1 (MQ=255)
                                                                                                     gTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                       >  1:1933168/1‑61 (MQ=255)
                                                                                                     gTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                       <  2:1933168/61‑1 (MQ=255)
                                                                                                       cccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCtt                                            <  2:1682535/150‑1 (MQ=255)
                                                                                                        ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTgg                                                                                                >  2:14391/1‑97 (MQ=255)
                                                                                                        ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTgg                                                                                                <  1:14391/97‑1 (MQ=255)
                                                                                                        ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                           >  1:2246500/1‑150 (MQ=255)
                                                                                                        ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                           >  2:428995/1‑150 (MQ=255)
                                                                                                        ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                           >  1:1446358/1‑150 (MQ=255)
                                                                                                         cAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGtt                                                         <  2:1189302/135‑1 (MQ=255)
                                                                                                         cAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGtt                                                         >  1:1189302/1‑135 (MQ=255)
                                                                                                         cAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctgc                                                   <  1:539082/141‑1 (MQ=255)
                                                                                                          aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTcta                                         <  2:2089133/150‑3 (MQ=255)
                                                                                                           gCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCACAATCAGCAGGCCCAGCCAGCCCTTCGGGGTTTGCTGCTACTCTTCGAg                                        <  2:149037/150‑1 (MQ=255)
                                                                                                            cccTTCTGCAGGCCCTGCATGGCGACCCCCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGc                                       <  2:754188/150‑1 (MQ=255)
                                                                                                            cccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                             <  2:1636411/80‑1 (MQ=255)
                                                                                                             ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAGAGCGGCAGGCCCAGCCAGCCCTGCGGGGCTTGCTGCCACTCTTCGAGcc                                      >  1:760457/1‑150 (MQ=255)
                                                                                                             ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgg                                                             <  1:2174450/127‑1 (MQ=255)
                                                                                                              cTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGGCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAACATCACCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGccc                                     <  2:1003285/150‑1 (MQ=255)
                                                                                                              cTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGcac                                     >  1:136106/1‑148 (MQ=255)
                                                                                                               ttCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGACCCAGCCAGCCCTGCTTGGTTTGCTGCCCATCATCGAGACCa                                    >  2:1599689/1‑150 (MQ=255)
                                                                                                                 cTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAcc                                  <  2:740732/150‑1 (MQ=255)
                                                                                                                  tGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAccc                                 <  2:2085445/150‑1 (MQ=255)
                                                                                                                    cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACCCTTCGAGCCCACCCGc                               <  2:640949/150‑1 (MQ=255)
                                                                                                                      ggCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCa                             <  1:1729865/150‑1 (MQ=255)
                                                                                                                        cccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAgg                           <  1:1131233/150‑1 (MQ=255)
                                                                                                                          cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                            <  1:701443/115‑1 (MQ=255)
                                                                                                                            gCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGt                       <  2:2242611/150‑1 (MQ=255)
                                                                                                                               tGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGt                                                          >  2:436318/1‑112 (MQ=255)
                                                                                                                               tGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGt                                                          <  1:436318/112‑1 (MQ=255)
                                                                                                                                  cGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGa                 <  2:2246500/150‑1 (MQ=255)
                                                                                                                                      caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcag                                                                                  <  1:363051/81‑1 (MQ=255)
                                                                                                                                      caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcag                                                                                  >  2:363051/1‑81 (MQ=255)
                                                                                                                                      caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg             >  2:1015926/1‑150 (MQ=255)
                                                                                                                                      caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                         >  2:669403/1‑122 (MQ=255)
                                                                                                                                      caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGa                                         <  1:669403/122‑1 (MQ=255)
                                                                                                                                       accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg            <  2:664627/150‑1 (MQ=255)
                                                                                                                                         cgcgacACGCCGCGCGCCTTCGAGGGCGACCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTc          >  2:2263667/6‑150 (MQ=37)
                                                                                                                                         cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagcca                                                                      >  2:1383074/1‑90 (MQ=255)
                                                                                                                                         cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCccagcca                                                                      <  1:1383074/90‑1 (MQ=255)
                                                                                                                                         cTGGGCACGCCGGACGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTc          >  1:1180951/1‑150 (MQ=255)
                                                                                                                                          gcgacACGCCGCGCGCCTTCGAGGGCGACCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtg                                                       <  1:2287815/100‑1 (MQ=21)
                                                                                                                                          gcgacACGCCGCGCGCCTTCGAGGGCGACCGGCGTGCCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtg                                                       >  2:2287815/5‑104 (MQ=21)
                                                                                                                                                 gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTATATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACa  >  1:1950065/1‑150 (MQ=255)
                                                                                                                                                 gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACa  >  2:1380569/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741791‑5742085

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