Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,289,431 T→C G150G (GGT→GGC yraK → putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,4310TC99.7% 917.1 / ‑5.8 286G150G (GGT→GGCyraKputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (158/127);  total (158/128)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 291 total aligned reads displayed.

GACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289342‑3289520
                                                                                         |                                                                                         
gACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGAt                                                                                        <  1:3457505/93‑1 (MQ=255)
gACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGAt                                                                                        <  1:4183064/93‑1 (MQ=255)
 aCCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGCATATTCTGCTGCAACCCACACGGGTGTTTGCCAGGACGATa                                                                                       >  1:2291670/1‑93 (MQ=255)
   cAGGTTGTCTCTGTTAACAACCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACg                                                                                     >  1:1324074/1‑93 (MQ=255)
    aGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGt                                                                                    >  1:4670016/1‑93 (MQ=255)
    aGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGt                                                                                    >  2:2975132/1‑93 (MQ=255)
     ggTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtc                                                                                   >  2:1107207/1‑93 (MQ=255)
     ggTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtc                                                                                   <  2:458304/93‑1 (MQ=255)
      gTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtct                                                                                  <  1:3236525/93‑1 (MQ=255)
      gTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtct                                                                                  <  1:3041087/93‑1 (MQ=255)
      gTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtct                                                                                  >  2:3178192/1‑93 (MQ=255)
      gTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtct                                                                                  <  1:4347282/93‑1 (MQ=255)
       ttGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtctc                                                                                 <  2:3707012/93‑1 (MQ=255)
       ttGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGtctc                                                                                 <  1:1524018/93‑1 (MQ=255)
        tGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCg                                                                                >  2:4984040/1‑93 (MQ=255)
        tGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCg                                                                                <  2:3580974/93‑1 (MQ=255)
         tcctctGTTAACCCCCGATACGCACAGATAAAAAACCACACGGTCCATTCTGCTGCAACCCAAGCGGGTGTTGGCCAGGGCCATACGTCGCGt                                                                               <  1:173348/91‑1 (MQ=255)
         gTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGt                                                                               >  1:3884016/1‑93 (MQ=255)
         gTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGt                                                                               <  1:5938446/93‑1 (MQ=255)
         gTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGt                                                                               >  1:6409344/1‑93 (MQ=255)
          tctcTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGtt                                                                              <  1:1902233/93‑1 (MQ=255)
            tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            <  2:1212119/93‑1 (MQ=255)
            tctGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTAt                                                                            <  1:6604878/93‑1 (MQ=255)
               gTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGa                                                                         >  2:2545876/1‑93 (MQ=255)
                 tAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAcc                                                                       >  2:2147874/1‑93 (MQ=255)
                  aaCACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGAccc                                                                      <  1:4526695/93‑1 (MQ=255)
                   acacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCg                                                                     >  1:1421402/1‑93 (MQ=255)
                    cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    <  1:972017/93‑1 (MQ=255)
                    cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    <  2:3318166/93‑1 (MQ=255)
                    cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    <  2:5297643/93‑1 (MQ=255)
                    cacCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGt                                                                    >  1:1285895/1‑93 (MQ=255)
                       ccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTaa                                                                 <  2:3663294/93‑1 (MQ=255)
                       ccGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTaa                                                                 <  2:6673350/93‑1 (MQ=255)
                         gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               >  2:721796/1‑93 (MQ=255)
                         gATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAata                                                               <  2:5491672/93‑1 (MQ=255)
                           ttCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATAtt                                                             <  2:537384/93‑1 (MQ=255)
                             cacacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTgg                                                           <  2:1654455/93‑1 (MQ=255)
                              acacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGt                                                          <  1:4302373/93‑1 (MQ=255)
                              acacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGt                                                          >  2:414237/1‑93 (MQ=255)
                               cacaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTg                                                         >  1:2410906/1‑93 (MQ=255)
                                acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        <  2:3643516/93‑1 (MQ=255)
                                acaGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGc                                                        >  1:5792038/1‑93 (MQ=255)
                                   gATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGaa                                                     <  2:5842521/93‑1 (MQ=255)
                                    aTAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAac                                                    >  2:4266339/1‑93 (MQ=255)
                                      aaaaaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAacac                                                  >  1:5062601/1‑93 (MQ=255)
                                       aaaaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcacc                                                 >  2:4154945/1‑93 (MQ=255)
                                        aaaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcacca                                                >  2:4214956/1‑93 (MQ=255)
                                         aaaCAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAAcaccac                                               >  2:6554397/1‑93 (MQ=255)
                                          aacaacACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACaccacc                                              >  2:1957709/1‑93 (MQ=255)
                                          aacaacACGGTATATTCTGCTGCAACCCAAACGGGCGGTTGCCAGGGCGATACGTCTCGTTATGGACCCGCTAATATTGGTGCGAACaccacc                                              >  1:1134530/1‑93 (MQ=255)
                                           acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCt                                             <  1:5575468/93‑1 (MQ=255)
                                           acaacaCGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTACGAACACCACCt                                             >  1:2584061/1‑93 (MQ=255)
                                              acacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTa                                          >  2:5029599/1‑93 (MQ=255)
                                              acacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTa                                          <  2:2040742/93‑1 (MQ=255)
                                               cacGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAc                                         <  1:177396/93‑1 (MQ=255)
                                                 cGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTAccc                                       <  1:4916164/93‑1 (MQ=255)
                                                  ggTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCt                                      >  1:3180915/1‑93 (MQ=255)
                                                   gTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTg                                     >  1:1049352/1‑93 (MQ=255)
                                                    tataTTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTga                                    >  1:858465/1‑92 (MQ=255)
                                                     atatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTa                                   <  2:1913986/93‑1 (MQ=255)
                                                      tatTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTAt                                  >  2:177468/1‑93 (MQ=255)
                                                       atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 >  1:528439/1‑93 (MQ=255)
                                                       atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 <  2:3433405/93‑1 (MQ=255)
                                                       atTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATg                                 >  2:556864/1‑93 (MQ=255)
                                                         tCTGCTGCAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTCGTGCGAACACCACCTTTACCCTGTATGTc                               >  1:6223492/1‑93 (MQ=255)
                                                            gctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAcc                            >  2:6261367/1‑93 (MQ=255)
                                                            gctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCAcc                            >  2:4795447/1‑93 (MQ=255)
                                                             ctgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCa                           >  2:3942965/1‑93 (MQ=255)
                                                              tgcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCaa                          <  2:2121786/93‑1 (MQ=255)
                                                               gcAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAg                         >  2:3156150/1‑93 (MQ=255)
                                                                cAACCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGc                        <  2:4055616/93‑1 (MQ=255)
                                                                  aCCCAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCa                      >  1:4731182/1‑93 (MQ=255)
                                                                   cccAAACGGGTGTTGGCCAGTGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAt                     <  1:1947794/93‑1 (MQ=255)
                                                                    ccAAACGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCAtt                    <  2:6470510/93‑1 (MQ=255)
                                                                      aaaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtc                  >  1:5206131/1‑93 (MQ=255)
                                                                       aaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtct                 >  1:3559199/1‑93 (MQ=255)
                                                                       aaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtct                 >  1:2899517/1‑93 (MQ=255)
                                                                       aaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtct                 >  1:2729311/1‑93 (MQ=255)
                                                                       aaCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtct                 >  2:159676/1‑93 (MQ=255)
                                                                        aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  2:4333331/93‑1 (MQ=255)
                                                                        aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  2:155506/93‑1 (MQ=255)
                                                                        aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  1:6258610/93‑1 (MQ=255)
                                                                        aCGGGTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTtctc                <  1:3810585/93‑1 (MQ=255)
                                                                           ggTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGc             >  1:396824/1‑93 (MQ=255)
                                                                           ggTGTTTGCCAGGGCGATACGTCTCGTTATGGACCCGGTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGc             >  1:482656/1‑93 (MQ=255)
                                                                           ggTGTTGGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGc             <  1:6707785/93‑1 (MQ=255)
                                                                            gtgtTTGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCt            >  2:3544556/1‑93 (MQ=255)
                                                                              gtTTGCCAGGGCCATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCg          >  2:1406071/1‑93 (MQ=255)
                                                                                ttGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGAt        <  1:5113109/93‑1 (MQ=255)
                                                                                 tGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACCCCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATg       <  1:3586657/93‑1 (MQ=255)
                                                                                 tGCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATg       <  1:3764572/93‑1 (MQ=255)
                                                                                  gCCAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGa      <  1:1028803/93‑1 (MQ=255)
                                                                                   ccAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGAc     >  1:5406475/1‑93 (MQ=255)
                                                                                    cAGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGAcc    >  2:473485/1‑93 (MQ=255)
                                                                                     aGGGCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCa   <  2:2270450/93‑1 (MQ=255)
                                                                                      gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  >  2:508485/1‑93 (MQ=255)
                                                                                      gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  >  2:5976367/1‑93 (MQ=255)
                                                                                      gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  >  1:4431695/1‑93 (MQ=255)
                                                                                      gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAt  >  1:6106687/1‑93 (MQ=255)
                                                                                      gggCGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGACATTTCTCGGCTCGCTTACCAt  >  2:6368397/1‑93 (MQ=255)
                                                                                         |                                                                                         
GACCAGGTTGTCTCTGTTAACACCCGATTCACACAGATAAAAAACAACACGGTATATTCTGCTGCAACCCAAACGGGTGTTTGCCAGGGTGATACGTCTCGTTATGGACCCGTTAATATTGGTGCGAACACCACCTTTACCCTGTATGTCACCAAGCCATTTCTCGGCTCGATGACCAT  >  NC_000913/3289342‑3289520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: