Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,508,818 T→C 5.0% intergenic (+76/+16) rlhA → / ← yncJ 23S rRNA 5‑hydroxycytidine C2501 synthase/protein YncJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,508,8180TC 5.0% 245.6 / 2.6 79intergenic (+76/+16)rlhA/yncJ23S rRNA 5‑hydroxycytidine C2501 synthase/protein YncJ
Reads supporting (aligned to +/- strand):  ref base T (44/31);  new base C (2/2);  total (46/33)
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.49e-01

TTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCCACA  >  NC_000913/1508671‑1508963
                                                                                                                                                   |                                                                                                                                                 
ttACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATcc                                                                                                                                                 <  1:745148/150‑1 (MQ=255)
  aCCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCgg                                                                                                                                               <  2:517090/150‑1 (MQ=255)
       gAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGTGTTACGCCGCATCCGGCATCa                                                                                                                                          <  1:911625/150‑1 (MQ=255)
        aaGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGCaggca                                                                                                                                         <  2:26696/150‑5 (MQ=255)
          gCGTTGTGCCATCCGTAAGTCCGGTTGAGCTGCTGAAAACATTGCCGAAGGGAAAAGGGTAAAACGCCAGTTTTCTGGCTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGATGCGGCGTAACGCCTTACCCGGCaggcaat                                                                                                                                       <  1:803102/150‑7 (MQ=37)
          gCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAAt                                                                                                                                       >  1:536278/1‑150 (MQ=255)
            gTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGa                                                                                                                                     <  2:1297008/150‑1 (MQ=255)
            gTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGa                                                                                                                                     >  1:792026/1‑150 (MQ=255)
             tAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGAc                                                                                                                                    <  1:937973/150‑1 (MQ=255)
                    aTCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGt                                                                                                                             <  1:1441631/150‑1 (MQ=255)
                            gTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGc                                                                                                                     >  1:1214104/1‑150 (MQ=255)
                             tCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCt                                                                                                                    >  2:137154/1‑150 (MQ=255)
                                     aGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGttt                                                                                                            >  2:1439979/1‑150 (MQ=255)
                                      gctgctGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTg                                                                                                           >  1:349384/1‑150 (MQ=255)
                                        tgctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTgcg                                                                                                         <  1:137154/150‑1 (MQ=255)
                                         tctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCAATTTTTCTGGTTACTCACAACTTATTGAATCTGAATGATATTTCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAAAGTTGTAACGCTTCTGGTTTgcgc                                                                                                        >  1:806121/2‑150 (MQ=255)
                                          ctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGcgcg                                                                                                       >  1:1504785/1‑150 (MQ=255)
                                          ctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGcgcg                                                                                                       >  1:944189/1‑150 (MQ=255)
                                          ctgAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGcgcg                                                                                                       >  2:1486422/1‑150 (MQ=255)
                                            gAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGAt                                                                                                     >  1:990232/1‑150 (MQ=255)
                                            gAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGAt                                                                                                     <  1:450514/150‑1 (MQ=255)
                                                 cATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTg                                                                                                <  1:583270/150‑1 (MQ=255)
                                                   ttGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGcc                                                                                              <  1:1065676/150‑1 (MQ=255)
                                                    tGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCa                                                                                             <  2:792026/150‑1 (MQ=255)
                                                     gCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAg                                                                                            >  1:92854/1‑150 (MQ=255)
                                                       cGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTa                                                                                          >  2:113289/1‑150 (MQ=255)
                                                               aaaGGGTAAAACGCCAGTATTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCt                                                                                  >  2:1409894/1‑150 (MQ=255)
                                                                 aGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCACATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTcg                                                                                >  2:198310/1‑150 (MQ=255)
                                                                    gTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCgcgc                                                                             >  2:1186072/1‑150 (MQ=255)
                                                                    gTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCgcgc                                                                             >  1:426143/1‑150 (MQ=255)
                                                                    gTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCgcgc                                                                             >  2:1350644/1‑150 (MQ=255)
                                                                    gTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCgcgc                                                                             >  2:67103/1‑150 (MQ=255)
                                                                      aaaaCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCaa                                                                           >  1:179065/1‑150 (MQ=255)
                                                                         aCGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAAcc                                                                        <  1:88216/150‑1 (MQ=255)
                                                                            ccAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTc                                                                     >  1:384249/1‑150 (MQ=255)
                                                                             cAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTcc                                                                    >  1:1119191/1‑150 (MQ=255)
                                                                                ttttCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTccgcc                                                                 <  2:384249/150‑1 (MQ=255)
                                                                                   tCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGAc                                                              <  2:944189/150‑1 (MQ=255)
                                                                                   tCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGAc                                                              <  1:387952/150‑1 (MQ=255)
                                                                                   tCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGAc                                                              <  1:1289538/150‑1 (MQ=255)
                                                                                   tCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGAc                                                              <  2:536278/150‑1 (MQ=255)
                                                                                             cacaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCca                                                    >  1:908369/1‑150 (MQ=255)
                                                                                             cacaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCca                                                    >  2:1333528/1‑150 (MQ=255)
                                                                                              acaACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCcac                                                   <  2:179065/150‑1 (MQ=255)
                                                                                                  cTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCa                                               >  2:607796/1‑150 (MQ=255)
                                                                                                  cTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCa                                               >  1:30794/1‑150 (MQ=255)
                                                                                                  cTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCa                                               >  2:602816/1‑150 (MQ=255)
                                                                                                        gAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCTCGCAAt                                         >  2:793402/1‑150 (MQ=255)
                                                                                                        gAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAAt                                         >  1:772797/1‑150 (MQ=255)
                                                                                                        gAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAAt                                         >  2:559232/1‑150 (MQ=255)
                                                                                                            cTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGACATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCAt                                     >  1:669227/1‑150 (MQ=255)
                                                                                                             tGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATc                                    >  2:214331/1‑150 (MQ=255)
                                                                                                              gCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCg                                   >  2:871592/1‑150 (MQ=255)
                                                                                                              gCAGGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGTTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACATCACGCAATTCATCg                                   <  2:1156743/150‑1 (MQ=255)
                                                                                                               cATGATATTGCGTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCAATTAAAAATGCTCGCGAAAACCTTACTCCTACTCTTCCTCCACATCACTAAATTAATcta                                  <  2:752003/150‑3 (MQ=255)
                                                                                                               cATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGc                                  <  1:1486422/150‑1 (MQ=255)
                                                                                                               cATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGc                                  <  1:602816/150‑1 (MQ=255)
                                                                                                               cATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGc                                  <  1:607796/150‑1 (MQ=255)
                                                                                                                 tGATATTGCCTGCCGGGTAAGGTGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTg                                >  1:1071325/1‑150 (MQ=255)
                                                                                                                 tGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTg                                >  1:1292465/1‑150 (MQ=255)
                                                                                                                    catttgatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGCCTCAACGTTGTCCCGCTTCTGGTTTGCGCAATTTTTGCCAGTAAAAATGCTCGCGCCAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCCTCGCTGTcc                             >  1:978022/8‑150 (MQ=37)
                                                                                                                     atTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCg                            >  2:1254573/1‑150 (MQ=255)
                                                                                                                      ttGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGc                           >  2:427205/1‑150 (MQ=255)
                                                                                                                         gatgccGGATGCGGCGTAACGCCTTACCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTTGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTc                        >  2:150769/3‑150 (MQ=255)
                                                                                                                          ctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCcac                                                   <  1:915928/122‑1 (MQ=255)
                                                                                                                          ctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCcac                                                   >  2:915928/1‑122 (MQ=255)
                                                                                                                          ctgccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCa                       <  1:1409894/150‑1 (MQ=255)
                                                                                                                            gccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATg                     <  1:1254573/150‑1 (MQ=255)
                                                                                                                            gccGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATg                     >  1:1279267/1‑150 (MQ=255)
                                                                                                                                ggTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGc                 >  1:1499776/1‑150 (MQ=255)
                                                                                                                                     ggCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACTCAATTCATCGCTGTCCGCTTCATGACGCATCtt            <  2:349384/150‑2 (MQ=255)
                                                                                                                                      gCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCa                                               <  2:131838/114‑1 (MQ=255)
                                                                                                                                      gCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCa                                               >  1:131838/1‑114 (MQ=255)
                                                                                                                                      gCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCtga           <  1:214331/150‑1 (MQ=255)
                                                                                                                                        gTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCtgatg         <  2:990232/150‑1 (MQ=255)
                                                                                                                                          tACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTgatgat       >  2:685863/1‑150 (MQ=255)
                                                                                                                                            cgccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCAAGCAATTCATCGCTGTCCGCTTCATGACGAAGCTGATGATcc     >  2:1173976/1‑150 (MQ=255)
                                                                                                                                              ccgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCcac   <  2:92854/150‑1 (MQ=255)
                                                                                                                                               cgcATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCcaca  >  1:791692/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTACCGGGAAGCGTAGTGGCATCCGTAAGTCCGGATGAGCTGCTGAAAACATTGCCGAAGCGAAAAGGGTAAAACGCCAGTTTTCTGGTTACTCACAACTTATTGAATCTGCATGATATTGCCTGCCGGGTAAGGCGTTACGCCGCATCCGGCATCAAATGACTCAACGTTGTCCCGCTTCTGGTTTGCGCGATTTTTGCCAGTAAAAATGCTCGCGCAAACCTTCCGCCGACTCTTCCGCCACAGCACGCAATTCATCGCTGTCCGCTTCATGACGCAGCTGATGATCCACA  >  NC_000913/1508671‑1508963

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