Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,577,413 A→G 6.7% L322P (CTG→CCG)  maeB ← malic enzyme: putative oxidoreductase/phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,577,4130AG 6.7% 177.0 / 3.1 60L322P (CTG→CCG) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (32/24);  new base G (2/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ATCCTGATCGCCATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTTTCTTCACCAT  >  NC_000913/2577265‑2577561
                                                                                                                                                    |                                                                                                                                                    
aTCCTGATCGCCATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGTAAGATGAACGGGAAGcaca                                                                                                                                                      <  8:272689/149‑1 (MQ=255)
  ccTGATCGCCATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAgg                                                                                                                                                    >  1:199765/1‑149 (MQ=255)
         gCCATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTGCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACgttg                                                                                                                                              <  3:48307/148‑1 (MQ=255)
           cATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACgttgtt                                                                                                                                            >  6:372095/1‑148 (MQ=255)
            aTACGCTGAAGCCACCACTTCTCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCa                                                                                                                                          >  7:395004/1‑149 (MQ=255)
                    aaGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTc                                                                                                                                  <  7:110376/149‑1 (MQ=255)
                         accacTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATa                                                                                                                             <  7:352744/149‑1 (MQ=255)
                             cTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTATCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGt                                                                                                                           >  8:105623/1‑147 (MQ=255)
                                cGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGaa                                                                                                                      >  1:199819/1‑149 (MQ=255)
                                cGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGa                                                                                                                       <  2:121548/148‑1 (MQ=255)
                                   tctGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTc                                                                                                                          >  5:292384/1‑142 (MQ=255)
                                   tctGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTc                                                                                                                          <  6:292384/142‑1 (MQ=255)
                                      gTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCg                                                                                                                >  2:79404/1‑149 (MQ=255)
                                       ttCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCg                                                                                                                <  2:199765/148‑1 (MQ=255)
                                       taccGCATGGGCGAGGTCTGCAATCGCACGTTCCGGCGCCAGTTTGATCTCTTCGTTGATGGCGGTTGCGCCAACGTGCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGGTGTTCACCTGGTTCGGATAGTCAGAACGAc                                                                                                                  <  7:361330/144‑1 (MQ=255)
                                           gCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGc                                                                                                            >  8:198776/1‑148 (MQ=255)
                                             aTGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCaaa                                                                                                         <  5:372965/149‑1 (MQ=255)
                                                ggCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACGGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAatga                                                                                                      >  4:62360/1‑149 (MQ=255)
                                                      ttCTGCAATCGAACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACATGTTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCAt                                                                                                <  5:387596/149‑1 (MQ=255)
                                                        cTGCACTCGCACGTAACGACGCCAGTTTCATCTCATCGTTGATGGCGGTTGCGCCCACGTCCAGCGCGCCACGGAAGATTAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGGCAGAACGACCGGTGCAAATGATGGCATc                                                                                               <  7:357054/148‑1 (MQ=255)
                                                         tGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCg                                                                                             >  4:235727/1‑149 (MQ=255)
                                                          gCAATCGCACGTACCGCCGCCAGTTTCATCTATTCGTTTATGGCGGTTGCGCCAACTTCCAGCGAGCCACGGATGATGAAGGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGATCGCCCGGGGCAAATGATGGCATCCg                                                                                             >  6:125567/1‑148 (MQ=255)
                                                             aTCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACGGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACg                                                                                         >  8:280078/1‑149 (MQ=255)
                                                              tCGCACGTACCGACGCCAGTTTCATCTCTTCGTTGATGGCGTTTGCGCCAACGTCCAGCGCCCCACGGAAGATGAACGGGAAGAACAGGACGTTGTTCACCTGGTTCGTATAGTCATAACCACCGGTGCAAATGATGGCATCCGGACGc                                                                                        >  6:50748/1‑149 (MQ=255)
                                                                    gTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCActtct                                                                                  <  2:199819/149‑1 (MQ=255)
                                                                               aGTTTCAGCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGCCGGTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAgcgg                                                                       <  2:267225/149‑1 (MQ=255)
                                                                                gTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGCCCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAgcggc                                                                      >  5:330254/1‑149 (MQ=255)
                                                                                   tCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGATCGGGAAGCACAGGACGTTGTTCACATGGTTCGTATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGcggcgg                                                                    >  1:74054/1‑148 (MQ=255)
                                                                                    cATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGcggcgg                                                                    >  7:53067/1‑147 (MQ=255)
                                                                                      tctcTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAg                                                                 >  1:93905/1‑148 (MQ=255)
                                                                                         ctTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAAtt                                                             >  6:335972/1‑149 (MQ=255)
                                                                                         ctTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAAtt                                                             >  5:22064/1‑149 (MQ=255)
                                                                                          ttCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAAttt                                                            >  2:61483/1‑149 (MQ=255)
                                                                                             gTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttcc                                                          <  8:158564/148‑1 (MQ=255)
                                                                                              ttGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCCCAGGACGTTGTGCCCCTGGTTCGGAGAGGCAGAACGCCCGGTGCAAATGATGGCAGCCGGACGCACGTCTTTCGCCAGCGGCGGCAGAAGttccgg                                                        >  1:371429/1‑149 (MQ=255)
                                                                                              ttGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttccgg                                                        <  3:235727/149‑1 (MQ=255)
                                                                                                gATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAatgat                                                                                                     <  8:16664/102‑1 (MQ=255)
                                                                                                gATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAatgat                                                                                                     >  7:16664/1‑102 (MQ=255)
                                                                                                 aTGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACGGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttccggttc                                                     <  7:280078/149‑1 (MQ=255)
                                                                                                  tGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAGATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttccggttcc                                                    >  8:203185/1‑149 (MQ=255)
                                                                                                    gCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttccggttccg                                                   <  7:198776/148‑1 (MQ=255)
                                                                                                    gCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTGCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATttccggttccgg                                                  <  5:318528/149‑1 (MQ=255)
                                                                                                       gTTGCGCCAACGTCCAGCGCGCCACGGAAGATGACCGGGAAGCACAGGACGTTGTTCACCTGGTGCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGGTCCGGGtt                                               <  5:50748/149‑1 (MQ=255)
                                                                                                            gcCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCggcggc                                                                   >  7:50220/1‑124 (MQ=255)
                                                                                                            gcCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCggcggc                                                                   <  8:50220/124‑1 (MQ=255)
                                                                                                            gcCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTcgc                                            >  7:273475/1‑147 (MQ=255)
                                                                                                            gcCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTcgc                                            >  8:379332/1‑147 (MQ=255)
                                                                                                              cAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTcgccagc                                        >  3:312596/1‑149 (MQ=255)
                                                                                                                  gTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTcgccagcgc                                      >  7:140384/1‑147 (MQ=255)
                                                                                                                    ccAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTcgccagcgccag                                   >  6:83326/1‑148 (MQ=255)
                                                                                                                          cgcCACGGAAGATGAACGGGAAGCACGGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCAtt                            <  3:62360/149‑1 (MQ=255)
                                                                                                                          cgcCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCAtt                            <  4:312596/149‑1 (MQ=255)
                                                                                                                            ccACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCAt                             >  6:224072/1‑146 (MQ=255)
                                                                                                                            ccACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCAt                             <  5:224072/146‑1 (MQ=255)
                                                                                                                             cACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGgc                         >  2:45895/1‑149 (MQ=255)
                                                                                                                                ggAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTGCGCCAGCGCCAGGATCATTGGCGCa                      <  4:360608/149‑1 (MQ=255)
                                                                                                                                   aGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGa                   <  5:372095/149‑1 (MQ=255)
                                                                                                                                    gATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAg                  >  7:99679/1‑149 (MQ=255)
                                                                                                                                       gAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCa               >  8:142168/1‑149 (MQ=255)
                                                                                                                                           gggAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCAtttt           <  5:83326/149‑1 (MQ=255)
                                                                                                                                             gAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTttct         >  7:355163/1‑149 (MQ=255)
                                                                                                                                                 cacaGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTttctt        <  5:157426/146‑1 (MQ=255)
                                                                                                                                                 cacaGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTttctt        >  6:157426/1‑146 (MQ=255)
                                                                                                                                                   caGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTTTCTTCACCa   <  8:395004/149‑1 (MQ=255)
                                                                                                                                                    aGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAGATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTTTCTTCACCAt  <  7:203185/149‑1 (MQ=255)
                                                                                                                                                    aGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTTTCTTCACCAt  >  5:242165/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCCTGATCGCCATACGCTGAAGCCACCACTTCGCTCTGTTCCGCATGGGCGAGTTCTGCAATCGCACGTACCGCCGCCAGTTTCATCTCTTCGTTGATGGCGGTTGCGCCAACGTCCAGCGCGCCACGGAAGATGAACGGGAAGCACAGGACGTTGTTCACCTGGTTCGGATAGTCAGAACGACCGGTGCAAATGATGGCATCCGGACGCACTTCTTTCGCCAGCGGCGGCAGAATTTCCGGTTCCGGGTTCGCCAGCGCCAGGATCATTGGCGCACGAGCCATTTTCTTCACCAT  >  NC_000913/2577265‑2577561

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: