Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,908 C→G 5.9% intergenic (+17/+22) yfeN → / ← yfeR putative outer membrane protein/transcriptional regulator of yefH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,9080CG 5.9% 593.6 / 20.1 169intergenic (+17/+22)yfeN/yfeRputative outer membrane protein/transcriptional regulator of yefH
Reads supporting (aligned to +/- strand):  ref base C (83/76);  new base G (5/5);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 174 total aligned reads displayed.

GGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525760‑2526056
                                                                                                                                                     |                                                                                                                                                    
ggCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTGTCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑Gc                                                                                                                                                      >  1:463222/1‑149 (MQ=255)
 gcgcAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCt                                                                                                                                                     >  2:523575/1‑149 (MQ=255)
        tGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCa                                                                                                                                              <  2:1442737/149‑1 (MQ=255)
          cACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATTGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATGTAACTTCTGATCGATTACTGGCGGG‑GCTCCCGGCAGt                                                                                                                                            <  2:55495/149‑1 (MQ=255)
           aCTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTa                                                                                                                                           <  4:1195541/149‑1 (MQ=255)
            cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTAt                                                                                                                                          <  3:1305681/149‑1 (MQ=255)
            cTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTaa                                                                                                                                          <  3:572051/149‑2 (MQ=255)
               tggtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                       <  1:1451919/149‑1 (MQ=255)
                 gtggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAAtt                                                                                                                                     <  1:859867/149‑1 (MQ=255)
                   ggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCa                                                                                                                                   >  1:832326/1‑149 (MQ=255)
                      aTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                                <  1:566562/149‑1 (MQ=255)
                      aTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                                <  1:1009385/149‑1 (MQ=255)
                      aTGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActc                                                                                                                                >  3:4617/1‑149 (MQ=255)
                       tGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActct                                                                                                                               >  4:156024/1‑149 (MQ=255)
                          aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGC‑CCCGGCAGTaatc                                                                                                                                        <  3:710499/137‑4 (MQ=255)
                          aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTaatc                                                                                                                                        >  4:710499/1‑134 (MQ=255)
                          aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACATCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCtata                                                                                                                            >  3:1199929/1‑149 (MQ=255)
                             aTCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCt                                                                                                                         <  4:951579/149‑1 (MQ=255)
                              tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGaa                                                                                                                                       >  2:1445432/1‑134 (MQ=255)
                              tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActct                                                                                                                               <  1:306535/142‑1 (MQ=255)
                              tCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCActct                                                                                                                               >  2:306535/1‑142 (MQ=255)
                               cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  1:28350/149‑1 (MQ=255)
                               cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  1:808172/149‑1 (MQ=255)
                               cacacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       >  3:1237683/1‑149 (MQ=255)
                                acacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGAt                                                                                                                      <  4:397260/149‑1 (MQ=255)
                                 cacaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCCGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATa                                                                                                                     <  2:770463/149‑1 (MQ=255)
                                  acaTTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATAc                                                                                                                    >  3:333239/1‑149 (MQ=255)
                                     ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                                 <  4:239713/149‑1 (MQ=255)
                                     ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAAc                                                                                                                 >  1:527750/1‑149 (MQ=255)
                                         tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGAt                                                                                                             >  1:218621/1‑149 (MQ=255)
                                          acGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                                            <  2:103935/149‑1 (MQ=255)
                                            ggggATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑AGCCCCGGCAGTaa                                                                                                                                          <  3:429486/117‑2 (MQ=255)
                                                   cAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTcc                                                                                                                                                   >  4:1076162/1‑101 (MQ=255)
                                                   cAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTcc                                                                                                                                                   <  3:1076162/101‑1 (MQ=255)
                                                    aGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGc                                                                                                  >  3:1261267/1‑149 (MQ=255)
                                                     gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTccc                                                                                                                                                  <  4:1028125/100‑1 (MQ=255)
                                                     gTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGcc                                                                                                 >  3:349290/1‑149 (MQ=255)
                                                          gATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGcc                                                                                            >  2:1370273/1‑149 (MQ=255)
                                                           atatGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCa                                                                                           <  2:1350028/149‑1 (MQ=255)
                                                              tGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                        <  3:233299/149‑1 (MQ=255)
                                                              tGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                        <  1:515826/149‑1 (MQ=255)
                                                              tGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTa                                                                                        <  4:536375/149‑1 (MQ=255)
                                                               gCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAg                                                                                       >  2:649859/1‑149 (MQ=255)
                                                                    tGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCAt                                                                                  >  1:739767/1‑149 (MQ=255)
                                                                        aaaCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAAGTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATGCCCTTCCCGCCCAGCCATTAGCGCATTGcc                                                                              >  1:958389/1‑149 (MQ=255)
                                                                          aCTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                            >  4:1154996/1‑149 (MQ=255)
                                                                           cTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                                           <  2:1169904/149‑1 (MQ=255)
                                                                             gggCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGt                                                                         <  1:1063339/149‑1 (MQ=255)
                                                                                 ccacctttttACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCg                                                                     >  4:85678/6‑149 (MQ=255)
                                                                                   aGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGTATTGCCCTGGTCTCGTa                                                                   <  4:492559/149‑1 (MQ=255)
                                                                                    gATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGAGCATTGCCCTGGTCTCGTAc                                                                  <  3:483501/149‑1 (MQ=255)
                                                                                          tACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTc                                                            <  3:121046/149‑1 (MQ=255)
                                                                                               ggATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATa                                                       >  3:102584/1‑149 (MQ=255)
                                                                                                 aTGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCCGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                     <  2:832326/149‑1 (MQ=255)
                                                                                                     gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                 >  1:92365/1‑149 (MQ=255)
                                                                                                     gATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCtt                                                 >  1:212940/1‑149 (MQ=255)
                                                                                                          tCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  4:254371/74‑1 (MQ=255)
                                                                                                          tCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcg                                            >  3:1404656/1‑149 (MQ=255)
                                                                                                          tCTACTCCAGCAAATTTAACGTCTGATTGATTACTGCCGCG‑GCTCCCGGCAGTCTGGATTCATTCTATATCTGATACAAAGGATCCCCTTCCCGCCCAGCCATTAGCTCATTGCCATGGTCTCGTACCACGTCCCATAAAGCTTCAgcg                                            >  4:633274/1‑149 (MQ=255)
                                                                                                            tACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAgcggc                                          >  4:45400/1‑149 (MQ=255)
                                                                                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGa                                                                                                                       <  3:1015884/70‑1 (MQ=255)
                                                                                                              cTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGcggcgg                                        >  4:961001/1‑149 (MQ=255)
                                                                                                                  aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCAct                                                                                                                                 >  4:790555/1‑56 (MQ=255)
                                                                                                                  aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCAct                                                                                                                                 <  3:790555/56‑1 (MQ=255)
                                                                                                                  aTCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGa                                    >  3:858475/1‑149 (MQ=255)
                                                                                                                   tCAAATTTAACTTCTGATTGATTACTTCCGGG‑GCTCCCCGCAGTATGAATTCACCCTATATCTGATTCAACGGATCCCCTTCCCCCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATACAGCCTCAGCGGCGGTGGaa                                   >  2:1234968/1‑149 (MQ=255)
                                                                                                                    cAAATTTACCTTCTTATTGAGTACTGCCGGT‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGCTCCCCTTCCCGCCCAGCCAGTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  2:463610/149‑1 (MQ=255)
                                                                                                                    cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  3:294179/149‑1 (MQ=255)
                                                                                                                    cAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  2:463222/149‑1 (MQ=255)
                                                                                                                      aaTTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGt                                >  1:420082/1‑149 (MQ=255)
                                                                                                                          tAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAAc                            <  4:676459/149‑1 (MQ=255)
                                                                                                                              ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCt                                                                            >  1:62592/1‑97 (MQ=255)
                                                                                                                              ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        >  3:689652/1‑149 (MQ=255)
                                                                                                                              ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        <  2:673865/149‑1 (MQ=255)
                                                                                                                              ttCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGtt                        <  4:1237683/149‑1 (MQ=255)
                                                                                                                               tCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGttc                       <  1:1137256/149‑1 (MQ=255)
                                                                                                                                 tgattgatTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  >  1:1237589/1‑136 (MQ=255)
                                                                                                                                 tgattgatTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  2:1237589/136‑1 (MQ=255)
                                                                                                                                  gattgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGttcttc                    >  1:1313998/1‑149 (MQ=255)
                                                                                                                                    ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                     >  1:1266436/1‑114 (MQ=255)
                                                                                                                                    ttgattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATaaa                                                     <  2:1266436/114‑1 (MQ=255)
                                                                                                                                      gattACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgc                <  4:102584/149‑1 (MQ=255)
                                                                                                                                        ttACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCA‑TAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa             >  2:1081418/1‑149 (MQ=255)
                                                                                                                                        atACTGCCGGGAGC‑CCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  >  1:1208262/2‑129 (MQ=255)
                                                                                                                                        atACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  2:1208262/128‑1 (MQ=255)
                                                                                                                                        atACTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTgcgc              >  2:154496/2‑149 (MQ=255)
                                                                                                                                         tACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCa             >  4:194651/1‑149 (MQ=255)
                                                                                                                                          aCTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCAc            >  4:592761/1‑149 (MQ=255)
                                                                                                                                           cTGCCGGG‑AGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGaaa                                  <  3:156930/126‑1 (MQ=255)
                                                                                                                                            tGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa          >  1:472348/1‑149 (MQ=255)
                                                                                                                                            tGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa          >  4:1318103/1‑149 (MQ=255)
                                                                                                                                            tGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCa          >  1:1204579/1‑149 (MQ=255)
                                                                                                                                            tGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGGGAACGGTTCTTCCTGCGCACCa          <  4:333239/149‑1 (MQ=255)
                                                                                                                                              ccGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc        >  1:1339629/1‑149 (MQ=255)
                                                                                                                                                ggg‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCAt      >  1:1440480/1‑149 (MQ=255)
                                                                                                                                                ggg‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCAt      >  2:1270712/1‑149 (MQ=255)
                                                                                                                                                 gg‑gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATc     <  4:1267795/149‑1 (MQ=255)
                                                                                                                                                    gCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCCTTAGCGCATTGCCATGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCaa   <  4:1219212/149‑1 (MQ=255)
                                                                                                                                                     cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACgg                          >  3:495376/1‑125 (MQ=255)
                                                                                                                                                     cTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAc  <  4:843187/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGCAGTTGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGG‑GCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCATCAAC  >  NC_000913/2525760‑2526056

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