Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 665,104 A→C 10.2% S29R (AGT→AGG rlpA ← septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913665,1040AC10.2% 212.3 / 3.1 80S29R (AGT→AGGrlpAseptal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (41/29);  new base C (3/7);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCG  >  NC_000913/664870‑665336
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
ctcgccttgttTCGGATGGAATGCTCCCCCGGGGGCCGTGAGGTTCGCGCCGGTTTGGGCCTTATGGGTTCCCCCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCGGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCATCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                                                                                                                                                                                      <  2:773742/232‑1 (MQ=255)
    aGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCATTGGTTCGAAACGCGGTTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGAAGGCTGCGGTACACTTACCg                                                                                                                                                                                                                                   >  1:138623/1‑238 (MQ=255)
     gCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGCGACGGATCCTGCACGCTTTTGTAGCTTTTACCGTCGCCCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGCCCCCCCCGCTAATTTCACCTATAGGGCGCGTACATACCCACGGCTCCGGCCCACTTACCgc                                                                                                                                                                                                                                  >  2:157616/1‑237 (MQ=255)
         tgtCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGCTTTTGTAGCTTTTACCGTCGCGCTGGTAATCCCGATTTGCCGTCACGTTCCTGGGTTCGAAACCGCGGTCCGCCCCGCTAATTTCAACTATAGGCCCGTTACATACCCCAGGCTCCGGTCACCTTACCGCCTGtt                                                                                                                                                                                                                             >  2:666444/1‑239 (MQ=255)
         tgtCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCCCAGGCTGCGGTACACTTACCGCCTGtt                                                                                                                                                                                                                             >  1:174811/1‑239 (MQ=255)
         tgtCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGGCCGCCCCGCTAATTTCAACTATAGGGCCGTTAAAAACCGCAGGCTGCGGTCCACTTACCGGCCGtt                                                                                                                                                                                                                             >  2:54107/1‑239 (MQ=255)
                tCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCCCAGGCTCCGGGACCATTACCGCCCGTTGCTGAcc                                                                                                                                                                                                                      >  2:257146/1‑239 (MQ=255)
                   aaTGTTTCCCCAGGGGCCGTCAGGTTACCGCCGGGTTCGGAATCAATGTTTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatc                                                                                                                                                                                                                   <  1:257146/239‑1 (MQ=255)
                   aaTGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACCCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCCGGCTGCGGTACACTTACCTTCTGTTGCTGACcatc                                                                                                                                                                                                                   >  1:955059/1‑239 (MQ=255)
                    ttGCTCCGCCGGAGGCCGTCAGGTTACGGCCGGGTTCGGCATCATAAGTTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatc                                                                                                                                                                                                                   <  1:371668/237‑1 (MQ=255)
                    atgccgcGCCGGGGGCCGTGGGGTTACGGCCGGGTTGGGCATAATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatc                                                                                                                                                                                                                   <  1:157616/232‑1 (MQ=255)
                     tgtctccccAGGGGGCGTCGGGTTGCCGCCGGGTTCGGCTCAATGGATTGCCCGCGGGCCCGCCTGGGTAAACCGAGCCGGATCCTGCCCGTTTGTGTAGCTTTTCCCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                                                                                                                                                                                 <  2:188926/232‑1 (MQ=255)
                        ttCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCCCAGGCTGCGGTCCACTTCCCGCCTGTTGCTGACCatcatc                                                                                                                                                                                                                >  1:968806/1‑237 (MQ=255)
                         tCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCt                                                                                                                                                                                                             >  1:662410/1‑239 (MQ=255)
                                    cGTCAGGTTACTGCCGAGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAAACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGAAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGCCTGTTGCTGACCATAATCCCTTGTACATTccc                                                                                                                                                                                                  >  1:673887/1‑238 (MQ=255)
                                      tCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTCTAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGcccc                                                                                                                                                                                                 >  1:548844/1‑236 (MQ=255)
                                                     gtgcGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAAACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc                                                                                                                                                                                  <  1:54107/235‑1 (MQ=255)
                                                                     tGCCCCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGcc                                                                                                                                                                 <  1:568741/239‑1 (MQ=255)
                                                                     gtcctccATTCCCGCCTGGTTAAATCGTGAGGGATCCTGCCCGATTTTGTAGCTTTCACCGGCGCGCTGGTAACCCTGTTTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATCCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGcc                                                                                                                                                                 <  2:132170/234‑1 (MQ=255)
                                                                       ccgccAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGCCCATCATCGCTTGTACATTCCGCGAGCATTCCTGCCGCGATGCGGATCCCGCGCCAc                                                                                                                                                               >  2:568741/1‑239 (MQ=255)
                                                                       ccgccAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCCCGATGCAGATCCCCGGCCa                                                                                                                                                                >  1:914248/1‑238 (MQ=255)
                                                                          ccAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACCATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGCCCATCATCCCTTGTACATGCCCCGGGCCTTCCTGCCCCGGTGCAAAACCCGCGCCAcg                                                                                                                                                              >  2:999243/1‑236 (MQ=255)
                                                                             gTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCCGGCTGCGGTACACTTACCGCCTGTTGCTGACCATCATCGCTTGTACATGCCGCGCGCCTTCCTGCCGCGCTGCAGATCCccacccacacctt                                                                                                                                                          >  2:291752/1‑233 (MQ=255)
                                                                                   ccTGGCTAAATCGAGACGGATCCTGCACGAGTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGCTTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTCATTTCAACTATAGGGCCGTTACATACCGCAGGCCTCGGTACACTTACCGCCTGTTGCTTGCCATCCTCGCTTGTACAAGACGCGCGCATTCCAGCCCCGGAGCAGAGCCCACCCCACTACTTACGCaat                                                                                                                                                   >  2:544393/1‑237 (MQ=255)
                                                                                      ggttAAATGGAGGCGGTTCCTTCCCGATTTTGTAGTTTTTACCGTCGGTCTGGTAACCCTGGTTGGCCTGCGGGTACGGTTGGTCGAAACGCGGGTCCGCCCCGTTAATTCCAACTATAGGGCCGTTACATACCGGAGGCTGCGGTACACTTACCGTTTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTgcg                                                                                                                                                <  2:931103/236‑1 (MQ=255)
                                                                                        ttAAATCGAGACGGTTCCGGCTCGATTTTGTAGCTTTTACCGTCGGGGTGGTAATCCTGATTTCCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCCGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTgcgc                                                                                                                                               <  1:918795/237‑1 (MQ=255)
                                                                                        cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTAGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACCCTTACCGTCGGTTGCTGACCATAATCCCTTGTAAATGCCCCACGCATTCCCGCCCCGTTCGACATCCCCCGCCACTCCTTACCCAATgcgc                                                                                                                                               >  2:1082889/1‑238 (MQ=255)
                                                                                        cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGGAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAAGCCGCGAGCATTCCTGCCCGGATGCAAACACCGAGACAATTCttcccaactcgcc                                                                                                                                               >  2:150991/1‑227 (MQ=255)
                                                                                        cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGAAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCCCCGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATACCTTGTACATGCCCCGCTCATTCCCGCCCCGTGGCAGATCCCGACCCCCTACTTaaccctcgccg                                                                                                                                               >  2:1006119/1‑228 (MQ=255)
                                                                                         tAAATCGGGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTGCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                                                                                                                             <  1:891441/239‑1 (MQ=255)
                                                                                          aaaccggggcGGATCCGGCACGATTTTGTAGCTTTTACCGTCGCGCGGGTAATCCGGATTGGCCGTCGCGTTAAGGGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                                                                                                                             <  1:291752/229‑1 (MQ=255)
                                                                                                 gaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatc                                                                                                                                                                                                                   <  2:302848/161‑1 (MQ=255)
                                                                                                 gaCGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatc                                                                                                                                                                                                                   >  1:302848/1‑161 (MQ=255)
                                                                                                           gCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                                                                                                                                                                                      <  2:336536/132‑1 (MQ=255)
                                                                                                           gCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                                                                                                                                                                                                      >  1:336536/1‑132 (MQ=255)
                                                                                                            cACGATTTTGTAGCTTTTCCCGTGGCGTTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGGGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATCCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGa                                                                                                                           <  1:666444/238‑1 (MQ=255)
                                                                                                                 ttttGTGCCTTTTCCCGTGGCGGTGGTAACCCTTGTTTGCCGTCGCGTTCAGTTGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATGGGGCCGTTACAATCCGCAGGCGGCGGTCCCCTTACCGTCGGTTGCGGCCCATCATCGCTTGTCCATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTCCACGCTTTTCGACAACAt                                                                                                                     <  2:228753/239‑1 (MQ=255)
                                                                                                                  tttGAAGCTTTTCCCAGCCGGTTGGTAAACTCGATTTGCCGCCGCTTCCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACCCTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAt                                                                                                                     <  1:238613/238‑1 (MQ=255)
                                                                                                                  ttgggaGCTTTTACCGTCGCTCTGGTAATCCTGATTTGCCGTCGCGTTCAGGGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAtt                                                                                                                    <  1:354856/234‑1 (MQ=255)
                                                                                                                    tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGGTGCAGATCCCCAGCCACTGCTTACCCATTGCGCACCCCTTACACGCTTTTCGAAAAAAttt                                                                                                                   >  2:367956/1‑238 (MQ=255)
                                                                                                                         cTTTTACCGTCGAGCTGGTAATCCTGTTTTGGCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtgt                                                                                                             <  1:150991/239‑1 (MQ=255)
                                                                                                                           tttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTCCCGTCTGTTGCCGACCATCCTCGCTTGTACATGCCGCGGGCATTCCAGCCGCGCTGCAGATCCCCAACCCCCCCTTAACCATTGGGCCCCCCTTCACCGCCTTTCGCCAACATTTTCCTgggg                                                                                                           >  2:741733/1‑239 (MQ=255)
                                                                                                                           tttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcg                                                                                                                                                                                   <  2:76158/167‑1 (MQ=255)
                                                                                                                           tttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcg                                                                                                                                                                                   >  1:76158/1‑167 (MQ=255)
                                                                                                                           tttACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCCATGCTTCCGCATTGGGCCCCCCTTCCACCCTTTTCGACAAAATTTTcccgggg                                                                                                           >  2:325975/1‑234 (MQ=255)
                                                                                                                                   cgcgCTGGTAATCCGTATTTGCCGTCGCTTTCAGTGGTTCGAACCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTCCCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGAt                                                                                                   <  1:999243/239‑1 (MQ=255)
                                                                                                                                       cTGGTAATCCTGATTTGCCGTCGCGTTCAGGGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                                                                                                              <  1:467999/224‑1 (MQ=255)
                                                                                                                                                  gATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTTCCGCGATGCAGATCCCGAGCCACTGCTTACCCATTGCGCCCCCCTTACACGCTTTTCGACAACATTTTCCTCTGGGTGGGGATTGAAATTAAAAAccc                                                                                    >  2:118864/1‑239 (MQ=255)
                                                                                                                                                    tttGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGc                                                                                                                                                      >  1:356408/1‑171 (MQ=255)
                                                                                                                                                    tttGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCCCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCCTTGCGCCCCTCTTACACCCTTTTCGCCAACACTTTCCTGTGGGGGTGGATTGACATTCAAATcccaa                                                                                  >  2:1031448/1‑237 (MQ=255)
                                                                                                                                                               gTTCAGTGGTTCGAAACGGGGGTCCGCCCCGCTAATTGCAACTATAGGGCCTTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTCCATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCAt                                                                        <  1:1082889/238‑1 (MQ=255)
                                                                                                                                                                                gcgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGa                                                      <  1:367956/239‑1 (MQ=255)
                                                                                                                                                                                gcgGGTCCCCCCCGCTAATTTCAACTATAGGGCCTTTACATACCGCAGGCTGCGGTACCCTTACCGTCTGTTGCTGACCATCAGCGCTTGTACGTGCAGCGAGCATTCCTGCCGCTATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGa                                                      <  1:878907/239‑1 (MQ=255)
                                                                                                                                                                                 cgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACa                                                                                         >  1:729669/1‑203 (MQ=255)
                                                                                                                                                                                 cgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGACCCACTGCTTACGCATTTCGCACCTCTTACACGCTTTTCGACAACCTTTTCCTTTGGGTGTGGTTTGACATTACAATCCCCAACCCAACCACCAGCCAAATTAgctcgca                                                      >  1:1028105/1‑233 (MQ=255)
                                                                                                                                                                                     tCCGCCCCTGTAACTTCAACTTATGGGCCTTTACCTCCCGCAGGCTGCGGTTCCCTTACCGCCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGAtcct                                                  <  1:344027/238‑1 (MQ=255)
                                                                                                                                                                                      ccGCCCCGGTAATTTCAACTTATGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGAtcctc                                                 <  1:741733/238‑1 (MQ=255)
                                                                                                                                                                                      ccGCCCCGCTAATTTCAACTATAGGGCCGTTCCCTCCCGCAGGCTGCGGTACCTTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGAtcctc                                                 <  1:794773/238‑1 (MQ=255)
                                                                                                                                                                                        gCCCCGCTAATTTCAACTTTGGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCAt                                              <  1:56474/239‑1 (MQ=255)
                                                                                                                                                                                              cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTACGCCCGACCCCCCATAACTg                                         >  1:1073075/1‑238 (MQ=255)
                                                                                                                                                                                              cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAAACCCGAACCCCGCCCTCAGCACAATTAGGCCCGAcccccacaaacga                                         >  2:1088690/1‑228 (MQ=255)
                                                                                                                                                                                                    tCAACTATGGGGCCGTTACCTCCCGAAGGCTGCGGTACACTTACCGGCTGTTGCTGACCATCACCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTg                                  <  1:325975/239‑1 (MQ=255)
                                                                                                                                                                                                           taGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGCGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCgg                           >  1:71089/1‑239 (MQ=255)
                                                                                                                                                                                                                   ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGCCATTACAATCCCGAACCCAGCCATCCGCACAATTAGCCCCGGCCCTCCAAAACTGACCAGGGGGAGCGGAAccccca                   >  2:571418/1‑237 (MQ=255)
                                                                                                                                                                                                                   ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATc                                                                                      >  1:939496/1‑172 (MQ=255)
                                                                                                                                                                                                                   ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAAtt                                                              >  1:912770/1‑196 (MQ=255)
                                                                                                                                                                                                                   ttACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCACCCATCCGCACAATTAGCGCCGATCCACCATAAACGACCCGGGGGAGCGGAGCCCCTa                   >  2:702926/1‑239 (MQ=255)
                                                                                                                                                                                                                       agaCCGCAGGCTGGTGTACATTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAAc               <  1:550185/237‑1 (MQ=255)
                                                                                                                                                                                                                        tACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATACGCCCAATTAGCGCCGATCCTCCTTACCt                                          >  1:591987/1‑211 (MQ=255)
                                                                                                                                                                                                                             cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCa                                                                             <  1:828454/171‑1 (MQ=255)
                                                                                                                                                                                                                             cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCa                                                                             >  2:828454/1‑171 (MQ=255)
                                                                                                                                                                                                                             cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAg         >  1:650626/1‑239 (MQ=255)
                                                                                                                                                                                                                             cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAACCATCAGCCCCATTTGCGCCGATCCTCCATAACTAACCAAGGGGAGCGGAACCCCAACAACCAGCaa         >  2:955989/1‑238 (MQ=255)
                                                                                                                                                                                                                               ggcttgggTACCCTTTCCGTTTGTTGCCTACCCTCAACCGTTGTACATGCCGCGAGCATTCCTGGCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGaa       <  1:908821/233‑1 (MQ=255)
                                                                                                                                                                                                                                gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                                                                                                              >  2:825957/1‑135 (MQ=255)
                                                                                                                                                                                                                                gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                                                                                                              <  1:825957/135‑1 (MQ=255)
                                                                                                                                                                                                                                gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAAc               >  1:158843/1‑230 (MQ=255)
                                                                                                                                                                                                                                   gaggTACCCTTACCGTCGTTTGCTGACCATCATCCCTTGTACATGCCGCGAGCATTCCTGCCGGGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATAc    <  1:1031448/236‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGTACACTTCCCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCg  <  1:1088690/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACATCTTACACGCTTTTCGACACACTTTTCCTGTGGGTATGGATTGCCATTACAATCCCGAACCCAGCCATCAGACCAATTAACGCCGATCCTCCCAAAATGACCAATGGAGACGGAAACCCAAAAAACGGCAAAAAAcac  >  2:550185/1‑237 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                        
CGTCAGCTGTGTCGGATCGAATGCTTCGCCAGAGGCCGTCAGGTTACTGCCGGGTTCGGCATCATAGATTGCCGCCAGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGACATTACAATCCCGAACCCAGCCATCAGCACAATTAGCGCCGATCCTCCATAACTGACCAGTGGGAGCGGAACCCCTACAACCGGCAGAATACCG  >  NC_000913/664870‑665336

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