Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,573 T→G 6.2% Q20P (CAG→CCG)  PP_3382 ← gluconate 2‑dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,829,5730TG 6.2% 237.3 / 6.6 81Q20P (CAG→CCG) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base T (35/41);  new base G (2/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGC  >  NC_002947/3829426‑3829722
                                                                                                                                                   |                                                                                                                                                     
gtggAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCtgt                                                                                                                                                     >  1:966064/1‑150 (MQ=255)
   gAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCa                                                                                                                                                  >  1:162078/1‑150 (MQ=255)
       ccacCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCg                                                                                                                                              <  2:1151887/149‑1 (MQ=255)
        cacCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCAGTCGCTACTGCCTGTGCcgccgc                                                                                                                                             >  1:236372/1‑150 (MQ=255)
         acCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCa                                                                                                                                            >  2:149812/1‑150 (MQ=255)
            gTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGc                                                                                                                                                >  1:1230170/1‑143 (MQ=255)
            gTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGc                                                                                                                                                <  2:1230170/143‑1 (MQ=255)
              gCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAGgc                                                                                                                                       <  1:461162/150‑1 (MQ=255)
                cAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGACGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                     >  2:802857/1‑150 (MQ=255)
                   tCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCa                                                                                                                                  >  2:380822/1‑150 (MQ=255)
                   tCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCa                                                                                                                                  >  1:345354/1‑150 (MQ=255)
                        gTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCAcc                                                                                                                             >  2:871418/1‑150 (MQ=255)
                           tCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAgg                                                                                                                          <  1:622491/150‑1 (MQ=255)
                           tCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAgg                                                                                                                          <  2:720490/150‑1 (MQ=255)
                            ccATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAggg                                                                                                                         >  1:1013172/1‑150 (MQ=255)
                              aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGc                                                                                                                                             >  2:857608/1‑128 (MQ=255)
                              aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGc                                                                                                                                             <  1:857608/128‑1 (MQ=255)
                              aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTAAGTGCTCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTg                                                                                                                       >  1:430652/1‑150 (MQ=255)
                               tCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCNTNATTCGCNACGGCCTGTGCANCCGCACCGGCGCCCAGCACCAGGGNgg                                                                                                                      <  1:1042174/150‑1 (MQ=255)
                                cGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                     <  2:102361/150‑1 (MQ=255)
                                   cAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                                  >  2:387974/1‑150 (MQ=255)
                                      cccGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCaa                                                                                                               <  1:1065283/150‑1 (MQ=255)
                                             gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACa                                                                                                             >  1:578477/1‑145 (MQ=255)
                                             gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                        >  2:604303/1‑150 (MQ=255)
                                             gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                        <  1:149812/150‑1 (MQ=255)
                                             gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGACGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                        <  1:802857/150‑1 (MQ=255)
                                             gCAAAGGGTTTGCCGCCTTTGGCGGGATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACa                                                                                                             <  2:578477/145‑1 (MQ=255)
                                               aaaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                        >  1:878403/1‑116 (MQ=255)
                                               aaaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                        <  2:878403/116‑1 (MQ=255)
                                                aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGTGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                     >  2:1262600/1‑150 (MQ=255)
                                                aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                     <  1:899430/150‑1 (MQ=255)
                                                aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                     >  2:1412365/1‑150 (MQ=255)
                                                aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                     >  2:31610/1‑150 (MQ=255)
                                                  gggTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCAt                                                                                                   >  2:969666/1‑150 (MQ=255)
                                                       tGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCa                                                                                              <  1:1407737/150‑1 (MQ=255)
                                                        gccgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCAt                                                                                             <  1:871418/150‑1 (MQ=255)
                                                          cgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                     >  2:412650/1‑124 (MQ=255)
                                                          cgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                     <  1:412650/124‑1 (MQ=255)
                                                           gccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGc                                                                                          <  1:31610/150‑1 (MQ=255)
                                                             cTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGTGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGctc                                                                                        <  1:1262600/150‑1 (MQ=255)
                                                                  gCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCaa                                                                                   <  2:285791/150‑1 (MQ=255)
                                                                         ggCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGa                                                                            <  2:1013172/150‑1 (MQ=255)
                                                                             ggCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                     <  1:1196222/89‑1 (MQ=255)
                                                                             ggCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                     >  2:1196222/1‑89 (MQ=255)
                                                                                 aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCAcag                                                                                                                                         <  1:1404559/81‑3 (MQ=255)
                                                                                 aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCAcag                                                                                                                                         >  2:1404559/1‑79 (MQ=255)
                                                                                    cAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCg                                                                 <  2:345354/150‑1 (MQ=255)
                                                                                      gccaacGGCGCGGGCCAGCTACTCGCCTAGGAGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACTGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCggg                                                               <  2:263130/145‑1 (MQ=255)
                                                                                         aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCaa                                                                        >  2:1525472/1‑138 (MQ=255)
                                                                                         aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCaa                                                                        <  1:1525472/138‑1 (MQ=255)
                                                                                            ggCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGttttt                                                         <  2:1492271/150‑1 (MQ=255)
                                                                                               gcgGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAg                                                      <  2:228860/150‑1 (MQ=255)
                                                                                                     cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTc                                                        >  1:91141/1‑142 (MQ=255)
                                                                                                     cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTc                                                        <  2:91141/142‑1 (MQ=255)
                                                                                                     cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTg                                                >  1:67935/1‑150 (MQ=255)
                                                                                                         cacgcGCCCATGCGTACCTTCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCa                                            <  2:354451/146‑1 (MQ=255)
                                                                                                          aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCAc                                           >  2:326822/1‑150 (MQ=255)
                                                                                                          aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCAc                                           <  1:387974/150‑1 (MQ=255)
                                                                                                              gCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAACAACGTTGGGTTCATGCGCATCGGTCCTTAAGGCTGGACGAACGGACCGGGGTTTTTCAGGTACTGCTCACGGAt                                       >  1:1102679/1‑150 (MQ=255)
                                                                                                                ccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATgg                                     >  1:1232491/1‑150 (MQ=255)
                                                                                                                      gTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCggg                               <  2:162078/150‑1 (MQ=255)
                                                                                                                      gTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCggg                               >  2:1557965/1‑150 (MQ=255)
                                                                                                                        aCCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGcc                             <  2:966064/150‑1 (MQ=255)
                                                                                                                         ccTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCg                            >  1:1580083/1‑150 (MQ=255)
                                                                                                                           tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                          <  2:67935/150‑1 (MQ=255)
                                                                                                                           tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                          <  2:1232491/150‑1 (MQ=255)
                                                                                                                             cgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                        <  2:1580083/150‑1 (MQ=255)
                                                                                                                             cgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                        <  1:1412365/150‑1 (MQ=255)
                                                                                                                               cATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATc                                                                                                                 <  2:1389054/59‑1 (MQ=255)
                                                                                                                               cATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATc                                                                                                                 >  1:1389054/1‑59 (MQ=255)
                                                                                                                                     cATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCCTGCTCCTCGCTCCTCACGCCTGGACCATCAATCCGGTGTAATTCAGGAACATCACTCTAATATGGCGGTCCGACCAGTACGTGa                <  2:430652/150‑1 (MQ=255)
                                                                                                                                      aTTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                           <  2:1151597/90‑1 (MQ=255)
                                                                                                                                      aTTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                           >  1:1151597/1‑90 (MQ=255)
                                                                                                                                          gCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATTCGCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCa           <  1:604303/150‑1 (MQ=255)
                                                                                                                                           cTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGt                                                             <  1:21331/99‑1 (MQ=255)
                                                                                                                                           cTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGt                                                             >  2:21331/1‑99 (MQ=255)
                                                                                                                                            tACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                        >  2:1318446/1‑135 (MQ=255)
                                                                                                                                            tACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGCCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGCTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGc         >  1:1569630/1‑150 (MQ=255)
                                                                                                                                            tACAGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                        <  1:1318446/135‑1 (MQ=255)
                                                                                                                                              cGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGCTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGcca       >  2:165275/1‑150 (MQ=255)
                                                                                                                                               ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTg                                                                              >  2:1473101/1‑78 (MQ=255)
                                                                                                                                               ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTg                                                                              <  1:1473101/78‑1 (MQ=255)
                                                                                                                                                   tgtgCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGc  >  1:1302217/1‑150 (MQ=255)
                                                                                                                                                   tgtgCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGc  >  2:1198434/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGC  >  NC_002947/3829426‑3829722

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