Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,180,168 T→A 6.4% L25H (CTT→CAT)  yegU → ADP‑ribosylglycohydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,180,1680TA 6.4% 202.1 / 3.7 63L25H (CTT→CAT) yegUADP‑ribosylglycohydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base A (2/2);  total (35/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 = 8.97e-01

CGCGAATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGATGCGT  >  NC_000913/2180021‑2180311
                                                                                                                                                   |                                                                                                                                               
cgcgAATCCGACAACGAAATTACGGCTATCAAGGTCGAGGATCGCGATATTGCGTTGCCACAAGGGGAAGTTAAATGAAAACAGACCGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCtt                                                                                                                                                >  1:248484/1‑149 (MQ=255)
 gcgAATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCttt                                                                                                                                               <  8:10476/149‑1 (MQ=255)
 gcgAATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCttt                                                                                                                                               >  3:77430/1‑149 (MQ=255)
  cgAATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTg                                                                                                                                              <  5:182082/149‑1 (MQ=255)
    aaTCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAATCAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGc                                                                                                                                            >  8:138370/1‑149 (MQ=255)
     aTCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGc                                                                                                                                            >  2:275784/1‑148 (MQ=255)
               gAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCcgc                                                                                                                                 >  3:175852/1‑149 (MQ=255)
                      cGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAg                                                                                                                          <  1:275784/149‑1 (MQ=255)
                      cGGCTATCAAGGTCGATGATCGCGATATTGCGGTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAg                                                                                                                          <  2:50807/149‑1 (MQ=255)
                          tATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGcaca                                                                                                                      >  2:270036/1‑149 (MQ=255)
                            tCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGTAGTTAAATGAAAACAGGACGTATTAGCGGTGCTCTTGATGGGCAGGGGTTAGGGGATGCGATGGGGATGACCTCCGAGCGGTGGCCACGCGGCCGCGTTAAAGCACACt                                                                                                                    >  4:65298/1‑149 (MQ=255)
                            tCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACt                                                                                                                    <  4:175852/149‑1 (MQ=255)
                               aGGTCGATGATCGCGATATTGCGTTGAAACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTg                                                                                                                 >  6:176538/1‑149 (MQ=255)
                                     atgatCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGa                                                                                                           >  4:45435/1‑149 (MQ=255)
                                        atCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGAtt                                                                                                         <  4:77430/148‑1 (MQ=255)
                                         tCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGa                                                                                                       >  1:210464/1‑149 (MQ=255)
                                              atatTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGTGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGtt                                                                                                  <  3:257461/149‑1 (MQ=255)
                                               tatTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGttt                                                                                                 >  5:221148/1‑149 (MQ=255)
                                                 ttGCGTTGACACAAGGGGAAGTTCACTGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGtttt                                                                                                <  5:242818/148‑1 (MQ=255)
                                                    cGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTttcttc                                                                                            >  3:26998/1‑149 (MQ=255)
                                                           acaAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGAc                                                                                      >  6:373784/1‑148 (MQ=255)
                                                             aaGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTGTGGCTGGATTGACCGTTTTCTTCCTGGACCa                                                                                    >  1:20628/1‑148 (MQ=255)
                                                              aGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCaaa                                                                                  <  3:45435/149‑1 (MQ=255)
                                                              aGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCaaa                                                                                  <  2:210464/149‑1 (MQ=255)
                                                                 ggAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAg                                                                                 >  6:362465/1‑147 (MQ=255)
                                                                 ggAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGga                                                                               <  7:138370/149‑1 (MQ=255)
                                                                  gAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCt                                                                                          >  6:250273/1‑137 (MQ=255)
                                                                  gAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCt                                                                                          <  5:250273/137‑1 (MQ=255)
                                                                    aGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGaa                                                                            <  4:26998/149‑1 (MQ=255)
                                                                      ttAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGCCGGTTTTCTTCCTGGACCAAAGGAGAAt                                                                           >  7:296667/1‑148 (MQ=255)
                                                                      ttAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCATTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGgag                                                                              <  1:221420/145‑1 (MQ=255)
                                                                      ttAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCATTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGgag                                                                              >  2:221420/1‑145 (MQ=255)
                                                                         aaTGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCaaa                                                                                  <  7:46618/138‑1 (MQ=255)
                                                                         aaTGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCaaa                                                                                  >  8:46618/1‑138 (MQ=255)
                                                                         aaTGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAAc                                                                        <  3:65298/148‑1 (MQ=255)
                                                                          aTGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAAcgc                                                                      <  5:373784/149‑1 (MQ=255)
                                                                             aaaaCAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGc                                                                   >  5:165889/1‑149 (MQ=255)
                                                                                     aCGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTAtt                                                           <  5:362465/149‑1 (MQ=255)
                                                                                         aTTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGgaga                                                                             >  6:7819/1‑127 (MQ=255)
                                                                                         aTTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGgaga                                                                             <  5:7819/127‑1 (MQ=255)
                                                                                                      tttATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTAtt                                                           <  4:164662/132‑1 (MQ=255)
                                                                                                      tttATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTAtt                                                           >  3:164662/1‑132 (MQ=255)
                                                                                                        tATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCAc                                        >  6:70936/1‑149 (MQ=255)
                                                                                                          tGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCAcc                                       >  3:240969/1‑148 (MQ=255)
                                                                                                            ggCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCATTGGCCACGCAGCCGCGTTAAAGCACACTTTTGCTGGATTGACCGTTTTCTTCCTg                                                                                         >  1:364379/1‑96 (MQ=255)
                                                                                                            ggCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCATTGGCCACGCAGCCGCGTTAAAGCACACTTTTGCTGGATTGACCGTTTTCTTCCTg                                                                                         <  2:364379/96‑1 (MQ=255)
                                                                                                              cAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACcgacga                                  <  4:240969/149‑1 (MQ=255)
                                                                                                                   gTTAGGGGATGCGATGGGGATGCCCTCCGGGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCt                             >  4:29016/1‑149 (MQ=255)
                                                                                                                        gggATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCgcgc                                                 >  8:140568/1‑124 (MQ=255)
                                                                                                                        gggATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCgcgc                                                 <  7:140568/124‑1 (MQ=255)
                                                                                                                              cGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCt                   <  1:270036/148‑1 (MQ=255)
                                                                                                                               gctGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCtgt                 <  3:283250/147‑1 (MQ=255)
                                                                                                                               gATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCtgt                 >  8:374580/1‑149 (MQ=255)
                                                                                                                                    ggATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACcgac                                    >  3:8043/1‑125 (MQ=255)
                                                                                                                                    ggATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACcgac                                    <  4:8043/125‑1 (MQ=255)
                                                                                                                                    ggATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTg            <  2:248484/149‑1 (MQ=255)
                                                                                                                                     gATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTgg           >  6:218964/1‑149 (MQ=255)
                                                                                                                                     gATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTgg           >  7:344139/1‑149 (MQ=255)
                                                                                                                                          ccTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAAcg                                                                       <  6:214563/84‑1 (MQ=255)
                                                                                                                                          ccTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAAcg                                                                       >  5:214563/1‑84 (MQ=255)
                                                                                                                                          ccTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGAt      >  8:131616/1‑149 (MQ=255)
                                                                                                                                           cTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGATg     >  3:260829/1‑149 (MQ=255)
                                                                                                                                            tCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGCCCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGATGc    >  5:257016/1‑149 (MQ=255)
                                                                                                                                            tCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCt                   >  1:47986/1‑134 (MQ=255)
                                                                                                                                            tCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCt                   <  2:47986/134‑1 (MQ=255)
                                                                                                                                               gAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGATGCGt  >  7:82857/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGCGAATCCGACAACGAAATTACGGCTATCAAGGTCGATGATCGCGATATTGCGTTGACACAAGGGGAAGTTAAATGAAAACAGAACGTATTCTCGGTGCTCTTTATGGGCAGGCGTTAGGGGATGCGATGGGGATGCCCTCCGAGCTTTGGCCACGCAGCCGCGTTAAAGCACACTTTGGCTGGATTGACCGTTTTCTTCCTGGACCAAAGGAGAATAACGCGGCCTGTTATTTTAACCGCGCCGAATTCACCGACGATACCTCGATGGCGCTGTGTCTGGCGGATGCGT  >  NC_000913/2180021‑2180311

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