Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,799,634 A→G 6.7% intergenic (+21/+16) 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,6340AG 6.7% 160.9 / 3.6 60intergenic (+21/+16)alaE/ygaCalanine exporter, alanine‑inducible, stress‑responsive/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base G (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 = 2.84e-01

CATCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCA  >  NC_000913/2799485‑2799781
                                                                                                                                                     |                                                                                                                                                   
cATCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCa                                                                                                                                                     >  1:1061584/1‑150 (MQ=255)
         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:532303/1‑150 (MQ=255)
                                                                                                   gTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTc                                                  >  2:455331/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)
                                                                                                                                                     |                                                                                                                                                   
CATCAGATTATGGCGGCGGTCAGTTCAAACATCGTTGTTTCGATGTTGATGGGGGCGGTTTATGGCTACTTCCTCGATTATTGCCGCCGACTGTTTAAAGTCAGCCGTTACCAGCAGGTAAAAGCCTGACTCTTCTTTTCGCGACTGGCATCGCCAGTCGCGTCTTTAACTGGCTTCGCCAAATAAACCATTCAAATAACGTTCAAGCGCAACACGCGAACTAAAGCCATGCGGGATACCGTAATGCTCGTGTCGCTCGCCTGCTAAATAGAGCGGAAATTGCTGGTAATGATCGCA  >  NC_000913/2799485‑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: