Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,680,639 T→C 5.1% S703G (AGC→GGC)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,680,6390TC 5.1% 345.0 / 6.1 98S703G (AGC→GGC) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base T (46/47);  new base C (2/3);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

GCTAGTAAATGGCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGC  >  NC_000913/2680494‑2680786
                                                                                                                                                 |                                                                                                                                                   
gcTAGTAAATGGCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTgg                                                                                                                                                   >  2:38608/1‑148 (MQ=255)
     tAAATGGCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCCAGCTTGc                                                                                                                                              <  3:1016053/148‑1 (MQ=255)
        aTGGCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAgg                                                                                                                                           <  1:50609/148‑1 (MQ=255)
              aGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGc                                                                                                                                     >  4:297260/1‑148 (MQ=255)
                 aaaCCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGc                                                                                                                                  <  1:38608/148‑1 (MQ=255)
                  aaCCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCa                                                                                                                                 >  4:1004990/1‑148 (MQ=255)
                  aaCCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCa                                                                                                                                 >  4:123783/1‑148 (MQ=255)
                     cAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCg                                                                                                                              >  1:467038/1‑148 (MQ=255)
                       gCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTGCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTg                                                                                                                            <  2:467038/148‑1 (MQ=255)
                         aCTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTgcg                                                                                                                          >  2:404371/1‑148 (MQ=255)
                           tctTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCt                                                                                                                        <  3:123783/148‑1 (MQ=255)
                           tcatCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCt                                                                                                                        <  3:1004990/145‑1 (MQ=255)
                            ctTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCg                                                                                                                              <  4:639193/141‑1 (MQ=255)
                            ctTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCg                                                                                                                              >  3:639193/1‑141 (MQ=255)
                              tCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGc                                                                                                                     >  4:1118514/1‑148 (MQ=255)
                                   aCTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAg                                                                                                                                      <  1:583437/126‑1 (MQ=255)
                                   aCTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAg                                                                                                                                      >  2:583437/1‑126 (MQ=255)
                                      ctcCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCCAGCTTGctg                                                                                                                                            <  2:566083/117‑3 (MQ=255)
                                      ctcCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCCAGCTTGctg                                                                                                                                            >  1:566083/1‑115 (MQ=255)
                                                gccACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCa                                                                                                   <  3:297260/148‑1 (MQ=255)
                                                    cttcttCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                                                                                               <  3:1118514/148‑1 (MQ=255)
                                                              cGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                                                                                               >  1:68305/1‑138 (MQ=255)
                                                              cGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                                                                                               <  2:68305/138‑1 (MQ=255)
                                                                 aTTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTg                                                                                  <  1:404371/148‑1 (MQ=255)
                                                                   tCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTg                                                                                  <  3:539823/146‑1 (MQ=255)
                                                                   tCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTg                                                                                  >  4:539823/1‑146 (MQ=255)
                                                                      ggAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCAt                                                                             >  1:1013872/1‑148 (MQ=255)
                                                                       gAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCgg                                                                                                                                   >  3:149716/1‑93 (MQ=255)
                                                                       gAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCgg                                                                                                                                   <  4:149716/93‑1 (MQ=255)
                                                                       gAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAg                                                                                                                                >  1:353003/1‑96 (MQ=255)
                                                                       gAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAg                                                                                                                                <  2:353003/96‑1 (MQ=255)
                                                                        aaaCGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCa                                                                                                      >  3:869277/1‑121 (MQ=255)
                                                                        aaaCGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCa                                                                                                      <  4:869277/121‑1 (MQ=255)
                                                                             gACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTg                                                                                             >  3:963579/1‑125 (MQ=255)
                                                                             gACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTg                                                                                             <  4:963579/125‑1 (MQ=255)
                                                                             gACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCa                                                                       >  3:261344/1‑147 (MQ=255)
                                                                             gACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCa                                                                       <  4:261344/147‑1 (MQ=255)
                                                                              aCAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                                                                                               >  1:830427/1‑122 (MQ=255)
                                                                              aCAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCAtctc                                                                                               <  2:830427/122‑1 (MQ=255)
                                                                              aCAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGt                                                                     >  4:97487/1‑148 (MQ=255)
                                                                               cAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGc                                                                                                    >  3:780011/1‑116 (MQ=255)
                                                                               cAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGc                                                                                                    <  4:780011/116‑1 (MQ=255)
                                                                                   cTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTgg                                                                                 <  4:196934/131‑1 (MQ=255)
                                                                                   cTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTgg                                                                                 >  3:196934/1‑131 (MQ=255)
                                                                                     gCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCAtt                                                                            <  3:815750/134‑1 (MQ=255)
                                                                                     gCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCAtt                                                                            >  4:815750/1‑134 (MQ=255)
                                                                                       gggAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGtt                                                            <  2:1013872/148‑1 (MQ=255)
                                                                                         gAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGc                                                                                                                                              <  1:1126809/64‑1 (MQ=255)
                                                                                         gAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGc                                                                                                                                              >  2:1126809/1‑64 (MQ=255)
                                                                                         gAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCg                                                                                                                                    >  4:291865/1‑74 (MQ=255)
                                                                                         gAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCg                                                                                                                                    <  3:291865/74‑1 (MQ=255)
                                                                                         gAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGc                                                          >  2:42516/1‑148 (MQ=255)
                                                                                            gACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTccc                                                                         <  3:357616/130‑1 (MQ=255)
                                                                                            gACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTccc                                                                         >  4:357616/1‑130 (MQ=255)
                                                                                              cATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTg                                                                                  <  2:911671/119‑1 (MQ=255)
                                                                                              cATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTg                                                                                  >  1:911671/1‑119 (MQ=255)
                                                                                              cATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGt                                                             >  4:142997/1‑140 (MQ=255)
                                                                                              cATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGt                                                             <  3:142997/140‑1 (MQ=255)
                                                                                                tCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCa                                                                                                                                             >  3:885330/1‑58 (MQ=255)
                                                                                                tCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCa                                                                                                                                             <  4:885330/58‑1 (MQ=255)
                                                                                                 cAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGt                                                  <  1:42516/148‑1 (MQ=255)
                                                                                                   aTGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCa                                                               <  1:258514/133‑1 (MQ=255)
                                                                                                   aTGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCa                                                               >  2:258514/1‑133 (MQ=255)
                                                                                                   aTGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtca                                                >  4:49829/1‑148 (MQ=255)
                                                                                                   aTGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtca                                                >  3:324052/1‑148 (MQ=255)
                                                                                                    tGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATc                                                                >  1:849592/1‑131 (MQ=255)
                                                                                                    tGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATc                                                                <  2:849592/131‑1 (MQ=255)
                                                                                                    tGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtcat                                               >  1:933260/1‑148 (MQ=255)
                                                                                                     ggCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTg                                                                                                                       <  3:480983/75‑1 (MQ=255)
                                                                                                     ggCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTg                                                                                                                       >  4:480983/1‑75 (MQ=255)
                                                                                                     ggCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtcatc                                              <  3:49829/148‑1 (MQ=255)
                                                                                                         tttGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCa                                                                       <  1:104176/119‑1 (MQ=255)
                                                                                                         tttGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCa                                                                       >  2:104176/1‑119 (MQ=255)
                                                                                                         tttGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtcat                                               <  1:992546/143‑1 (MQ=255)
                                                                                                         tttGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGtcat                                               >  2:992546/1‑143 (MQ=255)
                                                                                                         tttGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAggg                                          <  3:97487/148‑1 (MQ=255)
                                                                                                           tGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAg                                                                                                              <  1:895569/78‑1 (MQ=255)
                                                                                                           tGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAg                                                                                                              >  2:895569/1‑78 (MQ=255)
                                                                                                               aCCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCa                                                                                                               >  1:955284/1‑73 (MQ=255)
                                                                                                               aCCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCa                                                                                                               <  2:955284/73‑1 (MQ=255)
                                                                                                               aCCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGt                                       >  3:475889/1‑145 (MQ=255)
                                                                                                               aCCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGt                                       <  4:475889/145‑1 (MQ=255)
                                                                                                                 cAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCg                                        >  1:993754/1‑142 (MQ=255)
                                                                                                                 cAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCg                                        <  2:993754/142‑1 (MQ=255)
                                                                                                                 cAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACg                                  >  4:1139805/1‑148 (MQ=255)
                                                                                                                    ttCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAAt                                                                 <  2:894688/114‑1 (MQ=255)
                                                                                                                    ttCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAAt                                                                 >  1:894688/1‑114 (MQ=255)
                                                                                                                    ttCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTcg                               <  4:324052/148‑1 (MQ=255)
                                                                                                                    ttCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTcg                               <  2:933260/148‑1 (MQ=255)
                                                                                                                          aCGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTcg                               >  2:56806/1‑142 (MQ=255)
                                                                                                                          aCGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTcg                               <  1:56806/142‑1 (MQ=255)
                                                                                                                            gTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTg                       <  3:1139805/148‑1 (MQ=255)
                                                                                                                                    ttCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGt                                                                     <  4:377877/94‑1 (MQ=255)
                                                                                                                                    ttCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGt                                                                     >  3:377877/1‑94 (MQ=255)
                                                                                                                                       ggCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTc                           <  4:594970/133‑1 (MQ=255)
                                                                                                                                       ggCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTc                           >  3:594970/1‑133 (MQ=255)
                                                                                                                                          tgcaAGCCAGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAggg                                          <  2:324983/114‑1 (MQ=255)
                                                                                                                                          tgcaAGCCAGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAggg                                          >  1:324983/2‑115 (MQ=255)
                                                                                                                                                 tGGCTTTCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTgcgc  >  2:914948/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTAGTAAATGGCGAGCAAACCAGCACTCTTCAATACTCTCCAGCGCCGCCACTTCTTCCAGCGTATTCGGGAAACGGACAAACTGCGGGAAGACATCAATGGCTTTTGCGACCAGTTCGCCACGTTCGGCTTTCGGCAGCAAGCTGGCTTGCAGGTAAAGCGGCAGCGTGCGCTGGCAGTCCAGTGCGTTCAGCATCTCTGCCGCCAGGGTGGGCATTCCCCAGTTAATCAGTTGCCCGGCGGTCATCAGGGCGTTAACGTCGCGTCCCTGACACAGACCGCGCCACTGCGC  >  NC_000913/2680494‑2680786

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