Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 1,542,119 A→C F158V (TTT→GTT)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,542,1190AC96.5% 296.7 / ‑1.8 86F158V (TTT→GTT) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (39/44);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAAACCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTT  >  NC_000913_3_bme_pgi/1541972‑1542263
                                                                                                                                                   |                                                                                                                                                
gATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCAGCGACAGGCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAAccc                                                                                                                                                >  1:646690/1‑150 (MQ=38)
 aTAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGCAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAAcccc                                                                                                                                               >  1:113008/1‑150 (MQ=255)
           ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCNGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     <  1:541685/150‑1 (MQ=255)
           ccGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     <  1:1035799/150‑1 (MQ=255)
             ggcggcGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCccgcc                                                                                                                                            <  3:27632/141‑1 (MQ=255)
             ggcggcGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCccgcc                                                                                                                                            >  4:27632/1‑141 (MQ=255)
              gcggcgTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGc                                                                                                                                  >  1:755319/1‑150 (MQ=255)
              gcggcgTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGGCCATACGGTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGc                                                                                                                                  >  1:755615/1‑150 (MQ=255)
                 gcgTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACt                                                                                                                               >  1:1157539/1‑150 (MQ=255)
                   gTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     >  3:20611/1‑142 (MQ=255)
                   gTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     <  4:20611/142‑1 (MQ=255)
                    tAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGCAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgt                                                                                                                            >  3:507023/1‑150 (MQ=255)
                     aaGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     <  1:756088/140‑1 (MQ=255)
                     aaGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCa                                                                                                                                     >  2:756088/1‑140 (MQ=255)
                     aaGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgtt                                                                                                                           >  1:928365/1‑150 (MQ=255)
                      aGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgttt                                                                                                                          <  2:928365/150‑1 (MQ=255)
                       gAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgtttg                                                                                                                         >  4:678072/1‑150 (MQ=255)
                        aaaCCATCGACATCGCCGGGATCGGCGAGAAACCGTCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgtttgt                                                                                                                        >  2:237803/1‑150 (MQ=255)
                        aaaCCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgtttgt                                                                                                                        <  1:497124/150‑1 (MQ=255)
                         aaCCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTa                                                                                                                       <  1:920463/150‑1 (MQ=255)
                               cGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAg                                                                                                                 <  1:21488/150‑1 (MQ=255)
                                gACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAgg                                                                                                                >  1:153609/1‑150 (MQ=255)
                                     cGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACtt                                                                                                           >  1:215814/1‑150 (MQ=255)
                                      gCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACttt                                                                                                          <  2:484877/150‑1 (MQ=255)
                                          ggATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCg                                                                                                      <  1:646410/150‑1 (MQ=255)
                                               ggCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCg                                                                                                 >  1:807453/1‑150 (MQ=255)
                                               ggCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCg                                                                                                 >  3:367576/1‑150 (MQ=255)
                                                gCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCggtgtta                                                                                                >  3:640631/1‑145 (MQ=255)
                                                 cGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGtt                                                                                               >  4:451420/1‑150 (MQ=255)
                                                 cGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGtt                                                                                               <  2:755615/150‑1 (MQ=255)
                                                 cGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGtt                                                                                               <  2:755319/150‑1 (MQ=255)
                                                  gagaAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTc                                                                                              <  3:575804/150‑1 (MQ=255)
                                                  gagaAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTc                                                                                              >  2:1238872/1‑150 (MQ=255)
                                                   agaAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCa                                                                                             <  4:367576/150‑1 (MQ=255)
                                                   agaAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCa                                                                                             >  2:739207/1‑150 (MQ=255)
                                                       aCCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGc                                                                                                       <  1:857610/136‑1 (MQ=255)
                                                       aCCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGc                                                                                                       >  2:857610/1‑136 (MQ=255)
                                                       aCCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAAt                                                                                         >  2:537828/1‑150 (MQ=255)
                                                       aCCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAAt                                                                                         >  1:1135002/1‑150 (MQ=255)
                                                        ccGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATc                                                                                        >  4:133483/1‑150 (MQ=255)
                                                              aCGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGcc                                                                                  <  3:678072/150‑1 (MQ=255)
                                                               cgcgATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgttt                                                                                                                          >  1:401514/1‑109 (MQ=255)
                                                               cgcgATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACttgttt                                                                                                                          <  2:401514/109‑1 (MQ=255)
                                                                  gATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCa                                                                                >  1:709018/1‑148 (MQ=255)
                                                                  gATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCa                                                                                <  2:709018/148‑1 (MQ=255)
                                                                   aTCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGc                                                                             <  1:1172100/150‑1 (MQ=255)
                                                                   aTCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGc                                                                             <  1:806240/150‑1 (MQ=255)
                                                                       gggCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGt                                                                         >  1:338628/1‑150 (MQ=255)
                                                                           cGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGc                                                                                                       >  2:849466/1‑116 (MQ=255)
                                                                           cGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGc                                                                                                       <  1:849466/116‑1 (MQ=255)
                                                                           cGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCg                                                                     >  1:168383/1‑150 (MQ=255)
                                                                             tAGGTTTTGATGGTCCAGGCGTGAGTGGGGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGGGCGGATAACCCCGCCGCGCCCACGGACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGa                                                                   <  4:548895/150‑1 (MQ=255)
                                                                              aGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACtt                                                                                                                              >  3:499789/1‑90 (MQ=255)
                                                                              aGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACtt                                                                                                                              <  4:499789/90‑1 (MQ=255)
                                                                               ggTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATc                                                                 <  4:507023/150‑1 (MQ=255)
                                                                               ggTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATc                                                                 <  3:451420/150‑1 (MQ=255)
                                                                                    tGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCtgt                                                            <  4:213134/150‑1 (MQ=255)
                                                                                      aTGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAAccc                                                                                                                                                <  4:215274/64‑1 (MQ=255)
                                                                                          tCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCt                                                      >  4:218887/1‑150 (MQ=255)
                                                                                             aGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGa                                                   <  2:1135002/150‑1 (MQ=255)
                                                                                             aGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGa                                                   <  1:875111/150‑1 (MQ=255)
                                                                                              gACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCtt                                                     >  3:722660/1‑147 (MQ=255)
                                                                                              gACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCtt                                                     <  4:722660/147‑1 (MQ=255)
                                                                                              gACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAg                                                  <  2:1157539/150‑1 (MQ=255)
                                                                                                        gcggcAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTACTTGAGGGCTTCgcgc                                        <  2:646690/150‑1 (MQ=255)
                                                                                                          ggcAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTa                                                                                                                       >  1:226067/1‑69 (MQ=255)
                                                                                                          ggcAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTa                                                                                                                       <  2:226067/69‑1 (MQ=255)
                                                                                                            cAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCcaca                                    <  1:1238872/150‑1 (MQ=255)
                                                                                                            cAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCcaca                                    <  1:739207/150‑1 (MQ=255)
                                                                                                             aaTCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGa                                                                                      <  3:610895/99‑1 (MQ=255)
                                                                                                             aaTCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGa                                                                                      >  4:610895/1‑99 (MQ=255)
                                                                                                                  gCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCaa                              >  4:684262/1‑150 (MQ=255)
                                                                                                                   cTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCtgt                                                            >  3:320390/1‑119 (MQ=255)
                                                                                                                   cTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCtgt                                                            <  4:320390/119‑1 (MQ=255)
                                                                                                                   cTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAggg                                                >  4:746588/1‑131 (MQ=255)
                                                                                                                   cTGGTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAggg                                                <  3:746588/131‑1 (MQ=255)
                                                                                                                     ggTTTAGTTCCTGCCAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCa                           <  2:113008/150‑1 (MQ=255)
                                                                                                                        ttAGTTCCTGCCAGTTGGAGCGGATAAACCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATc                                                                 <  4:521252/109‑1 (MQ=255)
                                                                                                                        ttAGTTCCTGCCAGTTGGAGCGGATAAACCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATc                                                                 >  3:521252/1‑109 (MQ=255)
                                                                                                                                  ccAGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAAt              <  1:537828/150‑1 (MQ=255)
                                                                                                                                    aGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAgg                                                 <  1:920561/113‑1 (MQ=255)
                                                                                                                                    aGTTGGAGCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAgg                                                 >  2:920561/1‑113 (MQ=255)
                                                                                                                                           gCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTc                                <  1:935892/123‑1 (MQ=255)
                                                                                                                                           gCGGATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTc                                >  2:935892/1‑123 (MQ=255)
                                                                                                                                             ggATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTAtt   <  2:153609/150‑1 (MQ=255)
                                                                                                                                              gATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGc                                                       <  2:285877/97‑1 (MQ=255)
                                                                                                                                              gATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGc                                                       >  1:285877/1‑97 (MQ=255)
                                                                                                                                              gATAACCCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTAttt  >  4:140288/1‑150 (MQ=255)
                                                                                                                                              gATAAACCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTAttt  >  4:237603/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATAACGCGTTCCGGCGGCGTAAGAAACCATCGACATCGCCGGGATCGGCGAGAAACCGGCAACGCGATCCGGGCCGTAGGTTTTGATGGTCCAGACGTTAGCGGCGGCAATCAGCTGGTTTAGTTCCTGCCAGTTGGAGCGGATAAACCCGCCGCGCCCACGCACTTGTTTGTAGCTCAGGCACTTTTGCGGATCGTTCATAATCGATGCCCACGCCAGTACCGGATCGCTGTGTTGCTTGAGGGCTTCGCGCCACAGTTCAATCAGTCGTTTACGAATGAGCGGGTATTT  >  NC_000913_3_bme_pgi/1541972‑1542263

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