Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 646,110 A→G P157P (CCT→CCC)  citT ← citrate/succinate antiporter; citrate carrier

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,1100AG98.4% 209.9 / ‑4.0 61P157P (CCT→CCC) citTcitrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (32/28);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

GAACTCAGACTGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCA  >  NC_000913/645964‑646257
                                                                                                                                                  |                                                                                                                                                   
gAACTCAGACTGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGgtg                                                                                                                                                 >  1:358816/1‑151 (MQ=255)
 aaCTCAGACTGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGgtg                                                                                                                                                 <  2:551607/150‑1 (MQ=255)
   cTCAGACTGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTa                                                                                                                                               >  2:134408/1‑150 (MQ=255)
     cAGACTGGTCCTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTaaa                                                                                                                                             <  2:608361/150‑1 (MQ=255)
       gACTGGTACTAATGAACATCATCCACTTCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACgg                                                                                                                                          <  1:908073/151‑1 (MQ=255)
          tGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAAACGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGc                                                                                                                                       >  2:215359/1‑151 (MQ=255)
           ggTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGcc                                                                                                                                      >  1:48282/1‑151 (MQ=255)
                 aaTGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCcagcag                                                                                                                                 <  2:366636/150‑1 (MQ=255)
                 aaTGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCcagcag                                                                                                                                 <  2:774769/150‑1 (MQ=255)
                 aaTGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCcagc                                                                                                                                   >  2:797489/1‑148 (MQ=255)
                   tGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAAt                                                                                                                              >  1:285796/1‑151 (MQ=255)
                      ccATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTc                                                                                                                           <  2:121673/151‑1 (MQ=255)
                        atcatcCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCg                                                                                                                          >  1:482540/1‑150 (MQ=255)
                            tcCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAAt                                                                                                                     >  1:265613/1‑151 (MQ=255)
                               acaTCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGAc                                                                                                                  <  2:640664/151‑1 (MQ=255)
                                   cAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACaa                                                                                                                <  2:185795/149‑1 (MQ=255)
                                      aTAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCa                                                                                                           <  1:120365/151‑1 (MQ=255)
                                          ccgccAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATaa                                                                                                        >  1:909146/1‑150 (MQ=255)
                                            gccAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAAccc                                                                                                     >  1:84386/1‑151 (MQ=255)
                                               aaTACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAAc                                                                                                  >  1:463978/1‑151 (MQ=255)
                                                aTACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACg                                                                                                 >  1:832081/1‑151 (MQ=255)
                                                   cgacgCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCa                                                                                              <  1:810114/151‑1 (MQ=255)
                                                          cggacggaTCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAgcgtgcg                                                                                       >  1:747337/1‑151 (MQ=255)
                                                           ggccggaTCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTCAATGACCGGAAAAACCGTACCCCCGGTACGGGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAgcgtgcgt                                                                                      <  2:868718/148‑1 (MQ=255)
                                                             acggaTCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGttt                                                                                    >  1:212085/1‑151 (MQ=255)
                                                             acggaTCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGttt                                                                                    >  1:189637/1‑151 (MQ=255)
                                                               agaTCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAAACCAACGTCAGCGTGCGTTTGc                                                                                  >  2:576373/2‑151 (MQ=255)
                                                                   cGTTCGGGAATGATTTAAACAACGGCGGCGGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCAt                                                                              <  1:393507/151‑1 (MQ=255)
                                                                      tCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGa                                                                            >  1:61848/1‑150 (MQ=255)
                                                                       cGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAAt                                                                          <  2:394678/151‑1 (MQ=255)
                                                                        gggAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAAtt                                                                         >  1:739888/1‑151 (MQ=255)
                                                                         ggAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAAttt                                                                        >  2:868386/1‑151 (MQ=255)
                                                                              gATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCAcc                                                                    <  1:27552/150‑1 (MQ=255)
                                                                                 ttAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAgg                                                                 <  2:212059/150‑1 (MQ=255)
                                                                                 ttAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGAACGCGCGGTGTTGGCGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCAc                                                                     >  2:91103/1‑146 (MQ=255)
                                                                                  tAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAgg                                                                 >  1:744016/1‑149 (MQ=255)
                                                                                   acacaacGGCGGCAGGTTTTTAATGACCGGAAAAACCGTGACCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCGGAATGACGATAAAGACAATCGCATAACCCAACGTCAGCATGCGTTTGCCGATGAATATCACCAGGa                                                                <  2:430265/147‑1 (MQ=255)
                                                                                    aacaacGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAAc                                                                                                  >  2:266744/1‑114 (MQ=255)
                                                                                    aacaacGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAAc                                                                                                  <  1:266770/114‑1 (MQ=255)
                                                                                    aacaacGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGaaaa                                                             <  1:125999/151‑1 (MQ=255)
                                                                                      caacGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAgg                                                           >  2:521021/1‑151 (MQ=255)
                                                                                         cggcggCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCa                                                        <  1:105905/151‑1 (MQ=255)
                                                                                                 ggTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGcgacg                                                <  2:91907/151‑1 (MQ=255)
                                                                                                 ggTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGcgacg                                                <  2:549557/151‑1 (MQ=255)
                                                                                                     tttAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGAcc                                             <  2:417612/150‑1 (MQ=255)
                                                                                                       tAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTa                                           <  1:98499/150‑1 (MQ=255)
                                                                                                        aaTGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAAc                                         >  1:694952/1‑151 (MQ=255)
                                                                                                          tGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAAccc                                       >  2:231347/1‑151 (MQ=255)
                                                                                                           gACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAAccc                                       <  2:203637/150‑1 (MQ=255)
                                                                                                           gACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCATTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCg                                      >  1:430994/1‑151 (MQ=255)
                                                                                                             ccGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGa                                    <  2:683391/151‑1 (MQ=255)
                                                                                                                 aaaaaCCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACt                                >  2:828924/1‑151 (MQ=255)
                                                                                                                   aaaCCGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTc                              >  2:195866/1‑151 (MQ=255)
                                                                                                                      ccGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCAt                            >  2:294154/1‑150 (MQ=255)
                                                                                                                       cGTACCCCCGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATAc                          >  1:211353/1‑151 (MQ=255)
                                                                                                                             cccGGTACGCGCGGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAAt                    >  1:886345/1‑151 (MQ=255)
                                                                                                                                     gcgcgGTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAaaat            <  2:528135/151‑1 (MQ=255)
                                                                                                                                         ggTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAaaataaat        <  2:187515/151‑1 (MQ=255)
                                                                                                                                         ggTGTTGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAaaataaat        >  1:53718/1‑151 (MQ=255)
                                                                                                                                              tGGAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATgcgc    <  2:336086/150‑1 (MQ=255)
                                                                                                                                               ggAGGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCa  <  1:899660/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAACTCAGACTGGTACTAATGACCATCATCCACATCAAATAGCCGCCAATACGACGCGCGGACGGATCGTTCGGGAATGATTTAAACAACGGCGGCAGGTTTTTAATGACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCA  >  NC_000913/645964‑646257

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