Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,944 C→G 8.9% V118V (GTC→GTG)  PP_4304 → VIC family cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,9440CG 8.9% 192.1 / 6.2 64V118V (GTC→GTG) PP_4304VIC family cation transporter
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/4);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACGGTGTTCGGC  >  NC_002947/4895797‑4896088
                                                                                                                                                   |                                                                                                                                                
ggCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCAt                                                                                                                                                >  1:1411070/1‑150 (MQ=255)
         cATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATGCGGATGACCCGGAtgat                                                                                                                                       <  2:127717/150‑3 (MQ=255)
         cATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGATGACCCGGAtgat                                                                                                                                       <  2:127726/150‑3 (MQ=255)
          aTCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGAtgctg                                                                                                                                      >  2:1219612/1‑150 (MQ=255)
          aTCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGATGACCCGGatgatg                                                                                                                                      <  2:587263/150‑1 (MQ=255)
          aTCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCAGCATCCGGATGACCCGGatgatg                                                                                                                                      >  1:269013/1‑150 (MQ=255)
              cccGACTGTCCTGGTCGGCCAGGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGa                                                                                                                                  <  2:660908/150‑1 (MQ=255)
               ccGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGGCGCCCAGTACCTGCTCATCATCCGGATGACCCGGatgatgagatc                                                                                                                                 >  1:513460/1‑145 (MQ=255)
                  ccTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCg                                                                                                                                    <  1:1411917/144‑1 (MQ=255)
                  ccTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCg                                                                                                                                    >  2:1411917/1‑144 (MQ=255)
                    tGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTcc                                                                                                                            <  1:892492/150‑1 (MQ=255)
                      tACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCgt                                                                                                                          <  1:223992/150‑1 (MQ=255)
                         tGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTa                                                                                                                       <  1:1219612/150‑1 (MQ=255)
                           ttCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACt                                                                                                                     >  1:204374/1‑150 (MQ=255)
                              gCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCGAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCaa                                                                                                                  <  2:1131998/150‑1 (MQ=255)
                              gCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCaa                                                                                                                  >  1:746866/1‑150 (MQ=255)
                                ccAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATACTCGCGTCATGCTATAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGc                                                                                                                >  1:530630/1‑150 (MQ=255)
                                   aGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCa                                                                                                             <  2:572831/150‑1 (MQ=255)
                                   aGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCa                                                                                                             <  2:1220626/150‑1 (MQ=255)
                                      cGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTAc                                                                                                                      >  2:816478/1‑138 (MQ=255)
                                      cGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTAc                                                                                                                      <  1:816478/138‑1 (MQ=255)
                                              tATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGAtg                                                                                                                                         >  1:619389/1‑111 (MQ=255)
                                              tATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGAtg                                                                                                                                         <  2:619389/111‑1 (MQ=255)
                                                tGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCaa                                                                                                >  2:770263/1‑150 (MQ=255)
                                                tGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCaa                                                                                                <  1:1308421/150‑1 (MQ=255)
                                                 gCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAg                                                                                               >  1:1088070/1‑150 (MQ=255)
                                                  tcTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGACCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGc                                                                                              <  1:113912/149‑1 (MQ=255)
                                                  ccTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGc                                                                                              >  2:238242/1‑150 (MQ=255)
                                                  ccTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGc                                                                                              >  1:1070049/1‑150 (MQ=255)
                                                    ttCAGTTTCTACGGCCTGGTGGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGAGCATCATGCGTGGACGCAAGCTCAGCCCGTACCTCAAGCaa                                                                                            >  2:113912/1‑149 (MQ=255)
                                                    ttCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGatgatgagcag                                                                                                                                 >  1:919697/1‑108 (MQ=255)
                                                         tttCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCa                                                                                       >  1:886363/1‑150 (MQ=255)
                                                              aCGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTAcc                                                                                  <  2:1088070/150‑1 (MQ=255)
                                                                  ccTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTcc                                                                                                                            >  1:70900/1‑104 (MQ=255)
                                                                  ccTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTcc                                                                                                                            <  2:70900/104‑1 (MQ=255)
                                                                    tGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATgg                                                                            >  2:611264/1‑150 (MQ=255)
                                                                      gTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCg                                                                                                                           >  2:622586/1‑101 (MQ=255)
                                                                      gTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCg                                                                                                                           <  1:622586/101‑1 (MQ=255)
                                                                       tCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTGATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAgg                                                                         <  2:886363/150‑1 (MQ=255)
                                                                       tCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAgg                                                                         >  1:764820/1‑150 (MQ=255)
                                                                                tGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGcgcg                                                                >  2:1215921/1‑150 (MQ=255)
                                                                                   cGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCTgca                                                             <  2:764820/150‑1 (MQ=255)
                                                                                   cGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGgca                                                             <  1:355612/150‑1 (MQ=255)
                                                                                    gATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGgcag                                                            >  1:883309/1‑150 (MQ=255)
                                                                                        gTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCgt                                                                                                                          <  2:814877/84‑1 (MQ=255)
                                                                                        gTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCgt                                                                                                                          >  1:814877/1‑84 (MQ=255)
                                                                                         tGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCagcaagc                                                       >  1:926453/1‑150 (MQ=255)
                                                                                         tGCCGGGGACCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCagcaagc                                                       <  2:530630/150‑1 (MQ=255)
                                                                                             ggggATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCagaa                                                   <  2:1070049/150‑1 (MQ=255)
                                                                                                    cTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAg                                                                          <  1:813351/120‑1 (MQ=255)
                                                                                                    cTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAg                                                                          >  2:813351/1‑120 (MQ=255)
                                                                                                    cTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGatcatc                                            <  2:926453/150‑1 (MQ=255)
                                                                                                    cTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGatcatc                                            <  1:238242/150‑1 (MQ=255)
                                                                                                           tGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCATCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCCGCAAGCAGAAGATCATCGTGTTcc                                     >  2:649467/1‑150 (MQ=255)
                                                                                                              actacAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTgg                                  <  2:1411070/150‑1 (MQ=255)
                                                                                                                  cAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGAc                              <  1:451548/150‑1 (MQ=255)
                                                                                                                                     cTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACg           <  1:770263/150‑1 (MQ=255)
                                                                                                                                        atcatcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATg                                                                             <  1:671549/81‑1 (MQ=255)
                                                                                                                                        atcatcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATg                                                                             >  2:671549/1‑81 (MQ=255)
                                                                                                                                        atcatcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACGgtg        >  1:1208004/1‑150 (MQ=255)
                                                                                                                                         tcatcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCagcaagca                                                      >  1:1143105/1‑103 (MQ=255)
                                                                                                                                         tcatcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCagcaagca                                                      <  2:1143105/103‑1 (MQ=255)
                                                                                                                                            tcCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACGGTGTTCg    >  1:136369/1‑150 (MQ=255)
                                                                                                                                              cGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACGGTGTTCGGc  <  2:883309/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTACTACAGCGACGCCCAGTACCTGCTCATCATCCGGGTCATCCGGATGCTGCGGATCTTCCGTGTACTCAAGCTCAGCCCGTACCTCAAGCAGGCCCATTACCTGATGGAGGCGCTGCGCGGCAGCAAGCAGAAGATCATCGTGTTCCTGGTGACGGTATCGACGCTGGTAACGGTGTTCGGC  >  NC_002947/4895797‑4896088

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