Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,799,635 T→A 6.7% intergenic (+22/+15) alaE → / ← ygaC alanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,799,6350TA 6.7% 154.2 / 2.3 60intergenic (+22/+15)alaE/ygaCalanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (31/23);  new base A (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

ATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCA  >  NC_000913/2799494‑2799781
                                                                                                                                             |                                                                                                                                                  
aTGGCGGCGGTCAGATCCAACATCGTTGTAGCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTc                                                                                                                                            <  2:869667/150‑1 (MQ=255)
 tGGCGGCGGTCAGTGCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTcg                                                                                                                                           <  2:519002/150‑1 (MQ=255)
  ggcggcGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcgc                                                                                                                                          <  2:1234967/150‑1 (MQ=255)
  ggcggcGGGCAGTTCAAACATCGTTGATTCGATGTTGATGGGGGCGGGTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCGATGCCAGTcgc                                                                                                                                          <  1:775264/150‑1 (MQ=255)
   gcggcgGTCATTTCAAACATCGTTGTTTCGATTTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTcgcg                                                                                                                                         <  1:399643/150‑1 (MQ=255)
   gcggcgGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGAAAAATCCTGACTCTTCTTTTCGCGACTGGCATCGCACGTcgcg                                                                                                                                         >  2:662453/1‑150 (MQ=255)
         gTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGCCTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGCCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTa                                                                                                                                   >  2:261534/1‑150 (MQ=255)
         gTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTa                                                                                                                                   >  2:952484/1‑150 (MQ=255)
           cAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAAc                                                                                                                                 <  2:597757/150‑1 (MQ=255)
            aGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGATACTACCTCGATTATTGCCGCCGACTGTTTAAAGTCATCCTTTACCAGCAGGTATAAGCCTGACTCTTCTTTTCGCGACTGTCATCGCCCGTCGCGTCTTTAACt                                                                                                                                >  2:214544/1‑150 (MQ=255)
             gTTCAAACATCGTTGTTTCGCTGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTg                                                                                                                               <  1:103352/150‑1 (MQ=255)
                cAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCt                                                                                                                            >  2:65914/1‑150 (MQ=255)
                          tgttTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATa                                                                                                                  <  1:493331/150‑1 (MQ=255)
                            tttCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTAGCCAAATaaa                                                                                                                >  1:1092314/1‑150 (MQ=255)
                                gATGTTGATGGTGGCGGTTTATGGCTACTTGCTCGTTTATTGCGGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCGTGACGCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCTAATAAACCAt                                                                                                            <  2:321361/150‑1 (MQ=255)
                                 ttGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCAtt                                                                                                           <  2:465903/149‑1 (MQ=255)
                                   gTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCa                                                                                                         <  2:1254270/150‑1 (MQ=255)
                                   gTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGGCAGCCGTTACCACCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCa                                                                                                         <  1:653160/150‑1 (MQ=255)
                                     tGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTGTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTGaaa                                                                                                       >  1:134319/1‑150 (MQ=255)
                                       aTGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATa                                                                                                     >  1:1065154/1‑150 (MQ=255)
                                             gCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTAAGCCGTTACCAGCAGTTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCAAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTc                                                                                               <  2:1192060/150‑1 (MQ=255)
                                                gTTTATGGCTACTACCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAg                                                                                            <  2:179002/150‑1 (MQ=255)
                                                 tttATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAgc                                                                                           >  1:616419/1‑150 (MQ=255)
                                                    aTGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACAAGCAGGTAAAAGCCTGACTAGTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGGTTAGCCAAATTAACCATTCAAATAACTTTAAAgcgaa                                                                                        >  1:211381/1‑148 (MQ=255)
                                                         tACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACAcg                                                                                   >  2:1232301/1‑150 (MQ=255)
                                                           cTTCCTCGATTATTGCGGCGGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACAcgcg                                                                                 <  2:1179208/150‑1 (MQ=255)
                                                              ccTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAAc                                                                              <  1:1248224/150‑1 (MQ=255)
                                                               cTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGATAAAAGCCTGACTCTTCTTTTCGCGACTGGCATAGCCAGTCGCGGCTTTAACAGGCTTCGCCAAATTAACCCTTCAAATTACGTTGAAGCGCCCCAcgcggcct                                                                             >  1:1267581/1‑146 (MQ=255)
                                                                   actattGCCGCCGACTGTTTAAATTCAGCCGTTACCAGCAGGTAAAAGCCTGTCTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAg                                                                         <  2:135059/148‑1 (MQ=255)
                                                                    ttattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCAGCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGc                                                                        <  1:806345/150‑1 (MQ=255)
                                                                     tattGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGcc                                                                       <  2:145076/150‑1 (MQ=255)
                                                                         gccgccGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCTCCAGTCGCGTCTTTAACTGGCGTCGCCAAATAACCCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGc                                                                   >  1:881408/1‑150 (MQ=255)
                                                                         gccgccGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGc                                                                   >  1:1092742/1‑150 (MQ=255)
                                                                         gccgccGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGGCTTTAACTGGCTTCGCCAAATAAACCATTCATATTACGTTCAAGCGCAACACCCGAACTAAAGCCATGc                                                                   >  1:814726/1‑150 (MQ=255)
                                                                            gccgACTGTTTAAAGTCAGCAGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTAAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCggg                                                                <  2:475552/150‑1 (MQ=255)
                                                                             ccgACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCAGGCGGGa                                                               >  2:1287867/1‑150 (MQ=255)
                                                                               gACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCGTCGCCAAATAACCCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATa                                                             >  1:653978/1‑150 (MQ=255)
                                                                                     ttAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGGGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCCACACGCGAACTAAAGCCATGCGGGATACCGTaa                                                       >  2:480574/1‑150 (MQ=255)
                                                                                          gTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCACGCGCAACACGCGAACTTAAGCCATGCGGGAGACCGTAATGCTc                                                  >  2:53203/1‑150 (MQ=255)
                                                                                          gTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTc                                                  >  2:455331/1‑150 (MQ=255)
                                                                                          gTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTc                                                  >  2:532303/1‑150 (MQ=255)
                                                                                                cGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcg                                            >  1:48445/1‑150 (MQ=255)
                                                                                                   tACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgctc                                         <  2:1112765/150‑1 (MQ=255)
                                                                                                   tACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgctc                                         >  1:677282/1‑150 (MQ=255)
                                                                                                   tACCAGCAGCTCAAAGCCTGACTCTTCTTTTCTCGACTGGCATCGCGAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATCACGTTCCACCGCAACCCGCGAACTAACGCCAAGCGGGATACCGTAATGCTCGTGtcgctc                                         <  2:188223/150‑1 (MQ=255)
                                                                                                    aCCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgctcg                                        >  1:709213/1‑150 (MQ=255)
                                                                                                    aCCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGACAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGtcgctcg                                        <  1:810942/150‑1 (MQ=255)
                                                                                                                          cttcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt                  <  1:757513/150‑1 (MQ=255)
                                                                                                                          cttcttTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAAtt                  <  2:616419/150‑1 (MQ=255)
                                                                                                                              ttttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTg              <  1:261201/150‑1 (MQ=255)
                                                                                                                              ttttCGCGACTGGCATCGCCAGTCGCGGCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGTACGAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCGCGCCGGCAAAATAGAGCGGAAATTGGTg              >  1:1029596/1‑150 (MQ=255)
                                                                                                                               tttCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg             >  2:294491/1‑150 (MQ=255)
                                                                                                                               agacgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTgg             >  2:96263/4‑150 (MQ=255)
                                                                                                                                  cgcgACTGGCGATGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTaa          >  1:148338/1‑150 (MQ=255)
                                                                                                                                  cgcgACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCCAAATAGAGCGGAAATTGCTGGTat          >  2:895497/1‑149 (MQ=255)
                                                                                                                                   gcgACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAAt         >  2:3132/1‑150 (MQ=255)
                                                                                                                                   gcgACTGGCAGCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAAt         <  2:340704/150‑1 (MQ=255)
                                                                                                                                      aCTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGAt      >  1:1152309/1‑150 (MQ=255)
                                                                                                                                         ggCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGc   >  2:789321/1‑150 (MQ=255)
                                                                                                                                          gCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTATTGATCGCa  >  2:668115/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCA  >  NC_000913/2799494‑2799781

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