Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,836 G→A 7.0% pseudogene (33/198 nt) ypaA → pseudogene, C‑terminal fragment, Transposase_31 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,8360GA 7.0% 214.8 / 4.2 71pseudogene (33/198 nt)ypaApseudogene, C‑terminal fragment, Transposase_31 family
Reads supporting (aligned to +/- strand):  ref base G (38/28);  new base A (3/2);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

TATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATC  >  NC_000913/2357689‑2357981
                                                                                                                                                   |                                                                                                                                                 
tATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGc                                                                                                                                                  <  5:116053/149‑1 (MQ=255)
  tGCATGAACAAGCCATTTAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCgg                                                                                                                                                >  8:25351/1‑149 (MQ=255)
   gCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCgg                                                                                                                                                >  3:28960/1‑148 (MQ=255)
       gAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagg                                                                                                                                            >  8:253892/1‑148 (MQ=255)
        aaCAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagg                                                                                                                                            >  8:19137/1‑147 (MQ=255)
         aCAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCaggaag                                                                                                                                         >  7:133972/1‑149 (MQ=255)
            aGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGAc                                                                                                                                      <  6:119854/149‑1 (MQ=255)
                aTTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATc                                                                                                                                   >  6:49114/1‑148 (MQ=255)
                aTTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATc                                                                                                                                   >  8:155346/1‑148 (MQ=255)
                         gCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggctg                                                                                                                         >  8:138710/1‑149 (MQ=255)
                          cTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggctgg                                                                                                                        <  7:25351/149‑1 (MQ=255)
                                gcATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAgg                                                                                                                    >  8:301762/1‑147 (MQ=255)
                                  aTGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGCCAACCCAACTCAGGGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACGATTGCTGACCGGCTGCGGCAGGAAGGACATGAAATTGGCTGGCAGGAAg                                                                                                                 >  7:384700/1‑148 (MQ=255)
                                  aTGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                <  7:253892/149‑1 (MQ=255)
                                        gAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAAt                                                                                                           <  7:301762/148‑1 (MQ=255)
                                               gCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGt                                                                                                   <  7:155346/149‑1 (MQ=255)
                                                 tttGAACATGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGttt                                                                                                 >  6:304209/1‑149 (MQ=255)
                                                   tGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                                >  1:24238/1‑148 (MQ=255)
                                                    gAACATGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                               <  5:304209/148‑1 (MQ=255)
                                                     aaCGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCAt                                                                                             <  7:19137/149‑1 (MQ=255)
                                                      aCGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCAt                                                                                             <  5:49114/148‑1 (MQ=255)
                                                       cctGAGATAGTGCTGGCGACCACCCAACTCACTCAGGCACAAAAACCGAACTGGCATTAAAAGGAGAGAGAAATGCCTATTTCTGAACTGCTGCCACAGGAAGGACATCAAATTGCATCGCAGCAAGGAAAATAAGAAGGTTTGCAtc                                                                                            >  6:232841/3‑147 (MQ=255)
                                                       cGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGa                                                                                           <  4:28960/149‑1 (MQ=255)
                                                           agaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTATAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAGCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAAc                                                                                         >  3:137078/1‑147 (MQ=255)
                                                           agaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACaa                                                                                       <  7:138710/149‑1 (MQ=255)
                                                                 tgCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCAt                                                                                  <  2:24238/148‑1 (MQ=255)
                                                                      gCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTaaa                                                                              <  8:133972/147‑1 (MQ=255)
                                                                         aCAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAAtt                                                                                                          >  4:136034/1‑116 (MQ=255)
                                                                         aCAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAAtt                                                                                                          <  3:136034/116‑1 (MQ=255)
                                                                            aCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGAGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCttt                                                                      >  5:74099/1‑149 (MQ=255)
                                                                             cccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTg                                                                     >  6:124259/1‑149 (MQ=255)
                                                                              ccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGAGTATTGCTGAACGGCTGCTGCAGGAACCACATCAAATTGGCTGGAAGGAAGGTAAATTAGAAGGTTAGCATGAACAAGCCATTAAACTTGCtattc                                                                    >  2:90154/1‑145 (MQ=37)
                                                                              ccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgc                                                                    >  3:170364/1‑149 (MQ=255)
                                                                                       tGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                            >  1:364253/1‑148 (MQ=255)
                                                                                           gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                            <  7:246358/144‑1 (MQ=255)
                                                                                           gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                            >  8:246358/1‑144 (MQ=255)
                                                                                                  ttCCGAACTGTCATTATAAGGAGAGATTAATGAGTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt                                                 <  1:90154/148‑1 (MQ=255)
                                                                                                  ttCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCttt                                                <  4:170364/149‑1 (MQ=255)
                                                                                                     cGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAATGTAAATTAGAAGGTTTGCATGAACATTCCATTAAAATTGCTTTGCGCATGCTGGAACAGGCCTTTGAt                                             >  8:167719/1‑149 (MQ=255)
                                                                                                     cGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGAt                                             >  8:261707/1‑149 (MQ=255)
                                                                                                      gAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATc                                            <  2:364253/149‑1 (MQ=255)
                                                                                                         cTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGt                                          >  6:349783/1‑148 (MQ=255)
                                                                                                         cTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAGAAAAATGGCTGGCAGGAAGGTAAATTAGAAGGTTGGCATGAACATGCCATTAAAATTGCTTTGCGCGTGCTGGACCAGGGCTTTGATCGt                                          >  8:329816/1‑148 (MQ=25)
                                                                                                          tGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGa                                        >  4:56678/1‑149 (MQ=255)
                                                                                                            tCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAcc                                      >  8:360910/1‑149 (MQ=255)
                                                                                                              tatAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAAt                                                                                                                               <  1:344136/56‑1 (MQ=255)
                                                                                                              aTTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAg                                    >  8:376575/1‑149 (MQ=255)
                                                                                                                tAAAAGGAGAGATTAATGACTATTGCTGATCGGCTGCGGCTGGAAGGACTTGAAATTGGCTGGGAGGTAGGTATTTTAGAAGGTTTGGATGATGAAGCGATTAATATTGCTTTGCGGATGCTGGATGAGGGCTTTGATCGTGACCTGGt                                  >  5:84862/1‑149 (MQ=21)
                                                                                                                tAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGt                                  >  5:89899/1‑149 (MQ=255)
                                                                                                                caaaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                               <  1:161983/87‑1 (MQ=255)
                                                                                                                caaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGc                                                                                               >  2:161983/2‑88 (MQ=255)
                                                                                                                 aaaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTg                                 >  4:324446/1‑149 (MQ=255)
                                                                                                                   aaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCaggaag                                                                                                                                         >  8:383036/1‑43 (MQ=38)
                                                                                                                   aaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCaggaag                                                                                                                                         <  7:383036/43‑1 (MQ=38)
                                                                                                                   aaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCt                               >  3:3540/1‑149 (MQ=255)
                                                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGTTTAGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGCAACAGGGCTTTGATCGTGACCAGGTGCTc                              >  2:109090/1‑147 (MQ=255)
                                                                                                                       agagaTTAATGACTATTGCTGATCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGGTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTggcg                           <  1:109090/149‑1 (MQ=255)
                                                                                                                       agagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGATCGTGACCAGGTGCTcgcg                           >  4:308582/1‑149 (MQ=255)
                                                                                                                       agagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTcgcg                           <  6:89899/149‑1 (MQ=255)
                                                                                                                       acagaaaaaTGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGCTAAATAAGAAGGTTTGCATGAACAAGCCATAAAAATTGCATTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTcgcg                           <  6:84862/142‑1 (MQ=37)
                                                                                                                          gaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGcc                        <  7:376575/149‑1 (MQ=255)
                                                                                                                              aaTGACTATTGCTAAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGcccccc                    >  2:74637/1‑149 (MQ=255)
                                                                                                                                  aCTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGc                 <  4:3540/148‑1 (MQ=255)
                                                                                                                                    tATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCAt                                                                                             <  8:81534/70‑1 (MQ=255)
                                                                                                                                    tATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCAt                                                                                             >  7:81534/1‑70 (MQ=255)
                                                                                                                                           gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                     >  2:320543/1‑71 (MQ=255)
                                                                                                                                           gAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                     <  1:320543/71‑1 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg        <  5:124259/147‑1 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGcc      >  5:299099/1‑149 (MQ=255)
                                                                                                                                              cGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGa    <  3:324446/149‑1 (MQ=255)
                                                                                                                                                gCTGCGGCAGGAAGGACATGAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGATCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGGGACCAGGTGCTCGCGGCCGCCCAGCTACGCGAGGCCGCTc  >  6:163699/1‑149 (MQ=37)
                                                                                                                                                   |                                                                                                                                                 
TATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATC  >  NC_000913/2357689‑2357981

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: