Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,293,765 C→G 19.5% A123G (GCG→GGG)  lpoA → OM lipoprotein stimulator of MrcA transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,293,7650CG19.5% 115.5 / 5.8 75A123G (GCG→GGG) lpoAOM lipoprotein stimulator of MrcA transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (14/42);  new base G (6/9);  total (21/54)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

AGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAATATTC  >  NC_000913/3293528‑3293998
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
agggttgtttggCCTTTTTTCTTCAGGGAGTGCAGCAAAGCTCGGAGGAATCCAGGATCAACTGGCAATTACTCGCCATTCGGGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGc                                                                                                                                                                                                                                              <  1:197467/227‑1 (MQ=255)
       ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGGGAAAACACTGCTGGCGGGAAGGATTAAACTGGCGCCAGAAGAATTTTCTGGCGCGCAAAACATGCTTGGGGAAAtta                                                                                                                                                                                                                                      >  2:137933/1‑237 (MQ=255)
            ccTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTAACGCAAGAACCGAAACATGCGCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGAGAAAACACGGCGCGAAAGATTTTGCTGTGGCGCCAAAAATGTTGCGGAAAAAAACAccc                                                                                                                                                                                                                                >  1:721472/1‑238 (MQ=255)
              tttttctttCAGCATATTCAGAAAAGCTTGTATGTTACAAGGATCACCTGGAAATTACTCGCCTTTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                                                              <  1:216468/232‑1 (MQ=255)
                   tCTTCAGCAGATGCAGGAAAGCTCTGATGATACCAGGATCAACTGGCAATTCCCCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGAtt                                                                                                                                                                                                                         <  2:843895/239‑1 (MQ=255)
                   tCTTCAGCAGATGCAGCAAAGCTTTGATGATCCCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGAtt                                                                                                                                                                                                                         <  2:8246/239‑1 (MQ=255)
                   tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGc                                                                                                                                                                                                                              <  2:743126/234‑1 (MQ=255)
                                    aaaGCTTTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGc                                                                                                                                                                                                         <  2:207008/237‑1 (MQ=255)
                                      aGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGGGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAgc                                                                                                                                                                                                     <  2:233271/239‑1 (MQ=255)
                                        ctctGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCCCCGAAAGATTTTGATGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACCAAACCACCAAgccc                                                                                                                                                                                                   >  2:449753/1‑237 (MQ=255)
                                           tgatgaTACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCCGGCGGTTGAGTTGTTTAAACCACAACCGCAAGAACTGAACGATGCTCACAGTCGCGCGAAAACACTGCTGGCGGTAGGGGGTAAACTTGCGCCGAAAAAGTTTGTGTGCTCGCAAAACATGCTGAGGCAAAAACAACCAGCCCATTTAAAACaaaaaaacaaaaacccgt                                                                                                                                                                                                 >  1:372113/1‑225 (MQ=255)
                                                atcccAGGACCAACTGGAATTGTCCCCCCATTCGGGCACGGGTAAAAGAGGGTAAAACCGGGCAGGGTGGTGAGGTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCCGCGTCGCGAGAAAACCCTGCGTGCGGTAGAGATTAACCTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCAAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTgg                                                                                                                                                                                           <  1:792803/236‑1 (MQ=255)
                                                       gATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTAACGCAAGAACTGAACGATGCTCAGCGTTGCGGGAAAACAATGCTGGCGGTAAAGATTTAACCTGCGCCAAAAGAATTTTGTGGCGCGCCAAACTTTCTTGGGGAAATCCCACACCCCCGTTTTGAAAAAAAACAACAAACGCGCTAATgggaagaaa                                                                                                                                                                                    >  1:224150/1‑232 (MQ=255)
                                                           aaCTGGCAATTTCTCGCCATTCGTGCACTGGTGAAAGAAGGTAAACCCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAAt                                                                                                                                                                                <  2:712039/239‑1 (MQ=255)
                                                            aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGCTGCTCACCGTCGCGAGAAAACAATGCTGGCGGGAGAGATTAAAATAGGGCACAAAGATTTTGTTCGCGCCGCAAACATCGGTGGGCAAAACACACCCGCCCGTTTTGAAAAAAACCAAAAAAGCAGGTAAATGGCAGGaaaaaaa                                                                                                                                                                                >  1:805455/1‑238 (MQ=255)
                                                            aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGGCCCGAGAAAAAACCTCCTGGGGGAGGAGTTAAAACGGGCACGAAAGAATTTTGTGCGCCGCAAAACATGCTGGCGAAAAAACCACCCGCCGATTTAGAAAAAAAACCACAAACGCGGTACTggggaggaaaaaaa                                                                                                                                                                               >  1:458674/1‑227 (MQ=255)
                                                            aCTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAAATGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAAACTGCGCCGAAAGATTTTGCTGGCGCGCAAAAATTTCTGGGGAAAAACAACCCCTCCCATTTTAAAAAAAAACCAACAAACGCGTAACGGGAGGGCaaaaaa                                                                                                                                                                                >  1:743806/1‑238 (MQ=255)
                                                                 cAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATg                                                                                                                                                                          <  2:106803/238‑1 (MQ=255)
                                                                       cgcGCCATTGGTGCACGGGTAAAAGAAGGTAAAGCCGGGCAGGGGGTTGAGTTGTTTAACAAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATgccagcca                                                                                                                                                                   <  2:872951/237‑1 (MQ=255)
                                                                        tCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAAAGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGGGGAGGATAAAATTGCGAAGGAAAATTTTGTTTGGCCCCAAAACCTGCTGGGGAAAAAAACAACTGCCGATTTATAACAAAAACAAAAAGCCGGGTACTGGCAGGCaaaaaaaagaaacaaccag                                                                                                                                                                   >  1:644844/1‑225 (MQ=255)
                                                                           ttaTTCGGGCACGGGGGAAAGAAGGAAAAAGCGGGGAGGCGGTTGAGTTGTTTAACCACCTCGCGAAAAAACTGAAGGATGCCCGGCGTCCCGAAGAAACCTTGCTGGCGGGAGAGTTTAAACTGGCGGAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGGGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAgggg                                                                                                                                                               <  1:641815/237‑1 (MQ=255)
                                                                           ccttgcGTGCCCTGGTAAAAGAAGGTAAACCCGGGCAGGCGGTTGGGTTGTTTAACCAACTACCGCAAGAACTAAACGATGCCCGGGGTCGCGAAAAAACCCTGCGGGGGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGGGAAAATCCCACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAgggg                                                                                                                                                               <  1:20834/234‑1 (MQ=255)
                                                                               tCGTGCACTGGTGAAAGAAGGTAAAGCCGGGCAGGCGGTTGAGTTTTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTc                                                                                                                                                           <  2:863684/239‑1 (MQ=255)
                                                                               tCGTGCAATGGTGAAAGAAGGTAAAACCGGGCAGGGGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTc                                                                                                                                                           <  2:321545/239‑1 (MQ=255)
                                                                                   gCATTGGTGAAAGAAGGTAAAACCGGGAAGGCGGTTGAGTTTTTTAACCAACTACCGGAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTc                                                                                                                                                       <  2:843954/239‑1 (MQ=255)
                                                                                     ttggggtaaaaagggTAAAACCGGGAGGGGGGTGGGGGTTTTTAAACAATTTCCGCAAGAACGGAACGATGCGCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGGGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  1:71170/226‑1 (MQ=255)
                                                                                     tgtGGTGAAAAAAGGTAAAAGGGGGCAGGGGGTTTGGTTTTTTAACCATCTACCGCAAGACCTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGAGTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  1:253996/237‑1 (MQ=255)
                                                                                     ccgggTGAAAAAAGGTAAACCCGGGCAGGCGTTTGGGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTAAGCGTCGCGAGAAAACACGGCTGGAGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGGGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  1:300176/236‑1 (MQ=255)
                                                                                      cTGGTGAAAGGAGGTAAAACCGGGCCGGGGGGTGGGTTGTTTAACCAACTACCGGAAGAATTGAAGGATGCTCAGCGTAGCGAGAAAACAGTGCTGGGGGTAGAGATTAAACTGGCGCAGAAATATTTTGCGGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  2:737096/238‑1 (MQ=255)
                                                                                         gTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAATTGAACGATGCTCAGCGTCGCGAGAAAACACTGTTGGCGGTGGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAAc                                                                                                                                                                                                              <  2:511995/180‑1 (MQ=255)
                                                                                                   gTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCCGCCGATTTAGAACAAAACCAGCAAACGCGTTACTCGCAGGCAAAAAACGATGCCCCCCAAGGGCGTCCTTCCactttttttctacc                                                                                                                                        >  2:843835/1‑225 (MQ=255)
                                                                                                    ttaaaCGGGGCAGGGGGTTGAGTTTTTTAACCAACTCCCGAAAGAACTGAACGGTGCCACGTGTCCGGAAAAAACGCTGCGGGCGGAAGAGATTAAACGGGCGCAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  1:648640/235‑1 (MQ=255)
                                                                                                    caaaaCCGGGCCGGCGGTTGGGTTTTTTAACCACATACCGCAAAAACGGAACGATGCCAGAGGTCCCGAGAAAACCCGGCTGGCGGTAGAGATTAAACTGGCGCAAAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  1:777828/236‑1 (MQ=255)
                                                                                                    aaaaaCCGGGGAGGGGGTTGGTTTTTTAACCCAACTGCCGGAAGAAATGAACGGTGCCCAGGGCTGGAAGAAAACCCTGCTGGCGGTAGGATTTAAACTGGGGCAGAAAGATTTTGCTGGCGAGAAAAACTTGCTGGGGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  1:748880/236‑1 (MQ=255)
                                                                                                      aaaCCGGGCAGGCGGTTGAGTTGGTTAACCAACTCCCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcgcgc                                                                                                                                     <  2:774167/238‑1 (MQ=255)
                                                                                                          gcggcaggcttttttgtttgttAACCAAACTCCGGAAGAACTGAACGAGTCCGCGCGTCCGAAAAAAACTCTTGTGGGAGTCGATCTTAACCTGGCGCAGAAAGTTTTAGTTGGGCCGTAAACTTTGGTGGAAAAAAACGCACCTGCCGATTTAGACCAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTaa                                                                                                                                <  1:180424/219‑1 (MQ=255)
                                                                                                           gggCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGACCTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTAGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCGGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAAt                                                                                                                               <  2:566816/239‑1 (MQ=255)
                                                                                                            ggCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCCGCGTCGCGCGAAAACACTGCTGGCGGTAAAGATTAAAATTGCGCAGAAAAATTTTGCTGGGCCGCAAAAACTTCTGGCGAAAATTCAACCTGCCCGTTTAGAACAAAAACAACCAACGCGTTACTGGCAAGCAAAAAACGATTCCCAACCCAGGGGGCCCTTCCAtttttttcttcgtccggtaatt                                                                                                                               >  1:506773/1‑219 (MQ=255)
                                                                                                                        tggtttgtaaccccaCTCCCGCAAGAACGGAACGGTGCCCGGCGTGGGGAAGAAACACTGATGGGGGGAGAGTTGAAACGGGCGCAAAAAGATTTTGCTGGCGCGAAAAACTTGCGGGGGAAAAGCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAAcc                                                                                                                   <  1:351362/225‑1 (MQ=255)
                                                                                                                                     aaCTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        <  2:716029/205‑1 (MQ=255)
                                                                                                                                     aaCTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTgcg                                                                                                                                        >  1:716029/1‑205 (MQ=255)
                                                                                                                                       cccccGCAAGAGATGAACATGGCTCGGGGGCGCGAGAAAACCCCGCGGGCGGTAGAGTTAAAACTGCCGCAGAAAGATTTTGCGGGCGCGCAAAACTTGCTGGAGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCgaa                                                                                                    <  1:75633/236‑1 (MQ=255)
                                                                                                                                               aaaaCTGAACGATGCCCAGCGTCGGGAGAAAACCTGGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCGGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGaaaaa                                                                                             <  2:447540/236‑1 (MQ=255)
                                                                                                                                                aaaCTGAACGTGCCTCAGCGTGCGAAAAAACCCCGGCGGGCGGTAAAGATTAACCTGGCGCAAAAAGTTTTTGCTGGCCCGCAAAACTTGCTGGCGAAAAGCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAAca                                                                                           <  1:89338/238‑1 (MQ=255)
                                                                                                                                                 aactaaaggagtgtccgCGTCGCAAAAAAACCGGGCTGGCGGGAAGGATTAAACTGGCGCAGAAAGATTTTGCGGGCGCGCAAAACTGGCTGGTAAAAATCCCACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAAcag                                                                                          <  1:573144/223‑1 (MQ=255)
                                                                                                                                                               ccaggggcgagagaAACCCTTGCGGGCGGTAGAGATTAAACTGGCGCAAAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGAt                                                                              <  2:886240/228‑1 (MQ=255)
                                                                                                                                                               cAGCGTCGCGAGAAAACCCTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGAAAACCTGGCTGGGGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGc                                                                            <  2:89285/239‑1 (MQ=255)
                                                                                                                                                                   ggTCGCGAGAAAACCCGGCGGGCGGTAGAGATTAAACTGGCGGAGAAAGATTTTGCGGGCGCGAAAACCTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCa                                                                                                                                                     <  1:584370/162‑1 (MQ=255)
                                                                                                                                                                      cgcgAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAAATGCTCAGGAACCGCTGCTTGGCGCGCAAGAAAAAAAGAAGAatt                                                                                   >  2:532991/1‑225 (MQ=255)
                                                                                                                                                                           gaAAACCCTGCTGGCGGTAGAGTATAAACTGGCCCAGAAAGATTTGTCTGGCGCCCAAAACTTGCTGGCGAAAACCCCACCTGCCGATTTAGAGCAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAg                                                                  <  2:781085/238‑1 (MQ=255)
                                                                                                                                                                           gaAAACCCTGCGGGCGTGTGATATTAAACTGGCGAAGAAAGATTTTGCTGGCGCGAAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAg                                                                  <  1:137933/238‑1 (MQ=255)
                                                                                                                                                                             aaaCACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGGAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACt                                                                                                                                           >  2:584370/1‑163 (MQ=255)
                                                                                                                                                                              aaCCCTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGgcgc                                                              <  2:763244/239‑1 (MQ=255)
                                                                                                                                                                               acacTGCTGGCGGTAGAGATTAAACTGGCGAAGAAAGATTTTGCCGGCCCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGct                                                             <  2:707558/239‑1 (MQ=255)
                                                                                                                                                                                  ccctgctgGCGGTAGAAATTAAACTGGCGCAGAAAGATTTGTCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGctc                                                            <  2:857912/237‑1 (MQ=255)
                                                                                                                                                                                     tgctgGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCCAAAATCGATGCCACCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCCACGGGACCGCGCCTTGTGCCGGACAAAAAACCACAGAAAATTGTTGCCACCCGGccagggccccc                                                          >  1:309136/1‑228 (MQ=255)
                                                                                                                                                                                       ctgGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACCGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGGAAGAAAAACAGCAAAATATTGATGCCACCTGGCAGGCGCTcccct                                                        >  2:266227/1‑238 (MQ=255)
                                                                                                                                                                                        tgGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCGGCTTGGCGCGAAAGAAAAACAGCCAGATATTGATGCCACCTGGCAGGCGCCCtccccc                                                      >  2:359972/1‑239 (MQ=255)
                                                                                                                                                                                          gggttAGAGATTAAACTGGCGAAAAAAGTTTTTGCTGGCGCGAAAAATTTGCTGGGGAAAACACCCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  1:266227/235‑1 (MQ=255)
                                                                                                                                                                                           ggttAGAGATTAAACTGGCGAAGAAAGATTTTGCTGGCGCGCAAAGCTTGCTGGCGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  2:636199/235‑1 (MQ=255)
                                                                                                                                                                                           cGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAt                                                    <  2:116256/238‑1 (MQ=255)
                                                                                                                                                                                             gTAGAGTTTAAACGGGCGCAGAAAGATTTGGCTGGCCCGCAAACCTTCATGGCGAAAATCACCCCTGCCGATTTAAAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGa                                                  <  2:876308/238‑1 (MQ=255)
                                                                                                                                                                                              aagagATTAAACTGGCGAAGAAAGATTTTGCTGGCGGGAAAAACTTGCGGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGAc                                                 <  2:506773/237‑1 (MQ=255)
                                                                                                                                                                                                    ttAAACGGGCGCAGAAAGATTTGGCTGGCGCGAAAAACTTGCTGGCCAAAATCACCACGTCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAgg                                            <  2:825350/237‑1 (MQ=255)
                                                                                                                                                                                                        aCTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTTATTTACTGCGCGCGTTAATTTCCCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCAtaaacaaagaacaag                                      >  1:281651/1‑225 (MQ=255)
                                                                                                                                                                                                         acgggagcAAAAAGATTTGCTGGAGAGAAAACCTTGCTGGCGAAAATCACCCATGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGGCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGc                                     <  2:9137/233‑1 (MQ=255)
                                                                                                                                                                                                              gccggaaaacttttGCTGGAGCACAAACCGTGCTGGCAAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCTTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGaa                                  <  2:186299/228‑1 (MQ=255)
                                                                                                                                                                                                              cgcAAAAAGATTTTGCTGGCGCGAAAAACTTGCGGGCAAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGaa                                  <  2:644844/238‑1 (MQ=255)
                                                                                                                                                                                                                  gAAAGTTTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACg                              <  2:594326/238‑1 (MQ=255)
                                                                                                                                                                                                                           tGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGAATCAGGAACAGGCGACTAAGCTGGGGTTCa                    >  2:831162/1‑239 (MQ=255)
                                                                                                                                                                                                                           tGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGAAGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCCGGGGATCa                    >  2:14414/1‑239 (MQ=255)
                                                                                                                                                                                                                              gggCGCCGAAAACTTGCTGGCAAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACg                 <  2:721472/238‑1 (MQ=255)
                                                                                                                                                                                                                                 cgcgcAAAACTTGCTGGCGAAAATCCCCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGcc               <  2:743806/238‑1 (MQ=255)
                                                                                                                                                                                                                                 cgcgaAAATCTTGCGGGGAAAAACCCCCCCTGCCGATTTAAAAAAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAAGCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGcc               <  1:843835/238‑1 (MQ=255)
                                                                                                                                                                                                                                           cTTGCTGGCGAAAATCCCCCCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAata     <  2:224150/239‑1 (MQ=255)
                                                                                                                                                                                                                                              gCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTTATGCCACCTGGCAGGCGCTCCCCCCCATGACTCAGGAAAAGGGGGAAAACCGGGGGATAAAAGCCGACCaaaaaaatt  >  1:644498/1‑234 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            
AGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAACACTGCTGGCGGTAGAGATTAAACTGGCGCAGAAAGATTTTGCTGGCGCGCAAAACTTGCTGGCGAAAATCACACCTGCCGATTTAGAACAAAACCAGCAAGCGCGTTACTGGCAGGCAAAAATCGATGCCAGCCAGGGGCGTCCTTCCATTGATTTACTGCGCGCGTTAATTGCTCAGGAACCGCTGCTTGGCGCGAAAGAAAAACAGCAGAATATTGATGCCACCTGGCAGGCGCTCTCCTCCATGACTCAGGAACAGGCGAATACGCTGGTGATCAACGCCGACGAAAATATTC  >  NC_000913/3293528‑3293998

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