Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,653 A→G 5.6% L250L (CTA→CTG yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6530AG 5.6% 270.0 / 2.1 88L250L (CTA→CTGyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (36/46);  new base G (2/3);  total (39/49)
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

GCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACT  >  NC_000913/2357508‑2357795
                                                                                                                                                 |                                                                                                                                              
gCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACgg                                                                                                                                             >  8:216390/1‑149 (MQ=255)
 cTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGc                                                                                                                                            <  7:236368/149‑1 (MQ=255)
   tCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAg                                                                                                                                          >  8:282236/1‑149 (MQ=255)
      cTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAgg                                                                                                                                         >  8:285586/1‑147 (MQ=255)
       tAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAg                                                                                                                                      >  8:116779/1‑149 (MQ=255)
            aCAGCCCTCTACACACACTGTTAAATTATCTGCTGCCATGCGGCGATACCTCCCTTTTCACCCGTGTTATTGAGGAGCTTGCCGAACGTTCACCACTACAAAAGGAGATCTTAATGACTATTGCTGCACGGCTACGGCAGGAAGGGCag                                                                                                                                 <  5:242506/149‑2 (MQ=255)
               gccagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCa                                                                                                                               >  2:142237/1‑148 (MQ=255)
                 cagcTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAAt                                                                                                                            >  5:228081/1‑149 (MQ=255)
                    cTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtgg                                                                                                                         <  5:251789/149‑1 (MQ=255)
                     tACAAACACTGTTTAATTTTCTGCTGCAATGCGGCGATCCCTCACGTTTGACCCGTTTTATTGAGGAGTTTGCAGAACGTTCCCCACTACAAAAGGAGAGATTAATTACTATTTCTGAACGGATACTGCAGGTAGGGCTTCAAATtgga                                                                                                                        >  5:108313/1‑148 (MQ=255)
                       cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCTCTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggct                                                                                                                       <  3:182754/148‑1 (MQ=255)
                        aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                     >  1:271886/1‑149 (MQ=255)
                          acacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                      >  6:175749/1‑146 (MQ=255)
                          acacTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctg                                                                                                                      <  5:175749/146‑1 (MQ=255)
                           cacTGTTTAATTATCTGCTGCAATGCGGAGATACCTCCCGTTTCACCCTTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                  >  2:231240/1‑149 (MQ=255)
                                tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACTTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTTAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGTAAg                                                                                                              <  6:333174/148‑1 (MQ=255)
                                tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                               <  3:48093/147‑1 (MQ=255)
                                tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                               <  1:123681/147‑1 (MQ=255)
                                tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                             >  7:82325/1‑149 (MQ=255)
                                tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAg                                                                                                              <  7:216390/148‑1 (MQ=255)
                                  tAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                            <  8:283823/148‑1 (MQ=255)
                                  aaaTTATCTGCTGCAATGCCGCGATACCTCCCGTTTCACCCGTTTTATTGACGAGAGTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                            <  6:108313/147‑1 (MQ=255)
                                     ttATCTGCTGCTATCCGGCGATACCTACCTTTTCCACCGCTTAATCGAGGATACTGCCGAACGTTCACCCCTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGCCATCAAATTGGCTGGCAGGAAGGTaaat                                                                                                        <  2:77493/149‑4 (MQ=17)
                                       aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                      >  6:183216/1‑149 (MQ=255)
                                        tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACTTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                               <  6:246328/139‑1 (MQ=255)
                                        tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                               >  5:246328/1‑139 (MQ=255)
                                        tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                      <  1:203503/148‑1 (MQ=255)
                                           gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACTGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                   >  1:236741/1‑148 (MQ=255)
                                           gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                   >  4:326496/1‑148 (MQ=255)
                                            ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                   <  1:142237/147‑1 (MQ=255)
                                                   gcggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                           >  4:100351/1‑148 (MQ=255)
                                                   gcggcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCAt                                                                                           >  4:110190/1‑148 (MQ=255)
                                                      gcgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaa                                                                                       <  7:116779/149‑1 (MQ=255)
                                                       cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                      <  6:228081/149‑1 (MQ=255)
                                                       cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                      <  2:271886/149‑1 (MQ=255)
                                                       cgATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                      <  3:110190/149‑1 (MQ=255)
                                                             cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCtt                                                                                >  1:312357/1‑149 (MQ=255)
                                                             cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCtt                                                                                >  3:1679/1‑149 (MQ=255)
                                                             cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATGTATGAGTATTGCTGATCGGCTCCGGCAGGAAGGGCAGCAAATTGGATGGCAGGAAGGTATGCATGAACAAGCCATTTAAATTGCtt                                                                                >  7:263326/1‑149 (MQ=255)
                                                              tCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                               >  1:118330/1‑149 (MQ=255)
                                                               cccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTg                                                                              >  6:205447/1‑149 (MQ=255)
                                                                ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                               <  8:82325/147‑1 (MQ=255)
                                                                  gcttCACCCGTGTTATTTAGGAGATGGCCTCACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                           <  2:351095/147‑1 (MQ=25)
                                                                   tttCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCa                                                                          <  4:1679/149‑1 (MQ=255)
                                                                        cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                     >  8:230610/1‑149 (MQ=255)
                                                                        cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                     <  3:44282/149‑1 (MQ=255)
                                                                          cGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGACCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                     >  5:95466/1‑147 (MQ=255)
                                                                              ttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCa                                                                <  3:127412/148‑1 (MQ=255)
                                                                                aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                 <  5:217824/129‑1 (MQ=255)
                                                                                aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                                 >  6:217824/1‑129 (MQ=255)
                                                                                   gaggagATTGCCGAACGTTCAACACTACAAAAGGAGAGATTAATGACTATTTCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATAAAAATTGCTTTGCGAATGCTGGAGCAGGGCt                                                           <  8:263326/148‑1 (MQ=255)
                                                                                      gagaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGa                                                       <  4:257060/149‑1 (MQ=255)
                                                                                        gaTTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                     <  7:230610/149‑1 (MQ=255)
                                                                                          ttGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGGGCAGGGCTTTGAACGt                                                   <  7:285586/149‑1 (MQ=255)
                                                                                            gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTg                                                  <  2:118330/148‑1 (MQ=255)
                                                                                             ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGt                                                   >  5:66675/1‑146 (MQ=255)
                                                                                              cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAg                                                                                                                                      >  6:200457/1‑62 (MQ=255)
                                                                                              cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAg                                                                                                                                      <  5:200457/62‑1 (MQ=255)
                                                                                              cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTgaga                                               <  1:231240/149‑1 (MQ=255)
                                                                                                 aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                             >  8:119480/1‑84 (MQ=255)
                                                                                                 aCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                             <  7:119480/84‑1 (MQ=255)
                                                                                                  cGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAgg                                                                                                                                     >  6:316085/1‑59 (MQ=255)
                                                                                                  cGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAgg                                                                                                                                     <  5:316085/59‑1 (MQ=255)
                                                                                                      caccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCt                                        <  8:15574/148‑1 (MQ=255)
                                                                                                       accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTg                                       <  4:101527/148‑1 (MQ=255)
                                                                                                         cacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATTAATGAAGAATACAATAAAATAGCTTTTACCATGCTGTAGCACGTCATTGAAAGTGAGATtgatagggc                                     >  7:300280/1‑141 (MQ=17)
                                                                                                          acTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGa                                   >  6:115385/1‑149 (MQ=255)
                                                                                                           cTACATAAGGAGAGATTAATGACTATTGCTGTACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCTGGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAGGTGAGATTGTGCTGGCGAc                                  <  2:236741/149‑1 (MQ=25)
                                                                                                              cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGc                                                                                                                                   <  7:128478/49‑1 (MQ=255)
                                                                                                              cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGc                                                                                                                                   >  8:128478/1‑49 (MQ=255)
                                                                                                              cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGaa                                                                                                               >  2:242462/1‑69 (MQ=14)
                                                                                                               aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCtt                                                          <  8:313704/121‑1 (MQ=255)
                                                                                                               aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCtt                                                          >  7:313704/1‑121 (MQ=255)
                                                                                                               aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                              >  2:201215/1‑149 (MQ=255)
                                                                                                                  aGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCaa                           <  5:183216/149‑1 (MQ=255)
                                                                                                                   ggAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTAGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAAc                          <  4:338366/149‑1 (MQ=255)
                                                                                                                          ttAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGa                   <  7:282236/149‑1 (MQ=255)
                                                                                                                            aaTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGAt                  <  3:326496/148‑1 (MQ=255)
                                                                                                                             aTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                     <  8:251574/128‑1 (MQ=255)
                                                                                                                             aTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                     >  7:251574/1‑128 (MQ=255)
                                                                                                                               gACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                     <  2:42631/110‑1 (MQ=255)
                                                                                                                               gACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                     >  1:42631/1‑110 (MQ=255)
                                                                                                                                 cTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGat            <  1:201215/149‑1 (MQ=255)
                                                                                                                                  tATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGata           >  5:240887/1‑149 (MQ=255)
                                                                                                                                    ttGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATAtt         >  7:32945/1‑149 (MQ=255)
                                                                                                                                    ttGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCGGGCGACAACCCAACTCACTGATGCTGatat          <  3:100351/148‑1 (MQ=255)
                                                                                                                                    ttGCGGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGata           <  5:205447/147‑1 (MQ=255)
                                                                                                                                           aCGGATACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAATTTGCTTTGCGCATGCTGTAGCAGGGCTTTGTACGTGAGATTGTGCTGGCTACAACCCAACTCACTGATGCTGATATTCCGAACt  <  6:66675/149‑1 (MQ=38)
                                                                                                                                                 |                                                                                                                                              
GCTTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACT  >  NC_000913/2357508‑2357795

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: