Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,253,830 C→G 12.2% S708T (AGT→ACT)  ycgV ← putative adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,253,8300CG12.2% 312.6 / 30.5 98S708T (AGT→ACT) ycgVputative adhesin
Reads supporting (aligned to +/- strand):  ref base C (44/42);  new base G (10/2);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 6.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

TGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCA  >  NC_000913/1253682‑1253976
                                                                                                                                                    |                                                                                                                                                  
tGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAc                                                                                                                                                    >  1:67485/1‑149 (MQ=255)
tGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCAc                                                                                                                                                    >  4:640608/1‑149 (MQ=255)
  cGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                  <  2:679184/149‑1 (MQ=255)
    gTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGc                                                                                                                                                >  4:344527/1‑149 (MQ=255)
       ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                             <  4:611694/149‑1 (MQ=255)
       ccGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCaaa                                                                                                                                             <  1:391043/149‑1 (MQ=255)
        cGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            <  3:85872/149‑1 (MQ=255)
        cGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            <  2:67485/149‑1 (MQ=255)
          tCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGa                                                                                                                                          >  2:268950/1‑149 (MQ=255)
                ccGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAg                                                                                                                                    <  1:27024/149‑1 (MQ=255)
                 cGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAgg                                                                                                                                   <  3:503069/149‑1 (MQ=255)
                      taGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTcc                                                                                                                              <  1:403569/149‑1 (MQ=255)
                       aGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTccc                                                                                                                             <  3:679799/149‑1 (MQ=255)
                       aGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTccc                                                                                                                             <  4:477848/149‑1 (MQ=255)
                           cGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGc                                                                                                                                  >  4:383976/1‑140 (MQ=255)
                           cGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGc                                                                                                                                  <  3:383976/140‑1 (MQ=255)
                           cGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCAt                                                                                                                         >  3:351418/1‑149 (MQ=255)
                             gCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATaa                                                                                                                       <  3:640608/149‑1 (MQ=255)
                                 tGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                   <  2:20665/149‑1 (MQ=255)
                                   cATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                   >  3:545198/1‑147 (MQ=255)
                                   cATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTgc                                                                                                                   <  4:545198/147‑1 (MQ=255)
                                     tgtgTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCaa                                                                                                               >  4:240557/1‑149 (MQ=255)
                                       tgtTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            >  4:226611/1‑118 (MQ=255)
                                       tgtTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            <  3:226611/118‑1 (MQ=255)
                                        gtTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCa                                                                                                            >  4:193760/1‑149 (MQ=255)
                                           gAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGat                                                                                                         <  3:438996/149‑1 (MQ=255)
                                           gAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGat                                                                                                         >  4:295291/1‑149 (MQ=255)
                                             gCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGatat                                                                                                       >  2:222406/1‑149 (MQ=255)
                                              ccAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATAtt                                                                                                      >  1:273871/1‑149 (MQ=255)
                                                aaTATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTAc                                                                                                    >  1:512547/1‑149 (MQ=255)
                                                      cAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCtt                                                                                              <  2:251650/149‑1 (MQ=255)
                                                        gACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTa                                                                                            >  2:324818/1‑149 (MQ=255)
                                                         aCCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                  >  3:388922/1‑94 (MQ=255)
                                                         aCCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTg                                                                                                                                                  <  4:388922/94‑1 (MQ=255)
                                                         aCCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTAc                                                                                           >  3:564397/1‑149 (MQ=255)
                                                          ccGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTAct                                                                                          <  1:32050/149‑1 (MQ=255)
                                                            gACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGCTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTActct                                                                                        >  3:139720/1‑149 (MQ=255)
                                                               aTACAAGGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAt                                                                                     >  2:490855/1‑149 (MQ=255)
                                                               aTACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAt                                                                                     >  4:174928/1‑149 (MQ=255)
                                                                tACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAtt                                                                                    <  4:288838/149‑1 (MQ=255)
                                                                 aCAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGAttt                                                                                   <  4:658725/149‑1 (MQ=255)
                                                                  cAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTc                                                                                  <  3:295291/149‑1 (MQ=255)
                                                                   aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAACGCCGCTCAGTTTGCCACTGGCCACGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTCCTCTGCTTTCg                                                                                 >  4:115705/1‑149 (MQ=255)
                                                                     cGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                               <  3:193760/149‑1 (MQ=255)
                                                                        cATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGt                                                                            >  4:677471/1‑149 (MQ=255)
                                                                              aTCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCa                                                                      >  4:460663/1‑149 (MQ=255)
                                                                               tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            <  3:28022/78‑1 (MQ=255)
                                                                               tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAg                                                                                                                                            >  4:28022/1‑78 (MQ=255)
                                                                               tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCAt                                                                     >  2:429112/1‑149 (MQ=255)
                                                                                cAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATc                                                                    <  1:429112/149‑1 (MQ=255)
                                                                                    gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                <  3:677471/149‑1 (MQ=255)
                                                                                    gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                >  4:297386/1‑149 (MQ=255)
                                                                                         gTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATc                                                           >  2:107372/1‑149 (MQ=255)
                                                                                             aTCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCg                                                       >  1:187560/1‑149 (MQ=255)
                                                                                              tCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCTGTTGCATCAGCGt                                                      <  1:268950/149‑1 (MQ=255)
                                                                                                cccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATaa                                                                                                                       <  4:172490/82‑1 (MQ=255)
                                                                                                cccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATaa                                                                                                                       >  3:172490/1‑82 (MQ=255)
                                                                                                cccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTAc                                                    <  3:240557/149‑1 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCGTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATACCTGAGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCAt                                                            <  4:419780/138‑1 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                               >  1:504388/1‑119 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCa                                                                               <  2:504388/119‑1 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCAt                                                            >  3:419780/1‑138 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                 >  4:132133/1‑149 (MQ=255)
                                                                                                   accaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGt                                                 >  1:103983/1‑149 (MQ=255)
                                                                                                    ccaAACTGGATACCGCTGTAACCCATGTCAAACCCACTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGtt                                                <  3:276006/149‑1 (MQ=255)
                                                                                                        aCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCa                                            <  3:174928/149‑1 (MQ=255)
                                                                                                           ggATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACa                                         <  2:273871/149‑1 (MQ=255)
                                                                                                             aTACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                       <  1:490855/149‑1 (MQ=255)
                                                                                                                ccGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGtt                                    <  1:107372/149‑1 (MQ=255)
                                                                                                                ccGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGtt                                    <  4:139720/149‑1 (MQ=255)
                                                                                                                  gCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCa                                  <  4:351418/149‑1 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTc                                                                                                <  1:218884/79‑1 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTc                                                                                                >  2:218884/1‑79 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTc                                                                           >  4:78787/1‑100 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTc                                                                           <  3:78787/100‑1 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                          >  3:604658/1‑149 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                          >  2:338414/1‑149 (MQ=255)
                                                                                                                          ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTaa                          >  4:58035/1‑149 (MQ=255)
                                                                                                                                 aaaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAt                   >  3:88741/1‑149 (MQ=255)
                                                                                                                                       gCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGAt             >  4:305939/1‑149 (MQ=255)
                                                                                                                                       gCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGAt             >  4:481862/1‑149 (MQ=255)
                                                                                                                                         tCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATaa           <  4:604658/149‑1 (MQ=255)
                                                                                                                                         ctcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATaa           >  4:576155/5‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCCACATAGTTCAATAAGTAACCGACATTGAGGTAAt          >  3:226413/4‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                       >  1:348044/4‑120 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATa                                       <  2:348044/117‑1 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          >  1:453653/4‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          >  2:342823/4‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          >  3:447695/4‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          >  3:612109/4‑149 (MQ=255)
                                                                                                                                          tcctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          >  2:146939/4‑149 (MQ=255)
                                                                                                                                          cAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCATATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAt          <  2:512547/149‑1 (MQ=255)
                                                                                                                                           cctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAtt         <  4:612109/147‑1 (MQ=255)
                                                                                                                                            ctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTg        >  2:65694/2‑149 (MQ=255)
                                                                                                                                            ctttGCCAGTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTg        >  1:315125/2‑149 (MQ=255)
                                                                                                                                                 ccACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTgccgcc   >  2:566987/1‑149 (MQ=255)
                                                                                                                                                  cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa  >  4:375598/1‑149 (MQ=255)
                                                                                                                                                  cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa  >  3:103123/1‑149 (MQ=255)
                                                                                                                                                  cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa  >  3:646969/1‑149 (MQ=255)
                                                                                                                                                  cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa  >  4:151513/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCGGTACCGTCGCCTCCGCTATAGTCCGGCGATGCATGTGTTGAGCCAATATACAGACCGACATACAACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCA  >  NC_000913/1253682‑1253976

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