Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,836 G→A 8.5% 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 8.5% 218.9 / 6.4 70pseudogene (33/198 nt)ypaApseudogene, C‑terminal fragment, Transposase_31 family
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base A (3/3);  total (39/31)
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

TATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTG  >  NC_000913/2357689‑2357983
                                                                                                                                                   |                                                                                                                                                   
tATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTg                                                                                                                                                     >  3:67949/1‑148 (MQ=255)
   gCATGAACAAGTCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCgg                                                                                                                                                  >  6:170801/1‑148 (MQ=255)
    cATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGc                                                                                                                                                 <  4:67949/148‑1 (MQ=255)
       gAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCag                                                                                                                                               >  7:403192/1‑147 (MQ=255)
        aaCAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagg                                                                                                                                              >  7:15478/1‑147 (MQ=255)
                   aaaaTTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaa                                                                                                                                   <  6:337852/147‑1 (MQ=255)
                      aTTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtg                                                                                                                               <  7:169678/148‑1 (MQ=255)
                             tGCGCATGCTGGAGCAGGGCTCTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGATAGATTAATGACTATTGCTGACCGGCTGCGGCAGGAAGGACATCAAATtggctgg                                                                                                                          <  7:418633/146‑1 (MQ=255)
                                 cATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggct                                                                                                                            <  6:303957/140‑1 (MQ=255)
                                 cATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggct                                                                                                                            >  5:303957/1‑140 (MQ=255)
                                    gCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAgg                                                                                                                  >  5:377739/1‑147 (MQ=255)
                                       ggAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaaa                                                                                                              >  7:125140/1‑148 (MQ=255)
                                        gAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTaa                                                                                                               <  5:170801/146‑1 (MQ=255)
                                            aGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGa                                                                                                         <  8:403192/148‑1 (MQ=255)
                                                  ttGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGttt                                                                                                   >  6:131268/1‑148 (MQ=255)
                                                    gAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                                  <  8:125140/147‑1 (MQ=255)
                                                     aaCGTGAGATTGTGCTGGCGACAACCCAGCTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCa                                                                                                <  6:253730/148‑1 (MQ=255)
                                                          gagaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtg                                                                                                                               <  4:377527/112‑1 (MQ=255)
                                                          gagaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtg                                                                                                                               >  3:377527/1‑112 (MQ=255)
                                                           agaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAAc                                                                                           <  6:377739/147‑1 (MQ=255)
                                                                   cTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAt                                                                                                                                      <  1:107705/96‑1 (MQ=255)
                                                                   cTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAt                                                                                                                                      >  2:107705/1‑96 (MQ=255)
                                                                       cGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAt                                                                              >  7:419056/1‑148 (MQ=255)
                                                                            aCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagga                                                                                                                                             >  1:399534/1‑80 (MQ=255)
                                                                            aCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCagga                                                                                                                                             <  2:399534/80‑1 (MQ=255)
                                                                               cAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgc                                                                      >  8:331237/1‑148 (MQ=255)
                                                                                      cTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTTATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGATGGCAGGAAGGTAAATTAGAAGCTTTGCATGCACACGCCATTAAAATTGCTGTGCGCATGCt                                                               >  2:296657/1‑148 (MQ=255)
                                                                                      cTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCt                                                               >  2:155583/1‑148 (MQ=255)
                                                                                           gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           >  4:156329/1‑147 (MQ=255)
                                                                                           gCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           >  2:103521/1‑147 (MQ=255)
                                                                                              gATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                          >  8:189244/1‑145 (MQ=255)
                                                                                               atatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACa                                                         <  7:189244/145‑1 (MQ=255)
                                                                                                tatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAAc                                                                                           <  7:331165/110‑1 (MQ=255)
                                                                                                tatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAAc                                                                                           >  8:331165/1‑110 (MQ=255)
                                                                                                tatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                     >  1:182660/1‑132 (MQ=255)
                                                                                                tatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                     <  2:182660/132‑1 (MQ=255)
                                                                                                       aaCTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATc                                              >  7:181280/1‑148 (MQ=255)
                                                                                                           gTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGa                                          >  4:98870/1‑148 (MQ=255)
                                                                                                              aTTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAcc                                        <  7:110504/147‑1 (MQ=255)
                                                                                                              aTTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                       <  1:103521/148‑1 (MQ=255)
                                                                                                               ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGTCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGATTTGATCGTTACCAg                                      >  5:456/1‑148 (MQ=255)
                                                                                                               ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAg                                      <  8:419056/148‑1 (MQ=255)
                                                                                                               ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAg                                      >  1:216553/1‑148 (MQ=255)
                                                                                                               ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGAATGAACAAGCCAATAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAc                                         >  8:110504/1‑145 (MQ=255)
                                                                                                                tAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAgg                                     >  4:112430/1‑148 (MQ=255)
                                                                                                                   aaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                     >  6:139730/1‑113 (MQ=255)
                                                                                                                   aaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                     <  5:139730/113‑1 (MQ=255)
                                                                                                                       agagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                                  <  5:47100/80‑1 (MQ=255)
                                                                                                                       agagaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTg                                                                                                  >  6:47100/1‑80 (MQ=255)
                                                                                                                           aTTAATGACTATTGCTGAACGGCAGCGGCAGCTAGGACATCAAATTGGCTGGAAGGAATGTAAATTAGAATGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGGTGGAACAGGGCTTTGATCGTGACCAGCTGCTCGCGGcc                          <  6:456/148‑1 (MQ=37)
                                                                                                                               aTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACa                                                                                          >  4:154821/1‑80 (MQ=255)
                                                                                                                               aTGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACa                                                                                          <  3:154821/80‑1 (MQ=255)
                                                                                                                                tGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAcc                       >  2:198867/1‑146 (MQ=255)
                                                                                                                                tGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAcc                       >  1:324159/1‑146 (MQ=255)
                                                                                                                                 gACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAg                    >  7:84505/1‑148 (MQ=255)
                                                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCtt                                                                         >  4:196829/1‑91 (MQ=255)
                                                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCtt                                                                         <  3:196829/91‑1 (MQ=255)
                                                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTcgc                              >  2:315015/1‑134 (MQ=255)
                                                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTcgc                              <  1:315015/134‑1 (MQ=255)
                                                                                                                                         cTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCg             <  3:98870/147‑1 (MQ=255)
                                                                                                                                          tGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                                    >  8:68115/1‑59 (MQ=14)
                                                                                                                                          tGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                                    <  7:68115/59‑1 (MQ=14)
                                                                                                                                           gAACGGCTGCGGCATGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg          >  4:76664/1‑148 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt                                                   <  1:47180/106‑1 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt                                                   >  2:47180/1‑106 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTg                                           <  1:361138/114‑1 (MQ=255)
                                                                                                                                            aaCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTg                                           >  2:361138/1‑114 (MQ=255)
                                                                                                                                            aaCGGCTACGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTaaa                                                                                <  7:410822/77‑1 (MQ=255)
                                                                                                                                             aCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCATCTAAGCGAAGcc        >  5:52833/1‑148 (MQ=255)
                                                                                                                                             aCGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGcc        >  5:125439/1‑148 (MQ=255)
                                                                                                                                                 cTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATc    <  8:181280/148‑1 (MQ=255)
                                                                                                                                                   gCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTg  >  3:67322/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTG  >  NC_000913/2357689‑2357983

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