Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,791,876 C→T 6.0% P202S (CCG→TCG)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,791,8760CT 6.0% 218.4 / 3.5 66P202S (CCG→TCG) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (25/37);  new base T (2/2);  total (27/39)
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.82e-01

ACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAA  >  NC_000913/2791732‑2792021
                                                                                                                                                |                                                                                                                                                 
aCGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGc                                                                                                                                               >  1:190290/1‑149 (MQ=255)
 cGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCt                                                                                                                                              <  3:469753/149‑1 (MQ=255)
  gCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTg                                                                                                                                             <  3:293816/149‑1 (MQ=255)
    cGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTggg                                                                                                                                           >  1:666482/1‑149 (MQ=255)
       ggAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTa                                                                                                                                        <  2:795591/149‑1 (MQ=255)
        gAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTAt                                                                                                                                       >  1:796578/1‑149 (MQ=255)
        gAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTAt                                                                                                                                       >  4:458262/1‑149 (MQ=255)
        gAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTAt                                                                                                                                       >  2:92388/1‑149 (MQ=255)
              cccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACg                                                                                                                                 >  3:716286/1‑149 (MQ=255)
               ccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                <  1:531933/149‑1 (MQ=255)
               ccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                <  2:241891/149‑1 (MQ=255)
               ccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                <  2:796578/149‑1 (MQ=255)
               ccGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                <  4:29149/149‑1 (MQ=255)
                cggcggCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTg                                                                                                                               <  1:534624/149‑1 (MQ=255)
                cggcggCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTg                                                                                                                               <  2:117024/149‑1 (MQ=255)
                  gcggcgATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGt                                                                                                                             <  2:190290/149‑1 (MQ=255)
                   cggcgATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTc                                                                                                                            <  4:74791/149‑1 (MQ=255)
                          gatTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGtt                                                                                                                     >  2:81/1‑149 (MQ=255)
                           atTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTc                                                                                                                    <  2:476616/149‑1 (MQ=255)
                                ccGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTcggcgg                                                                                                               >  2:394686/1‑149 (MQ=255)
                                  gCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAAc                                                                                                                                  >  3:36937/1‑128 (MQ=255)
                                  gCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAAc                                                                                                                                  <  4:36937/128‑1 (MQ=255)
                                  gCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCggcgggc                                                                                                             >  1:216855/1‑149 (MQ=255)
                                     aaGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACGGGTTCGGCGGGcgcg                                                                                                          >  3:795076/1‑149 (MQ=255)
                                     aaGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGcgcg                                                                                                          <  1:27494/149‑1 (MQ=255)
                                             ccGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTaa                                                                                                  <  3:614229/149‑1 (MQ=255)
                                                gcgcTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGa                                                                                               >  2:236702/1‑149 (MQ=255)
                                                 cgcTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGaa                                                                                              <  1:92388/149‑1 (MQ=255)
                                                 cgcTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGaa                                                                                              <  1:46866/149‑1 (MQ=255)
                                                  gcTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAAc                                                                                             <  2:744763/149‑1 (MQ=255)
                                                   cTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACt                                                                                            <  2:216855/149‑1 (MQ=255)
                                                      gcagcaGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                >  4:538940/1‑110 (MQ=255)
                                                      gcagcaGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGt                                                                                                                                <  3:538940/110‑1 (MQ=255)
                                                      gcagcaGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGAc                                                                                         <  3:378291/149‑1 (MQ=255)
                                                        agcagGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCa                                                                                       <  4:391011/149‑1 (MQ=255)
                                                            ggCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTaa                                                                                   >  1:770178/1‑149 (MQ=255)
                                                             gCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAAc                                                                                  <  3:47200/149‑1 (MQ=255)
                                                                 cACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGc                                                                              <  2:666482/149‑1 (MQ=255)
                                                                  aCCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCt                                                                             >  4:429983/1‑149 (MQ=255)
                                                                    cATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTgg                                                                           >  3:239436/1‑149 (MQ=255)
                                                                                 cccGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTc                                                              <  1:236702/149‑1 (MQ=255)
                                                                                   cGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGt                                                            >  2:183866/1‑149 (MQ=255)
                                                                                    gCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTTCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGtt                                                           >  3:290401/1‑149 (MQ=255)
                                                                                    gCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTTCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGtt                                                           >  3:619610/1‑149 (MQ=255)
                                                                                     ccAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGttt                                                          <  1:394686/149‑1 (MQ=255)
                                                                                          cAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCgg                                                     >  3:357000/1‑149 (MQ=255)
                                                                                            gACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTTCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGtt                                                   <  4:619610/149‑1 (MQ=255)
                                                                                            gACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTTCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTAACGGtt                                                   <  4:290401/149‑1 (MQ=255)
                                                                                            gACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGtt                                                   >  4:802500/1‑149 (MQ=255)
                                                                                            gACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACCAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGtt                                                   <  3:458262/149‑1 (MQ=255)
                                                                                             aCGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTc                                                  <  1:790073/149‑1 (MQ=255)
                                                                                                ccGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGAc                                               <  1:196167/149‑1 (MQ=255)
                                                                                                  gTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGTGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCg                                             <  2:630225/149‑1 (MQ=255)
                                                                                                   ttCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGa                                            >  1:372110/1‑149 (MQ=255)
                                                                                                                 gctggcgGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGt                              >  3:81244/1‑149 (MQ=255)
                                                                                                                     gcgGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAAt                          <  2:75875/149‑1 (MQ=255)
                                                                                                                      cgGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATg                         <  4:357000/149‑1 (MQ=255)
                                                                                                                       ggAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGCAATTGGCCGCCAGTTAATgg                        >  3:658992/1‑149 (MQ=255)
                                                                                                                             ggCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGt                  >  2:1815/1‑149 (MQ=255)
                                                                                                                             ggCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGt                  >  1:61627/1‑149 (MQ=255)
                                                                                                                               cGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTgc                <  1:183866/149‑1 (MQ=255)
                                                                                                                                     gcgcgGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGaaa          <  4:239436/149‑1 (MQ=255)
                                                                                                                                         gggCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACa      <  4:716286/149‑1 (MQ=255)
                                                                                                                                         gggCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACa      >  4:258534/1‑149 (MQ=255)
                                                                                                                                          ggCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACAt     <  2:372110/149‑1 (MQ=255)
                                                                                                                                             gTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCaa  <  3:429983/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACGCCGTGGAACTTCCCGGCGGCGATGATTACCCGCAAAGCCGGTCCGGCGCTGGCAGCAGGCTGCACCATGGTGCTGAAGCCCGCCAGTCAGACGCCGTTCTCTGCGCTGGCGCTGGCGGAGCTGGCGATCCGCGCGGGCGTTCCGGCTGGGGTATTTAACGTGGTCACCGGTTCGGCGGGCGCGGTCGGTAACGAACTGACCAGTAACCCGCTGGTGCGCAAACTGTCGTTTACCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAA  >  NC_000913/2791732‑2792021

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