Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,554,731 C→T M313I (ATG→ATA maeA ← malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,554,7310CT98.9% 305.9 / ‑4.7 90M313I (ATG→ATAmaeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (50/39);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGT  >  NC_000913/1554635‑1554827
                                                                                                |                                                                                                
cATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCCCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGattatt                                                                                               <  1:1886087/100‑1 (MQ=255)
   cTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTAGCGCTGGGTCTGGGAGATGATTATTTCg                                                                                            >  1:2057750/1‑100 (MQ=255)
     tGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGc                                                                                          <  1:1078725/100‑1 (MQ=255)
     tGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGc                                                                                          <  2:2771255/100‑1 (MQ=255)
      gtcagtcaGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCa                                                                                         <  2:1684326/100‑1 (MQ=255)
        cagtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAAt                                                                                       >  1:434715/1‑100 (MQ=255)
        cagtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAAt                                                                                       >  2:199964/1‑100 (MQ=255)
         agtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATg                                                                                      <  2:2604064/100‑1 (MQ=255)
         agtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATg                                                                                      >  2:937413/1‑100 (MQ=255)
          gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGc                                                                                     >  2:1502018/1‑100 (MQ=255)
          gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGc                                                                                     >  1:2183103/1‑100 (MQ=255)
          gtcagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCGTAATACTTCGCGCTGGGTCTGGGAGATGAATATTTCGGCAATGc                                                                                     >  1:133553/1‑100 (MQ=255)
            cagCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCg                                                                                   <  2:2293267/100‑1 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGc                                                                                  >  1:1375883/1‑100 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGc                                                                                  <  1:593429/100‑1 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGc                                                                                  <  1:1009380/100‑1 (MQ=255)
             agCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGc                                                                                  >  1:866348/1‑100 (MQ=255)
              gCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 >  2:1934203/1‑100 (MQ=255)
              gCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCa                                                                                 >  2:1202196/1‑100 (MQ=255)
               cAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCAt                                                                                >  2:371917/1‑100 (MQ=255)
                    cAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGc                                                                           <  2:557756/100‑1 (MQ=255)
                     aaaGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCt                                                                          >  2:2570096/1‑100 (MQ=255)
                      aaGCGATCGACCATAAAGACTTTCCGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTg                                                                         >  2:1983822/1‑100 (MQ=255)
                      aaGCGATCGACCATAAAGACTTTCCGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTg                                                                         >  2:1983414/1‑100 (MQ=255)
                       aGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGa                                                                        <  1:2331441/100‑1 (MQ=255)
                        ncGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGaa                                                                       >  1:805925/2‑100 (MQ=255)
                           aTCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCt                                                                    <  2:1161284/100‑1 (MQ=255)
                           aTCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCt                                                                    >  2:1284887/1‑100 (MQ=255)
                              gACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTgcg                                                                 >  2:1859440/1‑100 (MQ=255)
                                 cATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCa                                                              <  2:1664972/100‑1 (MQ=255)
                                  aTAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCaa                                                             >  2:2894803/1‑100 (MQ=255)
                                  aTAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCaa                                                             >  2:1395540/1‑100 (MQ=255)
                                  aTAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCaa                                                             >  1:1501007/1‑100 (MQ=255)
                                   tAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAg                                                            >  2:234514/1‑100 (MQ=255)
                                   tAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAg                                                            <  2:2030675/100‑1 (MQ=255)
                                    aaaGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAgg                                                           >  2:1589054/1‑100 (MQ=255)
                                    aaaGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAgg                                                           <  2:1898436/100‑1 (MQ=255)
                                    aaaGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAgg                                                           <  1:1802502/100‑1 (MQ=255)
                                    aaaGACTTTCTGCCGCGCCGCTTCCTCGCTAAACCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACGTGCGCCAAgg                                                           >  1:1906844/1‑100 (MQ=255)
                                       gACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAg                                                        <  1:2159452/100‑1 (MQ=255)
                                        aCTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAgt                                                       >  1:1694231/1‑99 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      <  2:1127730/100‑1 (MQ=255)
                                         cTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAc                                                      <  2:977636/100‑1 (MQ=255)
                                          tttCTGCCTCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAg                                                        <  1:952075/97‑1 (MQ=255)
                                              tGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttt                                                 >  1:205841/1‑100 (MQ=255)
                                              tGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttt                                                 <  2:346534/100‑1 (MQ=255)
                                               gccgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAtttt                                                >  2:183587/1‑100 (MQ=255)
                                                ccgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttttt                                               >  1:445796/1‑100 (MQ=255)
                                                ccgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAttttt                                               <  2:1617043/100‑1 (MQ=255)
                                                 cgcgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGAtttttt                                              <  1:2608811/100‑1 (MQ=255)
                                                   cggcgcTTCCTCGCTTAATCCTTCGCGCTGGCGCTGGGAGATGATTATTTCGGCAATGCCGGATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                            <  1:1385608/97‑1 (MQ=255)
                                                   cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                            <  1:1817914/100‑1 (MQ=255)
                                                   cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtc                                            <  1:247906/100‑1 (MQ=255)
                                                   cgccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGAAGATTTTTTtc                                            <  1:2500373/100‑1 (MQ=255)
                                                    gccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtct                                           >  1:2878392/1‑100 (MQ=255)
                                                     ccgcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTtctc                                          >  1:1668636/1‑100 (MQ=255)
                                                       gcTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGc                                        >  2:2768208/1‑100 (MQ=255)
                                                             tCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACt                                  >  1:1281485/1‑100 (MQ=255)
                                                                cTTAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTg                                 <  1:1281756/98‑1 (MQ=255)
                                                                 ttAATCCTGCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGacc                              <  2:1375902/100‑1 (MQ=255)
                                                                  tAATCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGacc                              <  2:434720/99‑1 (MQ=255)
                                                                  tAAGCCTTCGTGCTGGGTCTGGGAGATGATTATTCCGGCAATGCCGCATCCCGCTCAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGacc                              <  1:2419247/99‑1 (MQ=255)
                                                                   aaTCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGaccac                            <  1:1934175/100‑1 (MQ=255)
                                                                     tCCTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCt                          >  2:2596338/1‑100 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCACGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                        >  2:1343167/1‑100 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTgn                        <  2:805931/100‑2 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                        >  2:2127340/1‑100 (MQ=255)
                                                                       cTTCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGc                        >  2:2127399/1‑100 (MQ=255)
                                                                        ttCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGcc                       <  1:2010393/100‑1 (MQ=255)
                                                                        ttCGCGCTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGcc                       <  1:1501994/100‑1 (MQ=255)
                                                                           gcgcTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgcg                    >  2:1791166/1‑100 (MQ=255)
                                                                           gcgcTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgcg                    <  1:914678/100‑1 (MQ=255)
                                                                             gcTGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCcgcgcg                  <  2:1627202/100‑1 (MQ=255)
                                                                               tGGGTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGc                <  1:526669/100‑1 (MQ=255)
                                                                                gggTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCt               >  2:340511/1‑100 (MQ=255)
                                                                                  gTCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtg             >  2:536291/1‑100 (MQ=255)
                                                                                   tCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgc            <  2:117313/100‑1 (MQ=255)
                                                                                   tCTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgc            >  1:489003/1‑100 (MQ=255)
                                                                                    cTGGGAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGCTTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgct           >  1:671842/1‑100 (MQ=255)
                                                                                      gggAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTtgctgc         >  2:320915/1‑100 (MQ=255)
                                                                                       ggAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCg        <  1:371912/100‑1 (MQ=255)
                                                                                       ggAGATGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCg        >  2:2025245/1‑100 (MQ=255)
                                                                                        gagaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGa       >  2:2650870/1‑100 (MQ=255)
                                                                                         agaTGATTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGAt      <  2:2878430/100‑1 (MQ=255)
                                                                                          gatgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATc     >  2:1024057/1‑100 (MQ=255)
                                                                                          gatgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATc     <  2:995077/100‑1 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  2:425196/1‑100 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  2:2323180/1‑100 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  1:1569345/1‑100 (MQ=255)
                                                                                            tgatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAg   >  2:1660358/1‑100 (MQ=255)
                                                                                             gatTATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAgt  >  1:1383932/1‑100 (MQ=255)
                                                                                                |                                                                                                
CATCTTGTCAGTCAGCAAGCCAAAGCGATCGACCATAAAGACTTTCTGCCGCGCCGCTTCCTCGCTTAATCCTTCGCGCTGGGTCTGGGAGATGATCATTTCGGCAATGCCGCATCCCGCTGAACCTGCGCCAAGGAAGACGATTTTTTTCTCGCTTAACTGACCACCTGCCGCGCGGCTTGCTGCGATCAGT  >  NC_000913/1554635‑1554827

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: