Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,565 C→T 8.2% A368V (GCG→GTG)  PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,268,5650CT 8.2% 163.0 / 4.3 60A368V (GCG→GTG) PP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (3/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGA  >  NC_002947/5268421‑5268705
                                                                                                                                                |                                                                                                                                             
ccGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCa                                                                                                                                          <  2:1039037/150‑1 (MQ=255)
  ggggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGc                                                                                                                                        >  1:931938/1‑150 (MQ=255)
   gggTCAAGGGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCTTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCg                                                                                                                                       <  1:274436/150‑1 (MQ=255)
   gggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCg                                                                                                                                       >  1:317837/1‑150 (MQ=255)
     gTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGt                                                                                                                                     >  1:945495/1‑150 (MQ=255)
        aaGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGTGGCCAGCGGTAT‑c                                                                                                                                 >  1:1145943/1‑150 (MQ=255)
              cGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGAGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCAcc                                                                                                                           <  2:1016161/150‑1 (MQ=255)
                  aCGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGt                                                                                                                       <  2:393638/150‑1 (MQ=255)
                    gtggtgGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTcc                                                                                                                     >  2:824311/1‑150 (MQ=255)
                            tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGc                                                                                                                                        >  1:572055/1‑124 (MQ=255)
                            tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGc                                                                                                                                        <  2:572055/124‑1 (MQ=255)
                                  cggcgTGATTCTGGTGCTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCATGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACaa                                                                                                       <  1:655822/150‑1 (MQ=255)
                                   ggcgTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCAACAAGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAg                                                                                                      <  1:1613591/150‑1 (MQ=255)
                                      gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                   >  2:148867/1‑150 (MQ=255)
                                      gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                                                   <  1:633507/150‑1 (MQ=255)
                                         aTTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGaacaac                                                                                                <  1:148867/150‑1 (MQ=255)
                                          ttCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                                                                                                              >  2:848751/1‑104 (MQ=255)
                                          ttCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                                                                                                              <  1:848751/104‑1 (MQ=255)
                                          ttCTGGTGGTACTCGGCTTGCTCCCAGTCATGGACCGGGTGATTGCCGCGGTTCCTACTCCGGTACGGGACCACGCTGGTATCGTGCTGTTCGGGACCGTGGCGGCCAGCGGTAT‑CCGCACGCTGTCCAAGGTCGGCTACAAGCACAACg                                                                                               >  2:220187/1‑150 (MQ=255)
                                             tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcgg                                                                                                                                             <  2:714866/102‑1 (MQ=255)
                                             tggtggTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcgg                                                                                                                                             >  1:714866/1‑102 (MQ=255)
                                                    aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑c                                                                                                                                 >  1:174464/1‑106 (MQ=255)
                                                    aCTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑c                                                                                                                                 <  2:174464/106‑1 (MQ=255)
                                                      tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                          >  2:649185/1‑143 (MQ=255)
                                                      tCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                                          <  1:649185/143‑1 (MQ=255)
                                                        ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                 >  2:1332464/1‑150 (MQ=255)
                                                        ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                                                 >  2:90832/1‑150 (MQ=255)
                                                         gCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcg                                                                                >  1:914885/1‑150 (MQ=255)
                                                          cTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCAtcgt                                                                               <  1:90832/150‑1 (MQ=255)
                                                                       aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCa                                                                                                                                          <  1:739270/79‑1 (MQ=255)
                                                                       aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCa                                                                                                                                          >  2:739270/1‑79 (MQ=255)
                                                                       aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACATGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTg                                                                  <  1:1332464/150‑1 (MQ=255)
                                                                       aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                                                      >  2:1555798/1‑130 (MQ=255)
                                                                       aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                                                      <  1:1555798/130‑1 (MQ=255)
                                                                        tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                 <  2:1386945/150‑1 (MQ=255)
                                                                        tGGGCCGGGTGATTGCCGCCGTGCCTACCACGGTCATGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                                                 <  2:1088800/150‑1 (MQ=255)
                                                                            ccGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTCTCGGGCGGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCtt                                                             >  1:1276437/1‑150 (MQ=255)
                                                                             cGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTc                                                            <  2:914885/150‑1 (MQ=255)
                                                                               ggTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCgg                                                          >  1:838692/1‑150 (MQ=255)
                                                                                    ttGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTc                                                            <  2:105180/143‑1 (MQ=255)
                                                                                    ttGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTc                                                            >  1:105180/1‑143 (MQ=255)
                                                                                     tGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGAt                                                    >  1:764223/1‑150 (MQ=255)
                                                                                      gccgccGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATc                                                   >  2:1560433/1‑150 (MQ=255)
                                                                                       ccgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATcc                                                  >  2:157914/1‑150 (MQ=255)
                                                                                        cgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATccc                                                 <  2:1146572/150‑1 (MQ=255)
                                                                                        cgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATccc                                                 <  2:1146701/150‑1 (MQ=255)
                                                                                                   ccccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGTCGGCCAGCCGCAT‑CAGCGCCCTGTTCGAGCGCATCGACAAGGACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCCGCATGATCCCGATTTCCGCAc                                      >  1:1019065/1‑150 (MQ=255)
                                                                                                   ccccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCAc                                      <  2:764223/150‑1 (MQ=255)
                                                                                                     ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCGCCCTGATAG‑CAGAGATGTCCCTCACCCTGTCCAAATTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCg                                    >  1:692091/1‑150 (MQ=255)
                                                                                                             gggCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTa                            >  2:1385601/1‑150 (MQ=255)
                                                                                                              ggcggcGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTac                           <  1:1385601/150‑1 (MQ=255)
                                                                                                                cggcgCCGGTATCGTGCTGTTCGGCACCGTGGTGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTacca                         <  2:1145943/150‑1 (MQ=255)
                                                                                                                 ggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTaccac                        <  2:317837/150‑1 (MQ=255)
                                                                                                                   cgcCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTAccacca                      <  1:824311/150‑1 (MQ=255)
                                                                                                                    gcCGGTATCGTGCTGTTCGGCACCGTGGTGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACcaccac                     >  1:23385/1‑150 (MQ=255)
                                                                                                                    gcCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACcaccac                     <  1:1480820/150‑1 (MQ=255)
                                                                                                                      cGGTATCGTGCTGTTCGTCATCGTGGCGGCCAGCGGTAT‑CCGCAGACTGTCCAACGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACtt                   <  2:692091/150‑1 (MQ=255)
                                                                                                                      cGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACtt                   >  2:200470/1‑150 (MQ=255)
                                                                                                                       ggTATCGTGCTGTTCGGCACCGTGGTGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTc                  <  2:23385/150‑1 (MQ=255)
                                                                                                                                 tGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGt        <  1:1498594/150‑1 (MQ=255)
                                                                                                                                   ttCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTc      <  1:157914/150‑1 (MQ=255)
                                                                                                                                      ggCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgag   >  1:60435/1‑150 (MQ=255)
                                                                                                                                       gcgccGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga  <  1:1136842/147‑1 (MQ=255)
                                                                                                                                       gCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCgaga  <  1:1560433/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAT‑CCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCACCACTTCCCGAACTGGTTCGAGA  >  NC_002947/5268421‑5268705

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