Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 910,972 G→A G26G (GGC→GGT poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913910,9720GA95.0% 194.4 / ‑2.2 62G26G (GGC→GGTpoxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (22/35);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGA  >  NC_000913/910827‑911119
                                                                                                                                                 |                                                                                                                                                   
cGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCtgt                                                                                                                                                <  1:648766/151‑1 (MQ=255)
cGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCtgt                                                                                                                                                <  2:897913/151‑1 (MQ=255)
cGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCtgt                                                                                                                                                <  2:144856/151‑1 (MQ=255)
             gcAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATg                                                                                                                                   <  2:267866/151‑1 (MQ=255)
              cAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGc                                                                                                                                  <  2:805999/151‑1 (MQ=255)
               aGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGc                                                                                                                                  <  1:1154153/150‑1 (MQ=255)
               aGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCg                                                                                                                                 <  2:1039780/151‑1 (MQ=255)
                    gCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCa                                                                                                                            <  2:10651/151‑1 (MQ=255)
                     ccAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGCATCCACTCCATGGGGCCCCCGCGATTAAGACTGTCACCAAAACCGTTCAGAGAGTCACCTGTGACTCCCCCGAAGCGGTTCAc                                                                                                                           >  1:1007559/1‑151 (MQ=255)
                     ccAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCAc                                                                                                                           <  2:1031040/151‑1 (MQ=255)
                          tctcCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCt                                                                                                                       <  2:320276/150‑1 (MQ=255)
                                      tgtgCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGCATCCCCCCCATGGGGGCCCTGCGGTTAAGACTGTCACTAAGACCGCTCAGAGAGTCCCCTGTGACTCCCCAGATGCGTTTTACCCCCGCCGATTCGGgtg                                                                                                          >  1:412463/1‑151 (MQ=255)
                                      tgtgCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCTCCATTGTGCCCATTCGAGTAAAACTGTCACTAAGACCGTTTAGAGAGTAACCTCTGGCTCCCCCGATGCGGTTCACCCCCGcccgatcgtg                                                                                                            >  2:138895/1‑141 (MQ=255)
                                       gtgGTTCGGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAgtg                                                                                                          <  1:822478/150‑1 (MQ=255)
                                         gCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGttt                                                                                                       <  2:430460/151‑1 (MQ=255)
                                            tCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGtttt                                                                                                      <  1:885458/149‑1 (MQ=255)
                                              aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGAATCCACTCGATGGGGCCCCTGCGATTAAGACTGTCACCAAAAACGTTCAGAGAGGCACCCGTGACTCCCCCGATGCGGTTTACCCCCGCCCATTCGAGTGtttt                                                                                                      >  1:785112/1‑147 (MQ=255)
                                                cgcgAGCGGCAAAGGCCGCCTCTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGGCTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGGGTTTCGCCCCTGCCGATTCGAGTGTTTTGGCGat                                                                                                <  1:1179917/151‑1 (MQ=255)
                                                     gCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAg                                                                                           <  2:260180/151‑1 (MQ=255)
                                                      cGGCAAAGGCCGCCACTTCTTCGTGGCGGGAGGGCATCCACTCGATGGTGGCCCTGCGATTAAAACCGTCCCTAAGAACGGTCAGAAAGTGACCTGTGACTCCCCACATGCGGTTCAACCCCCCCGATTTGAGGGTTTTTGCGGTATaa                                                                                            >  2:818621/1‑149 (MQ=255)
                                                         cAAAGGCCGCCCCTTCTTCGTGGCGGGTGGACATCCTCTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTg                                                                                        <  2:325755/150‑1 (MQ=255)
                                                         cAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACTCCCTGGGGCCCATGCGATTTAAACTGTCACTAAAACCGTTTAGAAAGGCACCTGTGACTCCCCAGATGCGGTTTACCCCCGCCGGTTCGAGTGTTTTGGCGATATAAACt                                                                                         >  2:520685/1‑149 (MQ=255)
                                                           aaGGCAGCCACTTCTGCGTGGCGGATGGCGATCCACTCGATGATGCACGAGCGATTAAGACTATCACTAAGCCCGTTAAGAGAGTCGCCTGTGACACCCCAGATGCGTTTCTCCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCa                                                                                      <  1:1043644/150‑1 (MQ=255)
                                                            aGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCCATGGTGGCCCCGCCATTTAGACTGTCACTCAGACCGTTCAGAGAGTCACCTGTGAACCCCCCGATGCGTTTCACCCCCGCCCGTTTGAGGGTTTTTGCGATATAAGCTGCa                                                                                      >  1:1050110/1‑149 (MQ=255)
                                                                 gccACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCgttt                                                                               <  1:224941/151‑1 (MQ=255)
                                                                  ccACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCgtttg                                                                              <  2:904376/151‑1 (MQ=255)
                                                                  ccACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCgtttg                                                                              <  1:455036/151‑1 (MQ=255)
                                                                             tGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGtt                                                                    >  2:175334/1‑150 (MQ=255)
                                                                              ggCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTtct                                                                  <  1:41787/151‑1 (MQ=255)
                                                                                cGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTtctc                                                                 >  1:806546/1‑150 (MQ=255)
                                                                                 gggTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCa                                                               <  2:38300/151‑1 (MQ=255)
                                                                                 gggTGGACATCCACTCGATGGCGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGCTTCCCCAGATGCGTTTCACCCCTGCCGATTCGTGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTtctc                                                                 >  1:691079/1‑149 (MQ=255)
                                                                                  ggTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCAt                                                              <  2:331232/151‑1 (MQ=255)
                                                                                    tGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCAt                                                              >  2:441791/1‑149 (MQ=255)
                                                                                        cATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTg                                                        >  2:521038/1‑151 (MQ=255)
                                                                                            cACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAAtg                                                    >  2:773921/1‑151 (MQ=255)
                                                                                             aCTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAAtgt                                                   <  2:189719/151‑1 (MQ=255)
                                                                                                cGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAAtgtg                                                  <  2:1187326/149‑1 (MQ=255)
                                                                                                cGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGAt                                                >  2:770928/1‑151 (MQ=255)
                                                                                                  aTGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATa                                               >  1:511298/1‑150 (MQ=255)
                                                                                                     gTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACgg                                           <  2:270390/151‑1 (MQ=255)
                                                                                                       gCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTa                                         <  2:278992/151‑1 (MQ=255)
                                                                                                        cccATGCGATTAAGACTGTCACTAAGACCGCTCAGAGAGTCCCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTaa                                        <  1:31838/151‑1 (MQ=255)
                                                                                                          cATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAAc                                       >  2:900164/1‑150 (MQ=255)
                                                                                                           aTGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACa                                      >  1:680384/1‑150 (MQ=255)
                                                                                                           aTGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTCGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACaa                                     <  1:268308/151‑1 (MQ=255)
                                                                                                              cGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGtt                                  >  1:549575/1‑151 (MQ=255)
                                                                                                                 ttAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGttt                                 >  1:489779/1‑149 (MQ=255)
                                                                                                                  tAAGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGt                              >  1:1139778/1‑151 (MQ=255)
                                                                                                                   aaGACTGTCCCTAAGACCGTTCAGAGAGTCACCTGTGACTCACCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGt                              <  1:594565/150‑1 (MQ=255)
                                                                                                                   aaGACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGtt                             <  1:258918/151‑1 (MQ=255)
                                                                                                                     gACTGTCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTc                            >  2:384694/1‑150 (MQ=255)
                                                                                                                          tCACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTg                       >  2:724885/1‑150 (MQ=255)
                                                                                                                           cACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGa                      <  2:349901/150‑1 (MQ=255)
                                                                                                                           cACTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGAc                     >  1:318560/1‑151 (MQ=255)
                                                                                                                            aCTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACg                    <  1:302924/151‑1 (MQ=255)
                                                                                                                            aCTAAGACCGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACg                    <  1:237331/151‑1 (MQ=255)
                                                                                                                                    cGTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGa            <  1:1009393/151‑1 (MQ=255)
                                                                                                                                     gTTCAGAGAGTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGaa           <  2:771698/151‑1 (MQ=255)
                                                                                                                                          gagagTCACCTGTGTCTCCCCAGATGCGTTTCACCCCTGACGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTAATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGAt      <  1:1234321/151‑1 (MQ=255)
                                                                                                                                            gagTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATgg    <  1:631530/151‑1 (MQ=255)
                                                                                                                                              gTCACCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGga  <  1:80489/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGCGATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGA  >  NC_000913/910827‑911119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: