Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,377,597 G→A 5.0% A256T (GCA→ACA)  ycjS → putative NADH‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,377,5970GA 5.0% 258.6 / 2.1 79A256T (GCA→ACA) ycjSputative NADH‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (31/44);  new base A (2/2);  total (33/46)
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.96e-01

AATGCGCATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCAACGGGAAA  >  NC_000913/1377450‑1377726
                                                                                                                                                   |                                                                                                                                 
aaTGCGCATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGAAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGc                                                                                                                                  <  4:236191/149‑1 (MQ=255)
 aTGCGCATAGCTTTCAAAAGATCGTCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTAGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCa                                                                                                                                 <  8:241325/149‑1 (MQ=255)
  tGCGCATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCAc                                                                                                                                >  3:141587/1‑149 (MQ=255)
     gcATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTACACTc                                                                                                                              >  1:81609/1‑148 (MQ=255)
     gcATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTACACTc                                                                                                                              >  1:76497/1‑148 (MQ=255)
      cATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCaa                                                                                                                            <  3:328581/149‑1 (MQ=255)
       aTAGCTTTCAAAAGATCGGCAAGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTGGCACTCaa                                                                                                                            >  3:165307/1‑148 (MQ=255)
          gCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACAt                                                                                                                         <  1:199089/148‑1 (MQ=255)
            tttCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCcg                                                                                                                      >  5:15671/1‑149 (MQ=255)
                aaaaGATCTGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGaa                                                                                                                  <  7:300634/149‑1 (MQ=255)
                aaaaGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGaa                                                                                                                  >  2:353359/1‑149 (MQ=255)
                 aaaGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCAGTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGaa                                                                                                                  >  7:299344/1‑148 (MQ=255)
                    gATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGt                                                                                                              <  7:367985/149‑1 (MQ=255)
                     aTCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGCGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTc                                                                                                             <  2:287600/149‑1 (MQ=255)
                      tCGGCACGCACAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCg                                                                                                            <  3:178341/149‑1 (MQ=255)
                        ggCACGCATAAGAGAGGTGGTAAATTTGGTGAGTGTGATACGGCAACTTACAGCGTAGAAGATTCGCTGTTTGGCACCATTGATTTTAATAACGGCGGCATTCTGTGGCTGGAAACGGCATTTGCACTCAACATCCGCGAACAGGCGa                                                                                                           <  6:246208/148‑1 (MQ=255)
                        ggCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGa                                                                                                           <  1:254792/148‑1 (MQ=255)
                        ggCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGa                                                                                                           <  3:186875/148‑1 (MQ=255)
                         gCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGAtt                                                                                                         <  8:97349/149‑1 (MQ=255)
                           aCGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACa                                                                                                                <  5:163807/140‑1 (MQ=255)
                           aCGCAAAAGAGATGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACa                                                                                                                >  6:163807/1‑140 (MQ=255)
                            cGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTAt                                                                                                       >  1:292012/1‑148 (MQ=255)
                            cGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTAt                                                                                                       >  1:287435/1‑148 (MQ=255)
                            cGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTACACTCAACATCCGCGAACAGTCGATTATg                                                                                                      <  2:76497/149‑1 (MQ=255)
                            cGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTACACTCAACATCCGCGAACAGTCGATTATg                                                                                                      <  2:81609/149‑1 (MQ=255)
                              cAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGaa                                                                                                    >  4:247678/1‑149 (MQ=255)
                                      cTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCt                                                                                            >  4:345127/1‑149 (MQ=255)
                                      cTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCt                                                                                            >  6:117147/1‑149 (MQ=255)
                                         tgGTCTATTTGGTGAGTGGGATGCGGCATCTTACAGCGTCGAAGTTTCGCTGTTTGTCTGCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCGGCGTACAGTGGATTATGAACGTCAGCtttt                                                                                         <  8:37217/149‑1 (MQ=255)
                                           gccAAGTTGGTGAGTTGGTCACGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCCTCCCGGCGGCATGCTGTGGCTGGAAACGGCATTTGCACTCCACATGCGCGAACCGTCGCTTATGACCGTCAGCTTTtgt                                                                                       <  4:283347/147‑1 (MQ=255)
                                              aaTTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTgtggt                                                                                    <  4:22533/149‑1 (MQ=255)
                                                        gTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTg                                                                          <  3:155873/149‑1 (MQ=255)
                                                           ggATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTg                                                                       >  3:198405/1‑149 (MQ=255)
                                                            gATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTgtgg                                                                                     <  4:136732/134‑1 (MQ=255)
                                                            gATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCGGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGc                                                                      <  6:227079/149‑1 (MQ=255)
                                                            gATCCGGCAAATTACAGCGTCGAAGATTCGCTGTTAGGCACCATTGAATTTCAGAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAAGAGTCGATTATGAACGTCAGCTTTTGTGGTGATATAGCtgctgc                                                                      >  3:223503/1‑149 (MQ=255)
                                                                cGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACg                                                                  >  1:265363/1‑149 (MQ=255)
                                                                     aCTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTg                                                                          <  1:216897/136‑1 (MQ=255)
                                                                     aCTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTg                                                                          >  2:216897/1‑136 (MQ=255)
                                                                     aCTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTGCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGtt                                                             <  5:342105/149‑1 (MQ=255)
                                                                      cTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGttt                                                            >  7:11923/1‑149 (MQ=255)
                                                                       ttACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGATACGTAATTTGCACTCAACATCAGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTc                                                           >  5:99026/1‑149 (MQ=255)
                                                                       ttACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTTGTGCGACGCTGTTTc                                                           >  6:28293/1‑149 (MQ=255)
                                                                         aCAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTcc                                                          <  5:246828/148‑1 (MQ=255)
                                                                         aCAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTc                                                           <  8:11923/147‑1 (MQ=255)
                                                                         aCAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTc                                                           <  8:150570/147‑1 (MQ=255)
                                                                          cAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCa                                                         <  3:319062/148‑1 (MQ=255)
                                                                          cAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAg                                                        <  1:89384/149‑1 (MQ=255)
                                                                          cAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAg                                                        <  2:292012/149‑1 (MQ=255)
                                                                            gCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAATGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGca                                                      <  3:247678/149‑1 (MQ=255)
                                                                               tCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACat                                                   <  6:15671/149‑1 (MQ=255)
                                                                                 gAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACatat                                                 >  5:309953/1‑149 (MQ=255)
                                                                                 gAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACatat                                                 >  7:108528/1‑149 (MQ=255)
                                                                                 gAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACatat                                                 >  8:206881/1‑149 (MQ=255)
                                                                                     aTTCGCTGTTTGTCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTa                                              >  2:215017/1‑148 (MQ=255)
                                                                                             tttGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGGCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATaa                                     <  2:287435/149‑1 (MQ=255)
                                                                                               tGGCACCATTGAATTTCATAACGGCGTCATGCTGTGGCTGGAAACGTCATTTGCACTCACCATCCGCGAACATGCGATTCTGACCGTCCGCTTTTGGTGTGATAAAGCTGGTGCGACGCGGTTTCCCGCACATATCTACACCGATaac                                    <  2:279689/148‑1 (MQ=255)
                                                                                                 gCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATaacaa                                  <  1:353359/148‑1 (MQ=255)
                                                                                                       ttGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGaa                           >  8:298302/1‑149 (MQ=255)
                                                                                                       ttGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGaa                           <  4:141587/149‑1 (MQ=255)
                                                                                                       ttGAATTTCAAACCGGCGGCATTCTGTGGCTTGAACAGGCCTATGCTCTCAACATCCGCGAACAGTCGATTATGAACGGCAGCTTTTGTGGTGATAAAGCTTGTTCGACGCTGTTTACAGCACATATCTACACCGATAACAACGGTGaa                           <  6:290001/149‑1 (MQ=255)
                                                                                                           aTTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTaa                       <  3:345127/149‑1 (MQ=255)
                                                                                                           aTTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTaa                       <  7:206881/149‑1 (MQ=255)
                                                                                                            tttCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGttt                                                            >  3:68747/1‑111 (MQ=255)
                                                                                                            tttCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGttt                                                            <  4:68747/111‑1 (MQ=255)
                                                                                                                 tAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGc                 >  6:72022/1‑149 (MQ=255)
                                                                                                                 tAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGc                 <  6:349135/149‑1 (MQ=255)
                                                                                                                    cggcggCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGa              <  7:298302/149‑1 (MQ=255)
                                                                                                                     ggcggcATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCTCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGAt             <  1:215017/149‑1 (MQ=255)
                                                                                                                      gcggcATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATg            <  6:70990/149‑1 (MQ=255)
                                                                                                                      gcggcATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATg            >  5:70990/1‑149 (MQ=255)
                                                                                                                         gcATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATa                                                                                <  5:38255/78‑1 (MQ=255)
                                                                                                                         gcATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATa                                                                                >  6:38255/1‑78 (MQ=255)
                                                                                                                           aTTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCAACg       >  5:2507/1‑149 (MQ=255)
                                                                                                                               tgtgGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCAACGGGa    >  2:233265/1‑148 (MQ=255)
                                                                                                                                 tgGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAATCTGGTTCGACGATGTTTCCAGCAGATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCAACGGGaaa  <  4:85171/148‑1 (MQ=255)
                                                                                                                                   gCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAAAGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCACCGGGaaa  >  1:138330/1‑146 (MQ=255)
                                                                                                                                     tGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGc                                                       <  8:137389/91‑1 (MQ=255)
                                                                                                                                     tGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGc                                                       >  7:137389/1‑91 (MQ=255)
                                                                                                                                                   |                                                                                                                                 
AATGCGCATAGCTTTCAAAAGATCGGCACGCAAAAGAGCTGTGGTCAATTTGGTGAGTGGGATCCGGCAACTTACAGCGTCGAAGATTCGCTGTTTGGCACCATTGAATTTCATAACGGCGGCATTCTGTGGCTGGAAACGTCATTTGCACTCAACATCCGCGAACAGTCGATTATGAACGTCAGCTTTTGTGGTGATAAAGCTGGTGCGACGCTGTTTCCAGCACATATCTACACCGATAACAACGGTGAATTAATGACGCTGATGCAACGGGAAA  >  NC_000913/1377450‑1377726

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