Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 809,038 C→T A74T (GCC→ACC)  bioA ← 7,8‑diaminopelargonic acid synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913809,0380CT96.0% 260.1 / ‑1.1 77A74T (GCC→ACC) bioA7,8‑diaminopelargonic acid synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (42/31);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCG  >  NC_000913/808892‑809184
                                                                                                                                                  |                                                                                                                                                  
cATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCa                                                                                                                                                <  2:737518/151‑1 (MQ=255)
  ttttCATCGCCACTTCCACCGCTACAGAACTGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCAACACAGCTCAATGGCTGGCGCATGGGTGATACCGCGAAACATCACATGCGACATGATATCAAt                                                                                                                                              >  1:170732/1‑151 (MQ=255)
      cATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAAttt                                                                                                                                            <  1:581069/149‑1 (MQ=255)
      cATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGc                                                                                                                                          >  1:977906/1‑151 (MQ=255)
           ccACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACtt                                                                                                                                     <  1:442850/151‑1 (MQ=255)
           ccACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACtt                                                                                                                                     >  2:1223851/1‑151 (MQ=255)
            cACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTc                                                                                                                                    >  2:634293/1‑151 (MQ=255)
            cACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTc                                                                                                                                    <  1:673397/151‑1 (MQ=255)
            cACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTc                                                                                                                                    >  2:739658/1‑151 (MQ=255)
                tCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAAGGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCAt                                                                                                                                  <  2:897237/149‑1 (MQ=255)
                  cACCGCTGCGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCTGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATcgcc                                                                                                                              <  1:517667/151‑1 (MQ=255)
                   aCCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTCCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGAATCAATTTGCGACTTCATcgcc                                                                                                                              >  1:650201/1‑150 (MQ=255)
                         aCGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTaa                                                                                                                       <  1:929439/151‑1 (MQ=255)
                                cGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCgg                                                                                                                >  2:244012/1‑151 (MQ=255)
                                  gAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCggg                                                                                                               >  1:1338030/1‑150 (MQ=255)
                                          gagaAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCGCCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAg                                                                                                      <  1:467865/151‑1 (MQ=255)
                                          gagaAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGt                                                                                                        >  2:367832/1‑149 (MQ=255)
                                             aaaaaCGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCg                                                                                                   >  1:193876/1‑151 (MQ=255)
                                              aaaaCGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGt                                                                                                  <  2:108235/151‑1 (MQ=255)
                                              aaaaCGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGt                                                                                                  >  2:327899/1‑151 (MQ=255)
                                                  cGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGAt                                                                                              <  1:1184767/151‑1 (MQ=255)
                                                  cGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGAt                                                                                              <  2:968351/151‑1 (MQ=255)
                                                      cTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATcgcc                                                                                          >  2:403264/1‑151 (MQ=255)
                                                          aGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGccc                                                                                      >  2:831821/1‑151 (MQ=255)
                                                          aGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGccc                                                                                      >  2:881937/1‑151 (MQ=255)
                                                             ggTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATCTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCccacc                                                                                   >  1:578196/1‑151 (MQ=255)
                                                                 gcggcgTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAgg                                                                                <  2:537499/150‑1 (MQ=255)
                                                                  cggcgTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGgac                                                                              <  2:145245/151‑1 (MQ=255)
                                                                     cgTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGgacgac                                                                           >  1:689741/1‑151 (MQ=255)
                                                                      gTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGgacgac                                                                           >  2:590263/1‑150 (MQ=255)
                                                                            gCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGt                                                                    >  2:271174/1‑151 (MQ=255)
                                                                             ccaccaGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGGGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTc                                                                   >  1:1171974/1‑151 (MQ=255)
                                                                              caccaGTTTGCGACACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTc                                                                   >  2:119959/1‑150 (MQ=255)
                                                                              caccaATTTTCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATAc                                                                       <  2:991662/146‑1 (MQ=255)
                                                                               accaGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCaa                                                                 <  2:951054/151‑1 (MQ=255)
                                                                                  aGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAAcc                                                               >  1:268944/1‑150 (MQ=255)
                                                                                      tGCGGCACAGCTCAATGGCTGGCGCATGGGGGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGCGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGc                                                           >  1:356894/1‑150 (MQ=255)
                                                                                       gCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGt                                                         <  2:648138/151‑1 (MQ=255)
                                                                                        cGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTc                                                        <  1:1175637/151‑1 (MQ=255)
                                                                                         ggCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCt                                                       >  1:458098/1‑151 (MQ=255)
                                                                                              aGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGt                                                  >  1:109300/1‑151 (MQ=255)
                                                                                                cTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTc                                                 >  2:362212/1‑150 (MQ=255)
                                                                                                 tCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTc                                                 >  1:1310543/1‑149 (MQ=255)
                                                                                                   aaTGGCTGGCGCATGGGTGATACCGCCAACCATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCGCCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGa                                              <  1:115418/150‑1 (MQ=255)
                                                                                                   aaTGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGa                                              <  2:1223015/150‑1 (MQ=255)
                                                                                                   aaTGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGAAATGGTATCAATTTGCGACTTCATGGCCGCATTAACCTACGGGTGATTGTAGCCGAGGATCCCCGCCCACCAGGCCGACATACCGTCAACCAGGCGTCTGCCGTc                                                 >  1:788119/1‑147 (MQ=255)
                                                                                                      ggctggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACaa                                           >  1:558277/1‑150 (MQ=255)
                                                                                                      ggctggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACaa                                           >  2:907358/1‑150 (MQ=255)
                                                                                                       gctggcGCATGCGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATGCCGTCAACCAGGCGTCTGCCGTCAGACaaaa                                         <  1:215096/151‑1 (MQ=255)
                                                                                                        ctggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAAt                                        >  2:552562/1‑151 (MQ=255)
                                                                                                             gcATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCt                                   >  1:1114590/1‑151 (MQ=255)
                                                                                                             gcATGGGTGATACCGCCAAACATCACATGAGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGc                                    >  1:521064/1‑150 (MQ=255)
                                                                                                              cATGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCGCCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTc                                  <  2:1081132/151‑1 (MQ=255)
                                                                                                               aTGGGTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCg                                 >  1:780725/1‑151 (MQ=255)
                                                                                                               aTGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCg                                 >  1:822341/1‑151 (MQ=255)
                                                                                                                 gggTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGc                                >  1:965002/1‑150 (MQ=255)
                                                                                                                 gggTGATACCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGc                                <  2:193885/150‑1 (MQ=255)
                                                                                                                      aTACCGCCAAACATCACATGCGACATGGTATCAATTTTCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAAcc                            <  1:722463/149‑1 (MQ=255)
                                                                                                                        aCCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTc                         <  2:456319/150‑1 (MQ=255)
                                                                                                                        aCCGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTc                         <  2:702085/150‑1 (MQ=255)
                                                                                                                          cGCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGgc                      >  1:1003713/1‑151 (MQ=255)
                                                                                                                           gCCAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGgcg                     <  1:979520/151‑1 (MQ=255)
                                                                                                                            ccAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGgcgc                    >  2:908687/1‑151 (MQ=255)
                                                                                                                             cAAACATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGgcgc                    >  2:1094505/1‑150 (MQ=255)
                                                                                                                              aaaCATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCCCCAGGACGACATCCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTc                  <  1:1308545/151‑1 (MQ=255)
                                                                                                                               aaCATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTca                 <  1:725394/151‑1 (MQ=255)
                                                                                                                                 cATCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTcacc               >  2:70579/1‑151 (MQ=255)
                                                                                                                                  aTCACATGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTcacc               >  2:570115/1‑150 (MQ=255)
                                                                                                                                    cacaTGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCaccacc            <  2:512226/151‑1 (MQ=255)
                                                                                                                                    cacaTGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCaccacc            <  2:1335630/151‑1 (MQ=255)
                                                                                                                                     acaTGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCg           >  2:830307/1‑151 (MQ=255)
                                                                                                                                       aTGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCgg          >  1:185057/1‑150 (MQ=255)
                                                                                                                                       aTGCGACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGa         <  2:916443/151‑1 (MQ=255)
                                                                                                                                           gACATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATaa      >  1:1275092/1‑150 (MQ=255)
                                                                                                                                            aCATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAAc    >  2:1190616/1‑151 (MQ=255)
                                                                                                                                             cATGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAAcc   >  1:436988/1‑151 (MQ=255)
                                                                                                                                               tGGTATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCg  <  1:1347387/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CATTTTCATCGCCACTTCCACCGCTACGGAACCGGAGTCCGCGAGAAAAACGCACTCCAGCGGTTGCGGCGTCATCGCCACCAGTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCG  >  NC_000913/808892‑809184

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