Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,689 T→C 7.8% G262G (GGT→GGC yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6890TC 7.8% 189.9 / 4.5 67G262G (GGT→GGCyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (25/35);  new base C (2/3);  total (29/38)
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

AATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATT  >  NC_000913/2357543‑2357824
                                                                                                                                                  |                                                                                                                                       
aaTTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAt                                                                                                                                       >  1:18379/1‑149 (MQ=255)
  ttATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATg                                                                                                                                      >  8:352625/1‑148 (MQ=255)
     tCTGCTGGAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                                                   <  8:320799/148‑1 (MQ=255)
     tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                                                   >  5:236679/1‑148 (MQ=255)
         ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGTGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCAGCAAATTGGCTGGCAGGAAGGTATGCATGa                                                                                                                                 <  5:61795/146‑1 (MQ=255)
           gcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAg                                                                                                                            >  3:160777/1‑149 (MQ=255)
            cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAGCAAGAACAAga                                                                                                                           >  5:249400/1‑148 (MQ=255)
             aaTGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGATATTGACGAACGTTAACCACTACAAAAGGAGAGATTAAAGACTATTGCTGAACGGATACGGCAGGAAGGGCATCAAATTCGCTGCCAGGAAGGAATGCATGAACACGc                                                                                                                           >  3:118006/1‑148 (MQ=255)
                         ccTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                                              <  4:160777/149‑1 (MQ=255)
                          cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCtt                                                                                                             >  2:82357/1‑149 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            >  1:421439/1‑149 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            >  1:415578/1‑149 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            >  8:229912/1‑149 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            >  2:334235/1‑149 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                                              <  2:218016/147‑1 (MQ=255)
                           tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                                              >  1:218016/1‑147 (MQ=255)
                             ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACTGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            <  1:257557/147‑1 (MQ=255)
                             ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            <  7:296237/147‑1 (MQ=255)
                             ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                                                          <  6:236679/149‑1 (MQ=255)
                                     caagTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGt                                                                                                                                         <  6:289031/107‑1 (MQ=255)
                                      ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTTGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                                                  >  4:187110/1‑148 (MQ=255)
                                      ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                                                  >  8:53637/1‑148 (MQ=255)
                                             aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                                                   <  8:43922/108‑1 (MQ=255)
                                             aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                                                   >  7:43922/1‑108 (MQ=255)
                                                gaggagATTGCCGAACGTTCACCACTACAAAAGTAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGACATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCt                                                                                        <  1:396523/148‑1 (MQ=255)
                                                gaggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCt                                                                                        <  3:326278/148‑1 (MQ=255)
                                                 aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                                                      <  3:117954/149‑1 (MQ=255)
                                                          ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTTCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGGTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              <  2:47872/148‑1 (MQ=255)
                                                          ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              <  5:331852/148‑1 (MQ=255)
                                                            gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATACGCTGGCTGGAAGGAATGAATGAACAAGCCCTTAAAAATGCTTTGCGGATGCTGGAGCAGGGCGTTGAACGTGAGAt                                                                           >  5:161498/1‑149 (MQ=25)
                                                               cGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCGTTGAACGTGAGATtgt                                                                        <  5:191183/149‑1 (MQ=255)
                                                                 ttCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGc                                                                      <  1:82357/149‑1 (MQ=255)
                                                                    accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGCATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              >  7:146755/1‑138 (MQ=255)
                                                                    accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGCATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              <  8:146755/138‑1 (MQ=255)
                                                                     ccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                                  <  7:352625/149‑1 (MQ=255)
                                                                      cacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                                                  >  8:1550/1‑148 (MQ=255)
                                                                         tACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                                                        >  5:191526/1‑107 (MQ=255)
                                                                         tACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                                                        <  6:191526/107‑1 (MQ=255)
                                                                             aaaTGAGAGATTAATGACTATTGCTGAACGGCTACGGCATGAAGGGCATCAAATTGGCGGGCAGGAAGGTATGCATGAACAAGCCTTTAAAATTGCTTTGCGCTTGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAccc                                                          <  4:118006/149‑1 (MQ=37)
                                                                              aaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                                  <  5:268659/124‑1 (MQ=255)
                                                                              aaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAtt                                                                          >  8:120754/1‑132 (MQ=255)
                                                                              aaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGAtt                                                                          <  7:120754/132‑1 (MQ=255)
                                                                              aaGGAGAGATAAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                                                  >  6:268659/1‑124 (MQ=255)
                                                                               aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                                         >  8:144532/1‑68 (MQ=40)
                                                                               aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGCCATCAAATTGGCTGGCAGGAAGGt                                                                                                                                         <  7:144532/68‑1 (MQ=255)
                                                                                    agaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCAc                                                   >  2:84962/1‑149 (MQ=255)
                                                                                     gaTTAATGACTATTGCTGAACGTCTACGGCAGGAAGGGCATCAACTTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGACCGTGAGATTGTGCTGGCGACAACCCAACTCACt                                                  <  5:336822/149‑1 (MQ=255)
                                                                                     gaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGCATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACt                                                  >  4:123490/1‑149 (MQ=255)
                                                                                           tGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTaaa                                                                                                                                      >  7:88173/1‑57 (MQ=255)
                                                                                           tGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                                   <  4:211914/126‑1 (MQ=255)
                                                                                           tGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                                                   >  3:211914/1‑126 (MQ=255)
                                                                                               tATTGCTGAGCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            <  2:134132/81‑1 (MQ=255)
                                                                                               tATTGCTGAGCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                                                            >  1:134132/1‑81 (MQ=255)
                                                                                               tATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGata                                        >  2:273345/1‑149 (MQ=255)
                                                                                                    cTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTcc                                    <  1:273345/148‑1 (MQ=255)
                                                                                                       aaCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAAc                                >  8:117093/1‑149 (MQ=255)
                                                                                                             tACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAt                          <  7:229912/149‑1 (MQ=255)
                                                                                                               cGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCAt                          <  7:117093/147‑1 (MQ=255)
                                                                                                                   aGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAg                    <  7:53637/149‑1 (MQ=255)
                                                                                                                     gAAGGGCATCAAATTGGCTTGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGga                  <  3:187110/149‑1 (MQ=255)
                                                                                                                     gAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGga                  <  1:84962/149‑1 (MQ=255)
                                                                                                                      aaGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgag                 <  1:334235/149‑1 (MQ=255)
                                                                                                                      aaGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgag                 <  7:1550/149‑1 (MQ=255)
                                                                                                                       aGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGgaga                <  2:18379/149‑1 (MQ=255)
                                                                                                                               aaaTTGGCTGGCAGGAAGGCATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAAt         <  3:123490/148‑1 (MQ=255)
                                                                                                                                    ggctggcAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGATATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTa    <  2:421439/148‑1 (MQ=255)
                                                                                                                                    ggctggcAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACGCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTa    <  2:415578/148‑1 (MQ=255)
                                                                                                                                     cccgccAGGAACGCATGCATGAACAAGCCATTAATTTTTCTTTGCTCTTGCTCGAGCTCGCTCTTCAAAGTCTTATTGTGTTTGCGACAACACTACTCACTTAAGCTGTTCTTCGTAAATGACATTAAAAGGAGAGATTAATGACTAtt  <  6:161498/144‑1 (MQ=255)
                                                                                                                                       tggcAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              <  4:146962/71‑1 (MQ=255)
                                                                                                                                       tggcAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                                                              >  3:146962/1‑71 (MQ=255)
                                                                                                                                                  |                                                                                                                                       
AATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATT  >  NC_000913/2357543‑2357824

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: