Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,674,896 G→C 5.4% L556L (CTG→CTC)  PP_4137 → outer membrane siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,674,8960GC 5.4% 240.5 / 2.6 74L556L (CTG→CTC) PP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (39/31);  new base C (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GTACTACTCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCGGGGTCAA  >  NC_002947/4674749‑4675042
                                                                                                                                                     |                                                                                                                                                   
gTACTACTCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGa                                                                                                                                                    <  2:116296/149‑1 (MQ=255)
 tactacTCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAg                                                                                                                                                   <  1:253050/149‑1 (MQ=255)
      ctCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTgc                                                                                                                                              >  5:299640/1‑149 (MQ=255)
       tCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTGGCGCAAGCTCAGCTTg                                                                                                                                               >  2:89617/1‑149 (MQ=255)
       tCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTgcg                                                                                                                                             >  5:310497/1‑149 (MQ=255)
          cAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTGGCGCAAGCTCAGCTTgcca                                                                                                                                            <  1:89617/149‑3 (MQ=255)
          cAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTCAGCTTgccagt                                                                                                                                          <  1:296335/149‑5 (MQ=255)
             gTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGc                                                                                                                                       <  2:252338/149‑1 (MQ=255)
                   cACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTacac                                                                                                                                 >  2:71341/1‑149 (MQ=255)
                   cACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GAAAGCTGAGCTTGCGCCAGGCCTaca                                                                                                                                  >  8:22874/1‑148 (MQ=255)
                    aCGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACAcc                                                                                                                                >  4:330122/1‑149 (MQ=255)
                        aCTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAg                                                                                                                                         >  8:142052/1‑136 (MQ=255)
                        aCTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAg                                                                                                                                         <  7:142052/136‑1 (MQ=255)
                        aCTGCTGGCCGGTGAGGTAGTGCCGAACTTCTCCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑CCAAGCTGAGCTTGCGCCAGGCCTACACCTTCa                                                                                                                            <  2:149351/149‑1 (MQ=255)
                            ctgGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGa                                                                                                                        <  6:86732/149‑1 (MQ=255)
                             tgGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGAc                                                                                                                       >  7:122860/1‑149 (MQ=255)
                                ccGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCAGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTc                                                                                                                    <  4:107564/149‑1 (MQ=255)
                                    tGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCAc                                                                                                                 >  4:144646/1‑148 (MQ=255)
                                           gTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGtga                                                                                                         <  6:299640/149‑1 (MQ=255)
                                           gTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGtga                                                                                                         <  1:71341/149‑1 (MQ=255)
                                            tGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTAcc                                                                                                             <  5:125138/144‑1 (MQ=255)
                                            tGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTAcc                                                                                                             >  6:125138/1‑144 (MQ=255)
                                               cGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATgac                                                                                                     >  6:71515/1‑149 (MQ=255)
                                                    ttCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCa                                                                                                                            <  8:33709/121‑1 (MQ=255)
                                                    ttCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCa                                                                                                                            >  7:33709/1‑121 (MQ=255)
                                                     tCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAg                                                                                                                                         <  3:182355/107‑1 (MQ=255)
                                                     tCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAg                                                                                                                                         >  4:182355/1‑107 (MQ=255)
                                                     tCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACg                                                                                               >  1:211585/1‑149 (MQ=255)
                                                      cACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGacac                                                                                                <  5:71515/147‑1 (MQ=255)
                                                      cACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGt                                                                                              >  5:309323/1‑149 (MQ=255)
                                                       tccTCGGAATTCAACGCCCGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGtt                                                                                             <  4:130997/148‑1 (MQ=255)
                                                       aCCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGtt                                                                                             <  3:296168/149‑1 (MQ=255)
                                                         cTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTg                                                                                           >  2:89866/1‑149 (MQ=255)
                                                         cTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTg                                                                                           >  4:41762/1‑149 (MQ=255)
                                                             gAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGa                                                                                      <  3:144646/149‑1 (MQ=255)
                                                                ttCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACaa                                                                                   >  6:288404/1‑149 (MQ=255)
                                                                   aaCGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAAcc                                                                                 <  6:310497/148‑1 (MQ=255)
                                                                   aaCGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCg                                                                                >  7:122965/1‑149 (MQ=255)
                                                                   aaCGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCg                                                                                >  3:291460/1‑149 (MQ=255)
                                                                    aCGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGt                                                                               <  3:330122/149‑1 (MQ=255)
                                                                      gCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCt                                                                             >  2:301835/1‑149 (MQ=255)
                                                                      gCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCt                                                                             >  6:244056/1‑149 (MQ=255)
                                                                         aGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGcc                                                                          <  8:122860/149‑1 (MQ=255)
                                                                                  gTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATccc                                                                 >  7:342321/1‑149 (MQ=255)
                                                                                    gCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCa                                                               <  5:288404/149‑1 (MQ=255)
                                                                                      aCCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGGCACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATg                                                             >  4:163496/1‑149 (MQ=255)
                                                                                      aCCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATg                                                             <  8:342321/149‑1 (MQ=255)
                                                                                        cAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACACCCGTCTGCCCGGCATCCCCATGCa                                                           >  5:378616/1‑149 (MQ=255)
                                                                                           ggCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCacta                                                        <  6:309323/149‑1 (MQ=255)
                                                                                                   aGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAgg                                                >  8:273632/1‑149 (MQ=255)
                                                                                                    gCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCGGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGc                                               >  8:250148/1‑149 (MQ=255)
                                                                                                          cTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCt                                         >  8:179193/1‑149 (MQ=255)
                                                                                                           tCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTg                                        >  1:125372/1‑149 (MQ=255)
                                                                                                             gACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCt                                         >  1:167101/1‑146 (MQ=255)
                                                                                                             gACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCt                                         <  2:167101/146‑1 (MQ=255)
                                                                                                             gACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GAGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCTAGCTgcg                                      >  6:156436/1‑149 (MQ=255)
                                                                                                              aCAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTGTGCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTgcg                                      <  1:276022/149‑1 (MQ=255)
                                                                                                              aCAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTGTGCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTgcg                                      <  1:280720/149‑1 (MQ=255)
                                                                                                                  cACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTAt                                  >  7:192639/1‑148 (MQ=255)
                                                                                                                  cACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTAt                                  >  7:240079/1‑148 (MQ=255)
                                                                                                                  cACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTCG‑GCGACAACCGTCTGCACGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTAt                                  >  4:221798/1‑148 (MQ=255)
                                                                                                                         tgGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTgg                            <  4:291460/147‑1 (MQ=255)
                                                                                                                         tgGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGc                           <  7:312650/148‑1 (MQ=255)
                                                                                                                             aGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAAc                      >  5:215344/1‑149 (MQ=255)
                                                                                                                                  gcgGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCtt                 <  7:179193/149‑1 (MQ=255)
                                                                                                                                    gggCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCt               <  3:41762/149‑1 (MQ=255)
                                                                                                                                        aCTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGc           >  4:281949/1‑149 (MQ=255)
                                                                                                                                          tggcgCAAGCTCAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGc           >  7:239143/5‑149 (MQ=255)
                                                                                                                                          tG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCg         >  3:165096/1‑149 (MQ=255)
                                                                                                                                          tG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCg         >  2:188532/1‑149 (MQ=255)
                                                                                                                                           g‑‑gCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACTACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGaa                       <  3:212822/134‑1 (MQ=255)
                                                                                                                                           g‑‑gCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACTACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGaa                       >  4:212822/1‑134 (MQ=255)
                                                                                                                                               cAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCggg       <  8:122965/148‑1 (MQ=255)
                                                                                                                                                   cTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCGGGGTCaa  >  3:330279/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTACTACTCGCAAGTACGCCACGAACTGCTGGCCGTTGAGGTAGTGCCGAACTTCACCTCGGAATTCAACGCCAGCGCTACCGTGCACCAAGGCGTGGAAGCCGGCCTCGACAGCACCCTGTGGGAGCGCGCGGGCACTG‑‑GCAAGCTGAGCTTGCGCCAGGCCTACACCTTCAGCGACTTCCACTACCGTGATGACGACACGTTTG‑GCGACAACCGTCTGCCCGGCATCCCCATGCACTACTACCAGGCCGAGCTGCGCTATGACTGGCCGAACGGCTTCTATGCCGGGGTCAA  >  NC_002947/4674749‑4675042

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