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.7% 327.7 / ‑1.7 87L5F (TTA→TTCyadNputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (34/51);  total (35/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

GATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGT  >  NC_000913/156722‑157015
                                                                                                                                                   |                                                                                                                                                  
gATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAt                                                                                                                                                 <  1:60403/151‑1 (MQ=255)
  tCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAttt                                                                                                                                               >  1:400991/1‑151 (MQ=255)
    aGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAttttt                                                                                                                                             >  2:1001984/1‑151 (MQ=255)
     gACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAtttttt                                                                                                                                            <  2:1201516/151‑1 (MQ=255)
       ccATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAAtttttt                                                                                                                                            >  2:532749/1‑149 (MQ=255)
       ccATCCTTGTGACCGCTGTCAACGCGGGTTGGACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTa                                                                                                                                           <  1:123046/150‑1 (MQ=255)
        cATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGa                                                                                                                                         <  2:164632/151‑1 (MQ=255)
        aaTCCTTGGTACTGCCGACAACGGGGGTTTCACAAGCGTGATCAACAACCAGACCACTGATATTTAACTGAGCGCGGACCAGATCAGCAAAAGCAGTACCAGAAACCATCGCGAGAAAAAGAACACACAGTTCAAAACCGAAttt                                                                                                                                               >  1:1066079/2‑145 (MQ=255)
           ccTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACAt                                                                                                                                      >  2:198124/1‑151 (MQ=255)
             ttGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATa                                                                                                                                     >  2:173530/1‑150 (MQ=255)
              tGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATa                                                                                                                                     <  2:1073590/149‑1 (MQ=255)
                    cgccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCt                                                                                                                             <  1:779869/151‑1 (MQ=255)
                     gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCtt                                                                                                                            >  1:226842/1‑151 (MQ=255)
                     gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCtt                                                                                                                            >  1:623965/1‑151 (MQ=255)
                     gccgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCtt                                                                                                                            >  2:737777/1‑151 (MQ=255)
                       cgTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTta                                                                                                                          <  2:147834/151‑1 (MQ=255)
                          cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaat                                                                                                                        >  1:542082/1‑150 (MQ=255)
                          cAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaat                                                                                                                        <  2:994844/150‑1 (MQ=255)
                           aaCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTtaata                                                                                                                       >  1:470774/1‑150 (MQ=255)
                            aCGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATaaa                                                                                                                     <  1:1102648/151‑1 (MQ=255)
                             cgcgGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATaaa                                                                                                                     >  1:60398/1‑150 (MQ=255)
                                gggTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTc                                                                                                                 >  2:1095658/1‑151 (MQ=255)
                                      cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTa                                                                                                            <  2:934537/150‑1 (MQ=255)
                                      cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTa                                                                                                            <  1:574610/150‑1 (MQ=255)
                                      cacaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAg                                                                                                           <  1:1252422/151‑1 (MQ=255)
                                       acaGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAg                                                                                                           <  2:1056307/150‑1 (MQ=255)
                                        caGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTc                                                                                                         >  1:425565/1‑151 (MQ=255)
                                          ggTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAg                                                                                                       <  2:701386/151‑1 (MQ=255)
                                           gtgtTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTc                                                                                                         >  1:754224/1‑148 (MQ=255)
                                               tATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTa                                                                                                   <  1:1038373/150‑1 (MQ=255)
                                               tATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTa                                                                                                   <  1:186090/150‑1 (MQ=255)
                                                aTCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACa                                                                                                 <  1:652113/151‑1 (MQ=255)
                                                  caacaaCCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACAt                                                                                                <  1:1142929/150‑1 (MQ=255)
                                                  caacaaCCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACAt                                                                                                <  1:186160/150‑1 (MQ=255)
                                                   aacaacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGt                                                                                              >  2:546573/1‑151 (MQ=255)
                                                    acaacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGtt                                                                                             <  2:1117767/151‑1 (MQ=255)
                                                    acaacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGACAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGtt                                                                                             <  1:1192747/151‑1 (MQ=255)
                                                      aacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTg                                                                                           <  1:533104/151‑1 (MQ=255)
                                                      aacCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTg                                                                                           <  2:827110/151‑1 (MQ=255)
                                                              cACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGa                                                                                   <  2:103016/151‑1 (MQ=255)
                                                                 tGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACt                                                                                 >  1:940314/1‑150 (MQ=255)
                                                                   atatTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTAtt                                                                              <  2:597849/151‑1 (MQ=255)
                                                                   atatTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTAtt                                                                              <  1:842909/151‑1 (MQ=255)
                                                                     atTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCa                                                                            <  1:717793/151‑1 (MQ=255)
                                                                      tttAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCAt                                                                           >  2:952679/1‑151 (MQ=255)
                                                                          aCTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGac                                                                        >  2:978374/1‑150 (MQ=255)
                                                                             gACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTTATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGAcaca                                                                     >  1:248213/1‑150 (MQ=255)
                                                                             gACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGacac                                                                      <  1:224453/149‑1 (MQ=255)
                                                                               ccgccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAAt                                                                  >  2:213557/1‑151 (MQ=255)
                                                                               ccgacgACCAAATCAGCAGAAGCAGTACCTGCAAACAACGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAAAAAAATTCCATTAGTCAGACTACATGTGTGAAGAATGACTATTCATGACACaa                                                                    <  1:565507/149‑1 (MQ=255)
                                                                                cgccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATa                                                                 <  1:1139168/151‑1 (MQ=255)
                                                                                  ccgTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAgg                                                               <  2:446524/151‑1 (MQ=255)
                                                                                    gTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAgg                                                               >  1:386749/1‑149 (MQ=255)
                                                                                       cATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGaaa                                                          <  2:1143086/151‑1 (MQ=255)
                                                                                           tCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACa                                                        >  1:352417/1‑149 (MQ=255)
                                                                                            cagcagAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAAt                                                     <  2:470787/151‑1 (MQ=255)
                                                                                             agcagAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATg                                                    <  1:960478/151‑1 (MQ=255)
                                                                                             agcagAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATg                                                    <  1:1053184/151‑1 (MQ=255)
                                                                                                agAAGCAGTACCTGCCACCATCGCCAGCCCAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCCTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGtt                                                  >  1:1104142/1‑150 (MQ=255)
                                                                                                agAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTa                                                 <  1:1028410/151‑1 (MQ=255)
                                                                                                 gAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGtt                                                  >  1:179502/1‑149 (MQ=255)
                                                                                                    gCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGata                                             <  1:555354/151‑1 (MQ=255)
                                                                                                     cAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGatat                                            >  2:528764/1‑151 (MQ=255)
                                                                                                     cAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGatat                                            <  2:994629/151‑1 (MQ=255)
                                                                                                     cAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGatat                                            <  2:101531/151‑1 (MQ=255)
                                                                                                         aCCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTaa                                         <  1:198122/150‑1 (MQ=255)
                                                                                                         aCCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAAt                                        >  1:109819/1‑151 (MQ=255)
                                                                                                         aCCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAAt                                        <  2:350315/151‑1 (MQ=255)
                                                                                                             gCAACCATCGCCAGCATAAGTCCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGc                                    >  2:378146/1‑151 (MQ=255)
                                                                                                                 ccATCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATg                                <  2:930861/151‑1 (MQ=255)
                                                                                                                    tCGCCAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGat                              >  2:1008155/1‑150 (MQ=255)
                                                                                                                       ccAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATAttt                          <  2:178797/151‑1 (MQ=255)
                                                                                                                       ccAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATAttt                          <  1:570045/151‑1 (MQ=255)
                                                                                                                       ccAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGCTAGATATTAATGAGCAATGATAttt                          <  1:364040/151‑1 (MQ=255)
                                                                                                                        cAGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTg                         >  1:764073/1‑151 (MQ=255)
                                                                                                                         aGCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGt                        >  2:817824/1‑151 (MQ=255)
                                                                                                                          gCATAAGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGtt                       <  1:1263077/151‑1 (MQ=255)
                                                                                                                               aGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTAccc                   >  2:370790/1‑150 (MQ=255)
                                                                                                                               aGACCAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCa                  <  2:600917/151‑1 (MQ=255)
                                                                                                                                  ccAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAAt               <  2:319289/151‑1 (MQ=255)
                                                                                                                                   cAGACAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAAtt              <  1:602261/151‑1 (MQ=255)
                                                                                                                                      aCAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTAc           >  1:587607/1‑151 (MQ=255)
                                                                                                                                      aCAGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTAc           <  2:221525/151‑1 (MQ=255)
                                                                                                                                        aGGGCAAAACCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACaa         <  2:434375/151‑1 (MQ=255)
                                                                                                                                             aaaaCCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCAtt    >  1:411748/1‑151 (MQ=255)
                                                                                                                                              aaaCCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTg   <  2:765031/151‑1 (MQ=255)
                                                                                                                                               aaCCGAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGt  >  1:483154/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATCAGACCATCCTTGTTACCGCCGTCAACGCGGGTTTCACAGGTGTTATCAACAACCAGACCACTGATATTTAACTGACCGCCGTCCATATCAGCAGAAGCAGTACCTGCAACCATCGCCAGCATAAGACCAGACAGGGCAAAACCTAATTTTTTAGACATAAAAATCCTTTAATAAAATTCCATTAGTCAGACTACATGTTTGAAGAATGACTATTCATGACACAAATAGGAGAAACAAATGTTAGATATTAATGAGCAATGATATTTGTTACCCAAATTTACAACCATTGT  >  NC_000913/156722‑157015

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