Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 156,869 T→G L5F (TTA→TTC yadN ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913156,8690TG97.3% 263.2 / ‑0.1 75L5F (TTA→TTCyadNputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (38/35);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

GATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGTT  >  NC_000913/156722‑157016
                                                                                                                                                   |                                                                                                                                                   
gATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAt                                                                                                                                                  <  1:918487/151‑1 (MQ=255)
      aCCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAAGCATCGCCAGCATAAGGCCAGACAGGGCAAAACCGAAtttt                                                                                                                                               <  2:202003/148‑1 (MQ=255)
      aCCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAtttttt                                                                                                                                             >  2:651461/1‑150 (MQ=255)
             ttGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATa                                                                                                                                      >  1:928010/1‑150 (MQ=255)
              tGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATaa                                                                                                                                     >  2:259761/1‑150 (MQ=255)
              tGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATa                                                                                                                                      <  1:1164362/149‑1 (MQ=255)
               gTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATaaa                                                                                                                                    <  1:93110/150‑1 (MQ=255)
               gTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATa                                                                                                                                      <  1:616861/148‑1 (MQ=255)
                ttACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATaa                                                                                                                                     <  2:836041/148‑1 (MQ=255)
                 tACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAAt                                                                                                                                 >  2:509331/1‑151 (MQ=255)
                  aCCGCCGTCAACGCGGGTTTCACAGGTGTTTTCAACAACCAGACCACTGATAATTAACTGACCGCCGTCCATATCAGAAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATGAAAATc                                                                                                                                <  1:683170/151‑1 (MQ=255)
                   ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATc                                                                                                                                <  1:1143391/150‑1 (MQ=255)
                   ccgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTAACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATcc                                                                                                                               <  2:500509/151‑1 (MQ=255)
                    cgccgTCAACGCGGGTTTCACAGGAGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACAAGACAGCGCAAAACCGAATTTTTTAGACATAAAAATCCt                                                                                                                              >  2:757712/1‑151 (MQ=255)
                       cgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCttt                                                                                                                            >  2:444491/1‑150 (MQ=255)
                       cgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTta                                                                                                                           <  1:679073/151‑1 (MQ=255)
                        gTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaa                                                                                                                          <  1:432926/151‑1 (MQ=255)
                          cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaata                                                                                                                        <  1:211717/151‑1 (MQ=255)
                          cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaat                                                                                                                         >  2:714790/1‑150 (MQ=255)
                             cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATaaa                                                                                                                      >  2:923735/1‑150 (MQ=255)
                              gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAAt                                                                                                                    <  1:550102/151‑1 (MQ=255)
                              gcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGCCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAAt                                                                                                                    >  2:616253/1‑151 (MQ=255)
                               cgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAAtt                                                                                                                   >  1:744399/1‑151 (MQ=255)
                                gggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTc                                                                                                                  >  2:664242/1‑151 (MQ=255)
                                       acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGt                                                                                                           >  1:647872/1‑151 (MQ=255)
                                       acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGt                                                                                                           >  1:836781/1‑151 (MQ=255)
                                        caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTTCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTc                                                                                                          <  2:752404/151‑1 (MQ=255)
                                              ttATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACt                                                                                                     >  2:60580/1‑150 (MQ=255)
                                                  caacaaCCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACAt                                                                                                 <  1:431868/150‑1 (MQ=255)
                                                     caacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGttt                                                                                             >  1:39340/1‑151 (MQ=255)
                                                      aacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGttt                                                                                             >  2:316406/1‑150 (MQ=255)
                                                      aacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTg                                                                                            >  2:1141872/1‑151 (MQ=255)
                                                        ccAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTga                                                                                           >  2:183659/1‑150 (MQ=255)
                                                           gACCACTGATATTTATCTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTCCATGTTTgaaga                                                                                        <  2:731622/150‑1 (MQ=255)
                                                                 tGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACt                                                                                  >  1:893926/1‑150 (MQ=255)
                                                                  gATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTAt                                                                                >  1:979901/1‑151 (MQ=255)
                                                                  gATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTAt                                                                                >  2:1003560/1‑151 (MQ=255)
                                                                    tatTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTc                                                                              >  2:302718/1‑151 (MQ=255)
                                                                     atTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTc                                                                              <  1:821748/150‑1 (MQ=255)
                                                                      tttAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCAt                                                                            >  2:232226/1‑151 (MQ=255)
                                                                       ttAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCAt                                                                            <  2:433594/150‑1 (MQ=255)
                                                                         aaCTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGac                                                                         >  1:1170737/1‑151 (MQ=255)
                                                                         aaCTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGa                                                                          <  1:1212339/150‑1 (MQ=255)
                                                                               ccgccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAAt                                                                   <  1:74524/151‑1 (MQ=255)
                                                                                cgccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCATACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAAt                                                                   >  1:172370/1‑150 (MQ=255)
                                                                                 gccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAAt                                                                   >  2:169852/1‑149 (MQ=255)
                                                                                      ccATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGaa                                                            <  1:158532/151‑1 (MQ=255)
                                                                                         tatCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAAc                                                          >  1:344284/1‑150 (MQ=255)
                                                                                               cagAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGtt                                                   >  2:442252/1‑151 (MQ=255)
                                                                                               cagAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGtt                                                   >  2:20657/1‑151 (MQ=255)
                                                                                               cagAAGCAGTACATGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGtt                                                   <  1:370843/151‑1 (MQ=255)
                                                                                                   aGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGat                                               <  1:1109364/151‑1 (MQ=255)
                                                                                                      aGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATAtt                                            <  2:1057073/151‑1 (MQ=255)
                                                                                                       gTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAAAAGGAGAAACAAATGTTAGATATTa                                           <  1:993485/151‑1 (MQ=255)
                                                                                                         aCCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACTCAAATAGGAGAAACAAATGTTAGATATTAAt                                         <  1:825344/151‑1 (MQ=255)
                                                                                                           cTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATg                                        <  2:317838/150‑1 (MQ=255)
                                                                                                               aaCCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCaa                                   <  2:408252/151‑1 (MQ=255)
                                                                                                               aaCCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCaa                                   <  1:243805/151‑1 (MQ=255)
                                                                                                               aaCCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCa                                    >  2:701823/1‑150 (MQ=255)
                                                                                                                aCCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAAt                                  >  1:715240/1‑151 (MQ=255)
                                                                                                                  aaTCGCCAGCAAAAGAGAAGGAAGGGCAAAACCGGAATTATTAGACATAAAAATTCTTTAATAAAATTCCATTAGTCAAACTACATGTTTGAAGCATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAAt                                  >  2:903315/2‑149 (MQ=255)
                                                                                                                   aTCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGat                               >  1:830292/1‑151 (MQ=255)
                                                                                                                     cGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGata                              <  1:494633/150‑1 (MQ=255)
                                                                                                                      gCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATAtt                            >  1:587689/1‑151 (MQ=255)
                                                                                                                      gCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATAtt                            >  1:931827/1‑151 (MQ=255)
                                                                                                                       ccAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATAttt                           <  1:923795/151‑1 (MQ=255)
                                                                                                                        cAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTg                          >  2:358967/1‑151 (MQ=255)
                                                                                                                         aGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGt                         <  2:219002/151‑1 (MQ=255)
                                                                                                                           cATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGtt                        >  2:439379/1‑150 (MQ=255)
                                                                                                                              aaGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTAccc                    >  2:25784/1‑151 (MQ=255)
                                                                                                                                  ccAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAAt                <  1:587398/151‑1 (MQ=255)
                                                                                                                                     gACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAAttt              >  1:1133467/1‑150 (MQ=255)
                                                                                                                                           gCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCa       <  1:1185807/151‑1 (MQ=255)
                                                                                                                                           gCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCa       <  1:283378/151‑1 (MQ=255)
                                                                                                                                                aCCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGtt  <  1:259772/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGTT  >  NC_000913/156722‑157016

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: