Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,555,783 A→G 6.9% intergenic (‑114/+43) maeA ← / ← sra malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,7830AG 6.9% 260.8 / 7.9 87intergenic (‑114/+43)maeA/sramalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base G (3/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

AAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTCCA  >  NC_000913/1555637‑1555931
                                                                                                                                                  |                                                                                                                                                    
aaGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCTGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAgg                                                                                                                                                    >  2:200852/1‑149 (MQ=255)
     aaCGCTGTTTTTTTGTTTTTGGGTCCATGGCACGCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGGCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATg                                                                                                                                               <  1:6965/149‑1 (MQ=255)
       cGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGgt                                                                                                                                             <  4:25248/149‑1 (MQ=255)
        gCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgtg                                                                                                                                            <  4:97288/149‑1 (MQ=255)
        gCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGACAGGGCGTCGCACTCAAGTATAAAGCAGATAAAAACAAAAACACAATTGCGCAGGCAATGgtg                                                                                                                                            >  5:208299/1‑149 (MQ=255)
         cTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgtgt                                                                                                                                           <  2:139217/149‑1 (MQ=255)
           gTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGtt                                                                                                                                          >  8:73006/1‑148 (MQ=255)
             ttttttGTTTTTGGTTACATTTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTATATCGTCTTTGCATTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGCTGttaa                                                                                                                                        >  6:313559/1‑146 (MQ=255)
              tttttgtttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGTGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGGCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGttttt                                                                                                                                       <  6:218223/148‑3 (MQ=255)
                tttgtttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATc                                                                                                                                     <  2:26450/148‑1 (MQ=255)
                      tttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAATCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                               >  7:251951/1‑148 (MQ=255)
                      tttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                               >  5:320732/1‑148 (MQ=255)
                       ttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTg                                                                                                                              <  3:167045/148‑1 (MQ=255)
                                tGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATg                                                                                                                    <  3:12434/149‑1 (MQ=255)
                                    actcactcTTTTTTGAATATCCATCCCTGTGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTGGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGtt                                                                                                                 <  5:330825/148‑1 (MQ=255)
                                    actcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTa                                                                                                                <  5:41668/149‑1 (MQ=255)
                                     ctcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAt                                                                                                               <  8:144104/149‑1 (MQ=255)
                                     ctcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAt                                                                                                               >  2:89605/1‑149 (MQ=255)
                                     ctcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAt                                                                                                               <  2:341345/149‑1 (MQ=255)
                                      tcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAt                                                                                                               >  3:292976/1‑148 (MQ=255)
                                      tcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCACTGAGGACTGATGGTTAt                                                                                                               >  1:64849/1‑148 (MQ=255)
                                         ctctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAAt                                                                                                           >  2:289282/1‑149 (MQ=255)
                                         ctctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAAt                                                                                                           >  5:149291/1‑149 (MQ=255)
                                          tctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATc                                                                                                                                     >  1:212592/1‑122 (MQ=255)
                                          tctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATc                                                                                                                                     <  2:212592/122‑1 (MQ=255)
                                              ttttGAATATCCATCCCTGGGTGGCTTTTATCGTCTTTGTTTTACCGCCAGGGCGTCGGCCTCTAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTa                                                                                                         <  1:358613/146‑1 (MQ=255)
                                              ttttGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATAGAGGACTGATGGTTATGAATTAc                                                                                                        <  8:135365/147‑1 (MQ=255)
                                                 tGAATATCCATCCCTGGCGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtt                                                                                                                                          >  7:201353/1‑110 (MQ=255)
                                                 tGAATATCCATCCCTGGCGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtt                                                                                                                                          <  8:201353/110‑1 (MQ=255)
                                                  gAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACtttt                                                                                                    >  7:93639/1‑147 (MQ=255)
                                                   aaTATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAg                                                                                                 >  7:28487/1‑149 (MQ=255)
                                                      atccatccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCgg                                                                                              >  5:266353/1‑149 (MQ=255)
                                                       tccatccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCggg                                                                                             <  7:339563/149‑1 (MQ=255)
                                                          atccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCg                                                                                          >  6:25761/1‑149 (MQ=255)
                                                            cccTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGtt                                                                                        >  6:195629/1‑149 (MQ=255)
                                                               tGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTAAGCTGGGCGTTTTc                                                                                     >  5:256443/1‑149 (MQ=255)
                                                                  ggggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTTCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATTAATTACTTTTCAGCGGGGCGTTTTCTGc                                                                                  <  1:184097/149‑1 (MQ=255)
                                                                  ggggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCACAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTTATGGTTATGAATTACTTTTAAGCGGGGCGATTTCTGc                                                                                  >  2:344737/1‑149 (MQ=255)
                                                                    ggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTAGAAAGCAGATAAAAACACACACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCg                                                                                <  5:95056/149‑1 (MQ=255)
                                                                     gCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCgg                                                                               >  3:310094/1‑149 (MQ=255)
                                                                       ttttATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt                                                                             >  4:175318/1‑149 (MQ=255)
                                                                           aTCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTggg                                                                          <  6:208299/148‑1 (MQ=255)
                                                                             cGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTa                                                                       >  6:212136/1‑149 (MQ=255)
                                                                               tCTTTGCTTTACCGCGAGGGCGTCGGGCTCAAGTATAAAGCAGATAAAAACAAAACCACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGTGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTAtt                                                                     <  3:266315/149‑1 (MQ=255)
                                                                                 tttGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                                   >  4:27106/1‑149 (MQ=255)
                                                                                 tttGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                                   >  8:263634/1‑149 (MQ=255)
                                                                                  ttGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTCATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTtac                                                                  >  3:179177/1‑149 (MQ=255)
                                                                                   tGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                                   <  3:6876/147‑1 (MQ=255)
                                                                                     cTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTtacta                                                                <  7:263634/148‑1 (MQ=255)
                                                                                      tttACCGCCAGGGGGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACg                                                              >  6:89782/1‑149 (MQ=255)
                                                                                      tttACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACg                                                              >  8:130162/1‑149 (MQ=255)
                                                                                         aCCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCt                                                            <  5:25761/148‑1 (MQ=255)
                                                                                         aCCGCCAGGGAGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCt                                                            <  1:301629/148‑1 (MQ=255)
                                                                                           cGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGa                                                         <  1:344737/149‑1 (MQ=255)
                                                                                            gCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAt                                                        <  6:149291/149‑1 (MQ=255)
                                                                                            gCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAt                                                        >  6:229335/1‑149 (MQ=255)
                                                                                              cAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGAttt                                                      <  8:28487/149‑1 (MQ=255)
                                                                                                  gCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCa                                                  >  5:166403/1‑149 (MQ=255)
                                                                                                   cGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCAc                                                 <  6:320732/149‑1 (MQ=255)
                                                                                                     tCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGt                                                    >  8:331735/1‑144 (MQ=255)
                                                                                                     tCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCt                                               >  7:82265/1‑149 (MQ=255)
                                                                                                     tCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGGTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGt                                                    <  7:331735/144‑1 (MQ=255)
                                                                                                        gCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCg                                            <  1:200852/149‑1 (MQ=255)
                                                                                                        gCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCg                                            <  3:175318/149‑1 (MQ=255)
                                                                                                         ccTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                           <  3:27106/149‑1 (MQ=255)
                                                                                                               gTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTa                                     >  7:112328/1‑149 (MQ=255)
                                                                                                                tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTAt                                    >  7:100210/1‑149 (MQ=255)
                                                                                                                tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTAt                                    >  3:69500/1‑149 (MQ=255)
                                                                                                                  taAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTAttt                                  <  1:289282/149‑1 (MQ=255)
                                                                                                                  taAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTAttt                                  >  4:220945/1‑149 (MQ=255)
                                                                                                                   aaaGCAGATAAAAACAAAAACACCATTGCGAAGGCAATGGTGATTAATCGTAATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTTGATTATTTACTACTCTGGAATTGTCACCTTCGGTAACTAtttt                                 >  6:227535/1‑149 (MQ=255)
                                                                                                                    aaGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTg                                <  5:265967/149‑1 (MQ=255)
                                                                                                                     aGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTg                                                                            <  6:218970/104‑1 (MQ=255)
                                                                                                                       cAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTgcgc                             <  8:112328/149‑1 (MQ=255)
                                                                                                                        aGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCt                            <  8:77534/149‑1 (MQ=255)
                                                                                                                             aaaacaaaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGtt                        <  5:195629/148‑1 (MQ=255)
                                                                                                                                gattaaacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGa                     >  3:185704/5‑148 (MQ=255)
                                                                                                                                   aaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAtt                 <  5:212136/149‑1 (MQ=255)
                                                                                                                                      acacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa              <  5:268588/149‑1 (MQ=255)
                                                                                                                                      acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTAGGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa              >  7:61829/1‑149 (MQ=255)
                                                                                                                                       cacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTAt             >  7:43310/1‑149 (MQ=255)
                                                                                                                                          cATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgt          >  2:138195/1‑149 (MQ=255)
                                                                                                                                           aTTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtc         <  4:179177/149‑1 (MQ=255)
                                                                                                                                           aTTGCGCAGGCAATGGGTTTTAATCGTCATTGAGGACGGATGGGTCTGAAGTACTTGTCAGCGGGGCGTTTTCTGCCGGTGGGGGTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtc         <  4:18275/149‑1 (MQ=255)
                                                                                                                                            ttGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtcc        <  8:195103/149‑1 (MQ=255)
                                                                                                                                               cgcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGTGAGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccag      <  5:227535/148‑1 (MQ=255)
                                                                                                                                                gcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccagtc    >  7:248930/1‑149 (MQ=255)
                                                                                                                                                  aGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccagtcca  <  8:82265/149‑1 (MQ=255)
                                                                                                                                                  aGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTAAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtcaagtcc   >  1:354571/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAGCGAACGCTGTTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTCCA  >  NC_000913/1555637‑1555931

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: