Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,228,476 C→T 5.3% L226F (CTT→TTT)  ygjJ → putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,228,4760CT 5.3% 239.2 / 2.1 76L226F (CTT→TTT) ygjJputative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base T (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTGGACGTACCGGGAAGTG  >  NC_000913/3228336‑3228623
                                                                                                                                            |                                                                                                                                                   
gCAACTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCg                                                                                                                                               >  6:212010/1‑147 (MQ=255)
 cAACTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAATGCAAAGCCTGGCTACATTTTGGCAATTACGATCGCTTTGACGTGAAATACACCTATCCGGCATTGATGTACGGTTACATCGAGTATGACCTTGCCGAAg                                                                                                                                            <  3:236352/149‑1 (MQ=255)
 cAACTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAg                                                                                                                                            <  2:348793/149‑1 (MQ=255)
  aaCTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAgg                                                                                                                                           >  4:296440/1‑149 (MQ=255)
  aaCTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAgg                                                                                                                                           <  7:45880/149‑1 (MQ=255)
   aCTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGc                                                                                                                                          <  5:212010/149‑1 (MQ=255)
     tttCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCAt                                                                                                                                        >  8:253096/1‑149 (MQ=255)
       tCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATc                                                                                                                                       <  5:28019/148‑1 (MQ=255)
          ttATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCt                                                                                                                                   >  1:288305/1‑149 (MQ=255)
           tataGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAAAACGATCGCTATGACGACAAATACACCTATACGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTg                                                                                                                                  >  2:33525/1‑149 (MQ=255)
            ataGCGTCACCGCGACTCCTGCCTGACATCATGATGAAAGCAAAGCGTGGCTACATTTTGGCAATTACGATCGCTATGACGACAATTACTCCTTTCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTg                                                                                                                                  <  4:68765/148‑1 (MQ=255)
             taGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGa                                                                                                                                <  2:288305/149‑1 (MQ=255)
             taGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGa                                                                                                                                <  3:12107/149‑1 (MQ=255)
                  tCACCGCGACTCCTTCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATCATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAAc                                                                                                                           <  6:128077/149‑1 (MQ=255)
                  tCACCGCGACTCCTGCTGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAAc                                                                                                                           <  3:258296/149‑1 (MQ=255)
                  tCACCGCGACTCCTGCCGGACATCATGATGATAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAAc                                                                                                                           <  4:16768/149‑1 (MQ=255)
                  tCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAAc                                                                                                                           <  7:253096/149‑1 (MQ=255)
                     ccGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCTCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGt                                                                                                                        <  7:330157/149‑1 (MQ=255)
                            tCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGaaa                                                                                                                 >  8:60834/1‑149 (MQ=255)
                              cTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAAt                                                                                                                >  1:222135/1‑148 (MQ=255)
                               tGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCa                                                                                                              >  4:186319/1‑149 (MQ=255)
                                 ccGGACCTCATGATGAACGCAACGCCTGGCTACAGTGTGGCAAATACGAGCGCTCGGCCGACAACTACACCTATGCGGCACTGATGACCGGTTACATCCTGTATGACCTTGCCGAAGGCATCTCCTGGATGAACGGTCTGGAAATCAcc                                                                                                            <  4:191434/149‑1 (MQ=255)
                                 ccGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCa                                                                                                              <  2:197490/147‑1 (MQ=255)
                                  cGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCAcc                                                                                                            >  1:86061/1‑148 (MQ=255)
                                   ggACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGa                                                                                                          <  5:125999/149‑1 (MQ=255)
                                       atcatGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGACCGGTTACATCCAGTATGACCTTTCAGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCTACgg                                                                                                       >  3:242389/1‑148 (MQ=255)
                                          atgatgAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGca                                                                                                     >  7:132026/1‑147 (MQ=255)
                                           tgatgaAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAg                                                                                                  >  2:263038/1‑149 (MQ=255)
                                               gaAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCTCTAAGACGACAAATACACCTAACCGGCAATGATGAACGGTTACATCCAGTATGACCTAGCCGAACGCATCACCTGGATGAACGGTCTCGAAATCACCGACGGCACAGGa                                                                                                >  8:317167/1‑147 (MQ=255)
                                               gaAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAg                                                                                                  >  1:68039/1‑145 (MQ=255)
                                               gaAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAg                                                                                                  <  2:68039/145‑1 (MQ=255)
                                               gaAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGa                                                                                                >  4:325632/1‑147 (MQ=255)
                                                  agcaaagcCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGct                                                                                           >  7:61832/1‑149 (MQ=255)
                                                    caaagcCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGctc                                                                                          >  7:348069/1‑148 (MQ=255)
                                                     aaagcCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTa                                                                                        >  3:231959/1‑149 (MQ=255)
                                                      aagcCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTAt                                                                                       >  2:248649/1‑149 (MQ=255)
                                                       agcCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTAtc                                                                                      <  3:296440/149‑1 (MQ=255)
                                                         ccTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTAtct                                                                                     >  3:327182/1‑148 (MQ=255)
                                                          cTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCa                                                                                   <  8:348069/149‑1 (MQ=255)
                                                             gCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACg                                                                                 <  3:186319/148‑1 (MQ=255)
                                                               tACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCAc                                                                                  >  7:289665/1‑145 (MQ=255)
                                                               tACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCAc                                                                                  <  8:289665/145‑1 (MQ=255)
                                                                  aTTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTg                                                                           >  7:179881/1‑149 (MQ=255)
                                                                  aTTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTg                                                                           >  4:342321/1‑149 (MQ=255)
                                                                   ttttGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACTTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTg                                                                           >  8:267546/1‑148 (MQ=255)
                                                                   ttttGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACTTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGc                                                                          >  8:267549/1‑149 (MQ=255)
                                                                    tttGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGc                                                                          >  8:242275/1‑148 (MQ=255)
                                                                       ggCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTaa                                                                       <  8:132026/148‑1 (MQ=255)
                                                                         cAATTACGATCGCTATGACGACAAATACCCCTATCCGGCAAAGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTc                                                                    <  8:196564/149‑1 (MQ=255)
                                                                             tACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTaa                                                                <  1:248649/149‑1 (MQ=255)
                                                                             tACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTaa                                                                <  1:263038/149‑1 (MQ=255)
                                                                               cGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACt                                                              >  6:289385/1‑149 (MQ=255)
                                                                                gATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGACTGGAAAGCACCGACGGCAAAGGACAGCTCTATCTCACCGGCCTGATAACTACTAACtt                                                             >  6:233318/1‑149 (MQ=255)
                                                                                 aTCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACttt                                                            <  1:104092/149‑1 (MQ=255)
                                                                                  tCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTg                                                           <  4:231959/149‑1 (MQ=255)
                                                                                  tCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTg                                                           <  8:179881/149‑1 (MQ=255)
                                                                                  tCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTg                                                           >  6:166716/1‑149 (MQ=255)
                                                                                    gCTATGACGACAAATGCACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACGTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGcc                                                         <  4:301087/149‑1 (MQ=255)
                                                                                      tATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTAtct                                                                                     >  6:93533/1‑119 (MQ=255)
                                                                                      tATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTAtct                                                                                     <  5:93533/119‑1 (MQ=255)
                                                                                          acgacAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACTTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTcg                                                    <  7:267546/148‑1 (MQ=255)
                                                                                          acgacAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACTTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTcg                                                    <  7:267549/148‑1 (MQ=255)
                                                                                               aaaTACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTg                                              <  7:60834/149‑1 (MQ=255)
                                                                                                     acCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCa                                        <  3:325632/149‑1 (MQ=255)
                                                                                                          tCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTgg                                                                                                                    >  3:320415/1‑68 (MQ=255)
                                                                                                          tCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTTACTTTGCCGCTCGCGCCTGGCACCATAc                                     <  4:327182/147‑1 (MQ=255)
                                                                                                            cGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAATGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGa                                 <  4:233625/149‑1 (MQ=255)
                                                                                                             ggCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATTAACGGGCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGTCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGAc                                <  4:242389/149‑1 (MQ=255)
                                                                                                                   gatgaACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCg                          >  2:335278/1‑149 (MQ=255)
                                                                                                                         cGGTTCCATCCAGTATGACCTTGCCGATGGCATCACCTGGATGAACGGGCTGGAAATCACCGTCGGCACAGGACTGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGc                    <  5:237988/149‑1 (MQ=255)
                                                                                                                          ggTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCt                   <  5:289385/149‑1 (MQ=255)
                                                                                                                            ttACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTg                  >  5:294187/1‑148 (MQ=255)
                                                                                                                                      tATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCTTACCGGACGCGCCGACGGGCTGGACGTACCggg       <  2:86061/149‑1 (MQ=255)
                                                                                                                                        tGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTCTCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTGGACGTACCGGGaa     <  5:233318/149‑1 (MQ=255)
                                                                                                                                        tGACCTTGCCGAAGCCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTCGACGTACCggg       >  2:21042/1‑147 (MQ=255)
                                                                                                                                         gACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTGGACGTACCGGGAAg    >  3:348424/1‑149 (MQ=255)
                                                                                                                                            cTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACggg                     >  4:72286/1‑129 (MQ=255)
                                                                                                                                            cTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACggg                     <  3:72286/129‑1 (MQ=255)
                                                                                                                                            cTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTGGACGTACCGGGAAGTg  >  7:369345/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GCAACTTTCGTTATAGCGTCACCGCGACTCCTGCCGGACATCATGATGAAAGCAAAGCCTGGCTACATTTTGGCAAATACGATCGCTATGACGACAAATACACCTATCCGGCAATGATGAACGGTTACATCCAGTATGACCTTGCCGAAGGCATCACCTGGATGAACGGTCTGGAAATCACCGACGGCACAGGACAGCTCTATCTCACGGGCCTGCTAACTCCTAACTTTGCCGCTCGCGCCTGGCACCATACCGGACGCGCCGACGGGCTGGACGTACCGGGAAGTG  >  NC_000913/3228336‑3228623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: