Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,497,025 C→A 9.4% intergenic (+33/+32) ypdK → / ← alaC inner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,497,0250CA 9.4% 265.3 / 16.7 98intergenic (+33/+32)ypdK/alaCinner membrane protein/glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base C (45/42);  new base A (3/6);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

TGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATC  >  NC_000913/2496877‑2497168
                                                                                                                                                    |                                                                                                                                               
tGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCa                                                                                                                                                <  1:998894/150‑1 (MQ=255)
 gAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGATATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAg                                                                                                                                               >  2:1049568/1‑150 (MQ=255)
   aGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                             <  1:1718351/150‑1 (MQ=255)
   aGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                             <  2:4197/150‑1 (MQ=255)
   aGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCAGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGc                                                                                                                                             <  2:112138/150‑1 (MQ=255)
    gTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCt                                                                                                                                            >  2:206869/1‑150 (MQ=255)
         tCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGt                                                                                                                                       <  2:1211942/150‑1 (MQ=255)
         tCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGt                                                                                                                                       <  2:716190/150‑1 (MQ=255)
          cAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtt                                                                                                                                      <  1:763014/150‑1 (MQ=255)
          cAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTACAATGAAAAACAGGAGCCATCAGGCTCCTGtt                                                                                                                                      >  1:862125/1‑150 (MQ=255)
           aGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttt                                                                                                                                     <  1:637100/150‑1 (MQ=255)
            ggCGTAACAGTAATTACGCGGAGAGATGTATAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttt                                                                                                                                    <  1:150178/150‑2 (MQ=255)
            ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGTAATATTTCTTTATTGGCATTTCTTTTATGGTCCTCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACCGGAGCCTGATGGCTCCTGtttt                                                                                                                                    <  1:766383/150‑2 (MQ=255)
            ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGtttt                                                                                                                                    <  1:945536/150‑2 (MQ=255)
            ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTg                                                                                                                                    <  1:816993/150‑1 (MQ=255)
            ggCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTg                                                                                                                                    <  2:912921/150‑1 (MQ=255)
             gCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCTGATGGCTCCTGttttt                                                                                                                                   <  2:1465658/150‑3 (MQ=255)
               gTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGttagaat                                                                                                                                 >  2:214837/1‑145 (MQ=255)
               gTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCAt                                                                                                                                 >  2:337553/1‑150 (MQ=255)
                  aCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTc                                                                                                                              >  2:466772/1‑150 (MQ=255)
                  aCAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTc                                                                                                                              >  1:1019072/1‑150 (MQ=255)
                   cAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCACTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGTATTTCt                                                                                                                             >  1:1503227/1‑150 (MQ=255)
                   cAGTAATTACGCGGAGAGATCTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGCTCATCGATTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGTATTTCt                                                                                                                             >  1:1503223/1‑150 (MQ=255)
                        aTTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTgacga                                                                                                                        >  1:882274/1‑150 (MQ=255)
                          tACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGacgacg                                                                                                                      >  1:841143/1‑150 (MQ=255)
                           aCGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACgg                                                                                                                     >  1:210671/1‑150 (MQ=255)
                           aCGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACgg                                                                                                                     <  2:14765/150‑1 (MQ=255)
                              cgGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGATTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAAGTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGgtgt                                                                                                                  >  1:1043093/1‑150 (MQ=255)
                               ggAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGtt                                                                                                                 <  1:122757/150‑1 (MQ=255)
                                    gaTGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAAATGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTc                                                                                                            >  1:1014835/1‑150 (MQ=255)
                                       gTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcgc                                                                                                         <  1:466772/150‑1 (MQ=255)
                                        tAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCcgcg                                                                                                        >  2:726462/1‑150 (MQ=255)
                                          aaGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGtt                                                                                                      <  1:1641560/150‑1 (MQ=255)
                                            gTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGtttt                                                                                                    >  1:1438972/1‑150 (MQ=255)
                                             tGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTc                                                                                                   >  2:1102378/1‑150 (MQ=255)
                                               aaaTATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGt                                                                                                 <  1:243153/150‑1 (MQ=255)
                                                 atatTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGAAATGAAAAACAGGAGCCATCAGGATCCTGTTTGCATTTCTGACGACGGTGTTTATTACGCGTTTTCGTGa                                                                                               >  1:18829/1‑150 (MQ=255)
                                                 atatTTCTTTATGGGCATTTCTTTTATGGTAATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGa                                                                                               <  2:938925/150‑1 (MQ=255)
                                                    tttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAata                                                                                            <  2:1099128/150‑1 (MQ=255)
                                                    tttCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAata                                                                                            <  1:496459/150‑1 (MQ=255)
                                                     ttCTTTATGGGCATTTCTTTTATGGTCATCGATTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTCTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTCAatat                                                                                           >  2:942781/1‑150 (MQ=255)
                                                        tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                >  2:277806/1‑142 (MQ=255)
                                                        tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTg                                                                                                <  1:277806/142‑1 (MQ=255)
                                                        tttATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGtt                                                                                        <  2:959268/150‑1 (MQ=255)
                                                          tATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtg                                                                                      <  2:882274/150‑1 (MQ=255)
                                                           aTGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgc                                                                                     <  2:1441019/150‑1 (MQ=255)
                                                            tGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgct                                                                                    >  1:1011840/1‑150 (MQ=255)
                                                            tGGGCATTTCTTTTATGGACATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCAATTCTGTCGACGGTGTTTATTCCGCGTTTTCGTCAATATGTTtgct                                                                                    >  1:259552/1‑150 (MQ=255)
                                                              ggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgctgc                                                                                  <  2:716546/150‑1 (MQ=255)
                                                              ggCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTtgctgc                                                                                  >  2:1014380/1‑150 (MQ=255)
                                                               gCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTgctgct                                                                                 >  1:472507/1‑150 (MQ=255)
                                                                cATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGctgctg                                                                                >  1:1666668/1‑150 (MQ=255)
                                                                 aTTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTgg                                                                               >  1:409543/1‑150 (MQ=255)
                                                                 aTTTCATTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTgg                                                                               <  1:1049568/150‑1 (MQ=255)
                                                                    tCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACCGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCgg                                                                            <  1:820477/150‑1 (MQ=255)
                                                                     cTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCggg                                                                           <  1:798823/150‑1 (MQ=255)
                                                                        ttATGGTCATCGTTTTGGCCGGTACTTTTGCCCTGATTATCTAAAAGCAAACTGTGCAATGAATAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGgtag                                                                        <  2:360448/150‑1 (MQ=255)
                                                                           tGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTaa                                                                     >  1:790005/1‑150 (MQ=255)
                                                                            ggTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTaaa                                                                    >  2:513451/1‑150 (MQ=255)
                                                                            ggTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTaaa                                                                    >  1:992265/1‑150 (MQ=255)
                                                                             gTCATCGTTGGGGCCGGTACTTTTGCCCTGATGATCAAAAAGCAAACTGAGCACTGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGGTTGCTGCTGGCGGGTAGTAAAc                                                                   <  2:280890/150‑1 (MQ=255)
                                                                                 tCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACATGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTc                                                               <  2:1440932/150‑1 (MQ=255)
                                                                                 tCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTTAATATGTTTGCTTCTGGAGGGTAGTAAACCggc                                                               <  2:263453/150‑3 (MQ=255)
                                                                                  cGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATTAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCg                                                              >  2:1348876/1‑150 (MQ=255)
                                                                                   gTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCgg                                                             <  2:1019072/150‑1 (MQ=255)
                                                                                   gTTTGGGACGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAATCAGGAGCCATCAGGCTACTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCgg                                                             <  2:97181/150‑1 (MQ=255)
                                                                                                 ttttGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCt                                               >  2:891465/1‑150 (MQ=255)
                                                                                                    tGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTa                                            >  1:135444/1‑150 (MQ=255)
                                                                                                     gCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTaa                                           <  1:1102378/150‑1 (MQ=255)
                                                                                                      cccTGATGATCTAAAAGCAAACTGTGCTATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                          >  2:874447/1‑150 (MQ=255)
                                                                                                      cccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTTTTTGTATTTCTGACGAGGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                          <  2:1503223/150‑1 (MQ=255)
                                                                                                      cccTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                          >  2:1674590/1‑150 (MQ=255)
                                                                                                      cccGGATGATCTAAAAGCAAACGGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGTATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAAt                                          <  2:1503227/150‑1 (MQ=255)
                                                                                                        cTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGTCACATCGCTTTAATcc                                        >  2:167143/1‑150 (MQ=255)
                                                                                                        cTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATcc                                        <  1:891465/150‑1 (MQ=255)
                                                                                                                cTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATTAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                <  1:1348876/150‑1 (MQ=255)
                                                                                                                cTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATg                                >  1:1530664/1‑150 (MQ=255)
                                                                                                                  aaaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGc                              >  2:736339/1‑150 (MQ=255)
                                                                                                                  aaaaGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGc                              <  1:1276872/150‑1 (MQ=255)
                                                                                                                     aGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTg                           <  2:333255/150‑1 (MQ=255)
                                                                                                                        aaaCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                        <  2:1666668/150‑1 (MQ=255)
                                                                                                                        aaaCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                        <  1:337553/150‑1 (MQ=255)
                                                                                                                        aaaCTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                        <  1:652551/150‑1 (MQ=255)
                                                                                                                            tgtgCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                        <  2:979981/146‑1 (MQ=255)
                                                                                                                            tgtgCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACg                        >  1:979981/1‑146 (MQ=255)
                                                                                                                                cAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGt                >  2:1280741/1‑150 (MQ=255)
                                                                                                                                   atgcaaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTcgc             >  2:48271/5‑150 (MQ=255)
                                                                                                                                    gAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTcgcg            >  1:646871/1‑150 (MQ=255)
                                                                                                                                    gAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTcgcg            >  1:1343671/1‑150 (MQ=255)
                                                                                                                                     aaaaaCAGGAGCCATCAGGGTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCgg           >  1:1065221/1‑150 (MQ=255)
                                                                                                                                      caaaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGATGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGt          >  1:976032/2‑150 (MQ=255)
                                                                                                                                      aaaaCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCTCTCTGATGGCCTGACGAATACGGTCGCGGt          >  1:1072380/1‑150 (MQ=255)
                                                                                                                                        aaCAGGAGCCTGATGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGttt        >  1:558905/1‑150 (MQ=255)
                                                                                                                                         aCAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGtttt       <  2:1124496/150‑1 (MQ=255)
                                                                                                                                          cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTATACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTc      <  1:214837/150‑1 (MQ=255)
                                                                                                                                          cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTc      >  1:707121/1‑150 (MQ=255)
                                                                                                                                          cAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTGTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTc      >  1:1011182/1‑150 (MQ=255)
                                                                                                                                              aGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATc  >  2:770992/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGAAGTCATTCAGGCGTAACAGTAATTACGCGGAGAGATGTAAAGTGAAATATTTCTTTATGGGCATTTCTTTTATGGTCATCGTTTGGGCCGGTACTTTTGCCCTGATGATCTAAAAGCAAACTGTGCAATGAAAAACAGGAGCCATCAGGCTCCTGTTTGCATTTCTGACGACGGTGTTTATTCCGCGTTTTCGTGAATATGTTTGCTGCTGGCGGGTAGTAAACCGTCGGCACGGAACATCGCTTTAATCCCTCTGATGGCCTGACGAATACGGTCGCGGTTTTCAATC  >  NC_000913/2496877‑2497168

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