Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,945 T→G 6.0% N1313T (AAC→ACC)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,074,9450TG 6.0% 230.2 / 6.7 67N1313T (AAC→ACC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base G (2/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCT  >  NC_000913/1074800‑1075092
                                                                                                                                                 |                                                                                                                                                   
aCGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTa                                                                                                                                                   >  2:748889/1‑148 (MQ=255)
  gTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                                                                 <  1:748889/148‑1 (MQ=255)
           tGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggca                                                                                                                                       <  3:160412/149‑1 (MQ=255)
            gAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcag                                                                                                                                      <  3:384014/149‑1 (MQ=255)
              ggTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagcg                                                                                                                                    <  4:122606/149‑1 (MQ=255)
              ggTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagc                                                                                                                                     <  2:613168/148‑1 (MQ=255)
               gTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagcg                                                                                                                                    >  2:648258/1‑148 (MQ=255)
                tATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcag                                                                                                                                            >  2:532883/1‑139 (MQ=255)
                tATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcag                                                                                                                                            <  1:532883/139‑1 (MQ=255)
                  tttttCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTa                                                                                                                                 >  2:352436/1‑148 (MQ=255)
                  tttttCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTAt                                                                                                                                >  3:467224/1‑149 (MQ=255)
                     ttCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTc                                                                                                                              <  2:24165/148‑1 (MQ=255)
                       caatcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTcac                                                                                                                            <  1:101677/148‑1 (MQ=255)
                          tcaaACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTc                                                                                                                        >  3:827058/1‑149 (MQ=255)
                             aaCGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAg                                                                                                                      >  2:175476/1‑148 (MQ=255)
                                 gTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGaa                                                                                                                  >  2:203905/1‑148 (MQ=255)
                                     aaaaCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                                                                     <  3:300455/141‑1 (MQ=255)
                                     aaaaCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                                                                     >  4:300455/1‑141 (MQ=255)
                                       aaCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCg                                                                                                            <  2:493661/148‑1 (MQ=255)
                                       aaCGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCg                                                                                                            <  1:352436/148‑1 (MQ=255)
                                          gTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGatat                                                                                                        <  4:22923/149‑1 (MQ=255)
                                               gCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGc                                                                                                   >  4:257992/1‑149 (MQ=255)
                                                          ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAg                                                                                        >  4:339793/1‑149 (MQ=255)
                                                            gAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGa                                                                                      <  3:257992/149‑1 (MQ=255)
                                                              ggggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGat                                                                                     >  1:139719/1‑148 (MQ=255)
                                                               gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGatat                                                                                   >  4:567098/1‑149 (MQ=255)
                                                                ggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGatat                                                                                   <  1:648258/148‑1 (MQ=255)
                                                                      gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTc                                                                            >  3:356605/1‑149 (MQ=255)
                                                                          aCGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGcc                                                                         >  1:97739/1‑148 (MQ=255)
                                                                             ttAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                                                                     <  3:469203/149‑1 (MQ=255)
                                                                                   aTGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaa                                                                >  2:681113/1‑148 (MQ=255)
                                                                                   aTGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaa                                                                >  1:224359/1‑148 (MQ=255)
                                                                                   aTGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaa                                                                >  1:140692/1‑148 (MQ=255)
                                                                                   aTGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaa                                                                >  1:448672/1‑148 (MQ=255)
                                                                                         gaggagGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGc                                                         >  3:491584/1‑149 (MQ=255)
                                                                                               gTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaa                                                   >  4:10151/1‑149 (MQ=255)
                                                                                               gTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaa                                                   >  4:439912/1‑149 (MQ=255)
                                                                                                   tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACa                                                >  2:382683/1‑148 (MQ=255)
                                                                                                     aCATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAAtt                                             <  4:581333/149‑1 (MQ=255)
                                                                                                       aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtg                                            <  1:415623/148‑1 (MQ=255)
                                                                                                       aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtg                                            <  2:139719/148‑1 (MQ=255)
                                                                                                         cctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGc                                         <  3:339793/149‑1 (MQ=255)
                                                                                                          ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGc                                         <  1:175476/148‑1 (MQ=255)
                                                                                                          ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGc                                         <  1:77430/148‑1 (MQ=255)
                                                                                                                  aCCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCcgcg                                >  3:189498/1‑149 (MQ=255)
                                                                                                                  aCCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCcgc                                 >  1:17901/1‑148 (MQ=255)
                                                                                                                        ttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCa                           >  2:588600/1‑148 (MQ=255)
                                                                                                                         tAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          <  2:310599/148‑1 (MQ=255)
                                                                                                                         tAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCATTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          <  1:588600/148‑1 (MQ=255)
                                                                                                                           aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACt                        <  1:681113/148‑1 (MQ=255)
                                                                                                                             cTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                     >  4:291235/1‑149 (MQ=255)
                                                                                                                             cTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                     >  3:225396/1‑149 (MQ=255)
                                                                                                                              tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                                                                        >  2:697625/1‑97 (MQ=255)
                                                                                                                              tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                                                                        <  1:697625/97‑1 (MQ=255)
                                                                                                                              tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                     >  2:403652/1‑148 (MQ=255)
                                                                                                                               ataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCtt                    <  1:403652/148‑1 (MQ=255)
                                                                                                                               ataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc                   <  4:827058/149‑1 (MQ=255)
                                                                                                                                   tCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcacac               >  4:443794/1‑149 (MQ=255)
                                                                                                                                     aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACaa             >  3:733239/1‑149 (MQ=255)
                                                                                                                                     aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACaa             >  4:422253/1‑149 (MQ=255)
                                                                                                                                      ttAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACaa             >  2:158444/1‑148 (MQ=255)
                                                                                                                                         cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCa                                       <  3:690084/118‑1 (MQ=255)
                                                                                                                                         cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCa                                       >  4:690084/2‑119 (MQ=255)
                                                                                                                                          gCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATgcgc        >  3:114757/1‑149 (MQ=255)
                                                                                                                                            tGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGcgcg       >  1:123731/1‑148 (MQ=255)
                                                                                                                                                gTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCt  >  3:787063/1‑149 (MQ=255)
                                                                                                                                                gTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCt  >  4:869484/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCT  >  NC_000913/1074800‑1075092

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