Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,336,010 C→G 5.9% S162T (AGC→ACC)  murA ← UDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,336,0100CG 5.9% 406.2 / 11.8 119S162T (AGC→ACC) murAUDP‑N‑acetylglucosamine 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (52/60);  new base G (5/2);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 119 total aligned reads displayed.

ATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAAC  >  NC_000913/3335862‑3336154
                                                                                                                                                    |                                                                                                                                                
atCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCt                                                                                                                                                 >  1:2095688/1‑150 (MQ=255)
 tCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTg                                                                                                                                                <  1:793730/150‑1 (MQ=255)
   ggTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGAc                                                                                                                                              >  1:1100547/1‑150 (MQ=255)
   ggTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGAc                                                                                                                                              <  2:1718346/150‑1 (MQ=255)
   ggTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAAAAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGAc                                                                                                                                              <  2:2288078/150‑1 (MQ=255)
      gCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCAGCGTTGGCGCAACGGTGAccat                                                                                                                                           <  2:1387278/150‑4 (MQ=255)
      gCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACttt                                                                                                                                           >  2:1187784/1‑150 (MQ=255)
       cccTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGGTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTa                                                                                                                                          <  2:2141167/150‑1 (MQ=255)
        ccTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCAGCGTTGGCGCAACGGTGAccatca                                                                                                                                         <  1:56606/150‑6 (MQ=255)
         cTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATc                                                                                                                                        >  1:2133237/1‑150 (MQ=255)
         cTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATc                                                                                                                                        <  2:498977/150‑1 (MQ=255)
          tGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATcc                                                                                                                                       <  1:24156/150‑1 (MQ=255)
             ccGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATc                                                                                                                                    >  1:1781045/1‑150 (MQ=255)
               gCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCAc                                                                                                                                  <  2:1641875/150‑1 (MQ=255)
               gCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATAGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCAc                                                                                                                                  >  2:710977/1‑150 (MQ=255)
                   aTTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGata                                                                                                                              <  2:845704/150‑1 (MQ=255)
                     tttCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATAtg                                                                                                                            <  1:666839/150‑1 (MQ=255)
                       tCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATAtgtg                                                                                                                          <  2:446304/150‑1 (MQ=255)
                          cACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCAc                                                                                                                       >  2:33918/1‑150 (MQ=255)
                          cACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCAc                                                                                                                       >  2:1338994/1‑150 (MQ=255)
                           aCCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCAcc                                                                                                                      <  2:2250332/150‑1 (MQ=255)
                           aCCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCAcc                                                                                                                      <  2:1232050/150‑1 (MQ=255)
                            cccAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCt                                                                                                                     >  2:396668/1‑150 (MQ=255)
                            cccAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCt                                                                                                                     >  1:1211276/1‑150 (MQ=255)
                              cAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCtta                                                                                                                   >  2:39948/1‑149 (MQ=255)
                               aGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTc                                                                                                                  <  2:660735/150‑1 (MQ=255)
                                     aTCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATATTCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGAAATGTGCACCTTTCAAACgc                                                                                                            >  1:2478868/1‑149 (MQ=255)
                                      tCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGAc                                                                                                           >  2:1464080/1‑150 (MQ=255)
                                      tCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGAc                                                                                                           <  1:2227107/150‑1 (MQ=255)
                                       cAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGAcc                                                                                                          >  1:420170/1‑150 (MQ=255)
                                       cAGGAAGTTCGCGGTATCGACGACTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTCCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGAcc                                                                                                          <  2:2222706/150‑1 (MQ=255)
                                        aGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCa                                                                                                         >  2:374902/1‑150 (MQ=255)
                                         ggAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCAt                                                                                                        >  1:579136/1‑150 (MQ=255)
                                          gAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATc                                                                                                       <  1:2230515/150‑1 (MQ=255)
                                              ttCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACg                                                                                                   <  1:1792924/150‑1 (MQ=255)
                                              ttCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACg                                                                                                   <  2:1729568/150‑1 (MQ=255)
                                                  cgGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAg                                                                                               <  2:2306915/150‑1 (MQ=255)
                                                   ggTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGc                                                                                              <  1:1089836/150‑1 (MQ=255)
                                                   ggTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGc                                                                                              >  2:2466905/1‑150 (MQ=255)
                                                         gacgaTTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAAc                                                                                        >  1:1471139/1‑150 (MQ=255)
                                                         gacgaTTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAAc                                                                                        <  2:729465/150‑1 (MQ=255)
                                                         cacgaTTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTCCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAAc                                                                                        <  2:459847/149‑1 (MQ=255)
                                                           cgaTTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGt                                                                                      >  2:1809222/1‑150 (MQ=255)
                                                              tttCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAc                                                                                   <  2:2173442/150‑1 (MQ=255)
                                                              tttCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAc                                                                                   <  1:1203555/150‑1 (MQ=255)
                                                               ttCCGGTTAACTCGCTGCGTTTTCAATAATCGTGTTTCCTTCAGCCAGGTTTGCATCACACATGATGTTAACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAAcc                                                                                  >  2:1672197/1‑150 (MQ=255)
                                                                tCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACct                                                                                 <  1:1187784/150‑1 (MQ=255)
                                                                 ccGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACctt                                                                                >  2:1920035/1‑150 (MQ=255)
                                                                   ggTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACcttct                                                                              >  2:2218189/1‑150 (MQ=255)
                                                                   ggTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACcttct                                                                              >  1:1561973/1‑150 (MQ=255)
                                                                   ggTTCACGCGCGGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACTCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACcttct                                                                              <  2:2337306/150‑1 (MQ=255)
                                                                    gTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACcttctt                                                                             <  1:1984839/150‑1 (MQ=255)
                                                                    gTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACcttctt                                                                             >  1:1067297/1‑150 (MQ=255)
                                                                         cgcgCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGt                                                                        >  2:1998378/1‑150 (MQ=255)
                                                                           cgcTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGttt                                                                      <  2:1812923/150‑1 (MQ=255)
                                                                            gcTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTg                                                                     <  2:1425451/150‑1 (MQ=255)
                                                                            gcTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTg                                                                     <  2:1434051/150‑1 (MQ=255)
                                                                             cTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGa                                                                    >  2:1511556/1‑150 (MQ=255)
                                                                               gCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATg                                                                  >  1:1563124/1‑150 (MQ=255)
                                                                                cGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATgg                                                                 <  2:934618/150‑1 (MQ=255)
                                                                                 gTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCTCACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGt                                                                >  2:1595135/1‑150 (MQ=255)
                                                                                      cAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCCCCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCgcgc                                                           >  1:2353860/1‑150 (MQ=255)
                                                                                         taatCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTa                                                        <  1:778247/150‑1 (MQ=255)
                                                                                         taatCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTa                                                        >  2:427324/1‑150 (MQ=255)
                                                                                               tggtgCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGtt                                                  >  1:108790/1‑150 (MQ=255)
                                                                                               tggtgCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGtt                                                  >  2:1825253/1‑150 (MQ=255)
                                                                                                ggtgCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTc                                                 <  2:606765/150‑1 (MQ=255)
                                                                                                  tgCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGa                                               >  2:956357/1‑150 (MQ=255)
                                                                                                   gCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAg                                              >  1:1747912/1‑150 (MQ=255)
                                                                                                   gCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAg                                              <  1:215238/150‑1 (MQ=255)
                                                                                                           ccAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGaaa                                      <  1:721734/150‑1 (MQ=255)
                                                                                                             aGGGTTGCAGCACACATGATGGTCACCGTTGCTCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCGAGAAAtg                                    <  1:1805239/150‑1 (MQ=255)
                                                                                                             aGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtg                                    <  2:2258281/150‑1 (MQ=255)
                                                                                                               ggTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtgtg                                  <  1:442434/150‑1 (MQ=255)
                                                                                                               ggTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAAtgtg                                  <  1:1526624/150‑1 (MQ=255)
                                                                                                                gTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATgtgt                                 <  1:299340/150‑1 (MQ=255)
                                                                                                                 ttGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTa                                <  1:430411/150‑1 (MQ=255)
                                                                                                                  tGCAGCACACATGATGGTCACCTTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAATCGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAg                               <  2:435489/150‑1 (MQ=255)
                                                                                                                   gcagcaCACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGa                              <  1:2222821/150‑1 (MQ=255)
                                                                                                                         cacaTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAAcc                        >  1:163603/1‑150 (MQ=255)
                                                                                                                           caTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGCCCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCgg                      >  1:108454/1‑150 (MQ=255)
                                                                                                                           caTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCgg                      >  1:2006513/1‑150 (MQ=255)
                                                                                                                           caTGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCgg                      >  1:1338742/1‑150 (MQ=255)
                                                                                                                             ataaagTCAGCGTTGGCGCAACGGTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGAc                    >  1:603825/6‑150 (MQ=255)
                                                                                                                                tgGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGAcgcg                 <  1:1338994/150‑1 (MQ=255)
                                                                                                                                aagTCAGCGTTGGCGCAACGGTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTTACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGAcgcg                 >  1:822002/3‑150 (MQ=255)
                                                                                                                                aagTCAGCGTTGGCGCAACGGTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCATAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGaaccc                 >  2:310102/3‑146 (MQ=255)
                                                                                                                                  gTCAGCGTTGGCGCAACGGTGACTTTATCCATCACGATATGTGCACCTTTCAAACTACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCTGaacccca               >  2:1175556/1‑144 (MQ=255)
                                                                                                                                  gTCAGCGTTGGCGCAACGGTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCa               >  2:2133800/1‑150 (MQ=255)
                                                                                                                                  gTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCa               <  2:579136/150‑1 (MQ=255)
                                                                                                                                  gTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCa               <  1:1947152/150‑1 (MQ=255)
                                                                                                                                   tCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCAc              >  1:1081134/1‑150 (MQ=255)
                                                                                                                                      ccGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCCCCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGa           >  1:1719878/1‑150 (MQ=255)
                                                                                                                                       cGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGAt          <  2:951015/150‑1 (MQ=255)
                                                                                                                                         ttGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCg        <  1:1998378/150‑1 (MQ=255)
                                                                                                                                           gcgcCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTa      <  1:2466905/150‑1 (MQ=255)
                                                                                                                                            cgcCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTAc     >  2:1055022/1‑150 (MQ=255)
                                                                                                                                            cgcCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTAc     >  1:1492528/1‑150 (MQ=255)
                                                                                                                                               cAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAAc  <  2:1067297/150‑1 (MQ=255)
                                                                                                                                               cAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAAc  <  2:913811/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCGGTGCCCTGACCGCTAATTTTCGCACCCAGCGTAATCAGGAAGTTCGCGGTATCGACGATTTCCGGTTCACGCGCTGCGTTTTCAATAATCGTGGTGCCTTCCGCCAGGGTTGCAGCACACATGATGGTCACCGTTGCGCCAACGCTGACTTTATCCATCACGATATGTGCACCTTTCAAACGACCATCGACGGAAGCTTTAACGTAACCTTCTTCCAGTTTGATGGTCGCGCCTAATTGTTCGAGGCCAGAAATGTGTAGATCAACCGGACGCGCACCGATCGTACAAC  >  NC_000913/3335862‑3336154

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