Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,200,185 (A)7→6 intergenic (+91/‑96) yjaG → / → hupA DUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,200,1790A.99.2% 597.4 / ‑2.0 128intergenic (+85/‑102)yjaG/hupADUF416 domain protein/HU, DNA‑binding transcriptional regulator, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (54/73);  total (54/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 134 total aligned reads displayed.

AGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAG  >  NC_000913/4200031‑4200323
                                                                                                                                                    |                                                                                                                                                
aGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaa                                                                                                                                                <  4:986071/151‑1 (MQ=255)
    cTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaa                                                                                                                                              >  4:1085886/1‑149 (MQ=255)
    cTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaa                                                                                                                                              >  4:391249/1‑149 (MQ=255)
    cTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaa                                                                                                                                              >  1:313349/1‑149 (MQ=255)
      tAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaagt                                                                                                                                           <  4:726013/150‑3 (MQ=255)
       aGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGaaaaaagtg                                                                                                                                          <  3:384560/150‑4 (MQ=255)
            aGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATc                                                                                                                                   >  3:42609/1‑151 (MQ=255)
             gACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCg                                                                                                                                  <  4:582856/151‑1 (MQ=255)
             gACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCg                                                                                                                                  <  3:372430/151‑1 (MQ=255)
             gACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCCATCg                                                                                                                                  >  1:237947/1‑151 (MQ=255)
               ccTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGt                                                                                                                                <  1:385218/151‑1 (MQ=255)
                 tGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGc                                                                                                                              >  2:1054267/1‑151 (MQ=255)
                  gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTg                                                                                                                               <  1:307636/149‑1 (MQ=255)
                  gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGc                                                                                                                              <  4:1531185/150‑1 (MQ=255)
                  gCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCg                                                                                                                             <  4:410427/151‑1 (MQ=255)
                    gTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCgtg                                                                                                                           <  4:1363366/151‑1 (MQ=255)
                        ggCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATg                                                                                                                       <  4:355189/151‑1 (MQ=255)
                           gggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGc                                                                                                                      >  4:369297/1‑149 (MQ=255)
                           gggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAg                                                                                                                    <  3:5373/151‑1 (MQ=255)
                            ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                                                                   <  4:128438/151‑1 (MQ=255)
                            ggTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                                                                   <  4:928630/151‑1 (MQ=255)
                               gagaGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAgg                                                                                                                   >  4:107299/1‑148 (MQ=255)
                                 gagCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGagag                                                                                                               <  2:538798/150‑1 (MQ=255)
                                   gCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGc                                                                                                            <  4:1509196/151‑1 (MQ=255)
                                    cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTg                                                                                                             <  1:83222/149‑1 (MQ=255)
                                    cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCt                                                                                                           <  1:397024/151‑1 (MQ=255)
                                    cAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCt                                                                                                           >  4:680489/1‑151 (MQ=255)
                                       tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                                                                        >  2:966933/1‑151 (MQ=255)
                                       tatTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTAtt                                                                                                        >  3:460568/1‑151 (MQ=255)
                                          tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCt                                                                                                      <  2:719336/150‑1 (MQ=255)
                                          tGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTg                                                                                                     >  1:443665/1‑151 (MQ=255)
                                                 aTTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTGCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTc                                                                                              >  1:216676/1‑151 (MQ=255)
                                                   ttttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCAGTTcc                                                                                             <  1:162385/150‑1 (MQ=255)
                                                    tttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                                                                            <  2:1120500/150‑1 (MQ=255)
                                                     ttCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGGGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCg                                                                                            <  4:1247071/149‑1 (MQ=255)
                                                      tCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTc                                                                                          <  2:1387491/150‑1 (MQ=255)
                                                       cAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGc                                                                                        >  4:1137596/1‑151 (MQ=255)
                                                        aGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCa                                                                                       <  3:366035/151‑1 (MQ=255)
                                                          cAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACt                                                                                     >  2:1005296/1‑151 (MQ=255)
                                                            aTAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCg                                                                                   >  2:1238469/1‑151 (MQ=255)
                                                              aaGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGAt                                                                                 >  4:139090/1‑151 (MQ=255)
                                                               aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                <  3:313366/151‑1 (MQ=255)
                                                               aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                <  2:606567/151‑1 (MQ=255)
                                                               aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                >  2:230016/1‑151 (MQ=255)
                                                               aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                <  4:328450/151‑1 (MQ=255)
                                                               aGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGCGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATg                                                                                >  4:572278/1‑151 (MQ=255)
                                                                 aCCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCt                                                                              <  2:147415/151‑1 (MQ=255)
                                                                   cAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTa                                                                            <  4:1023813/151‑1 (MQ=255)
                                                                    aGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCtt                                                                             >  2:1320176/1‑149 (MQ=255)
                                                                      aaaaCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTACGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGc                                                                          >  1:150885/1‑150 (MQ=255)
                                                                         aCGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGc                                                                      <  2:1208164/151‑1 (MQ=255)
                                                                             gATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                                                                   <  2:970276/150‑1 (MQ=255)
                                                                             gATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGAt                                                                   >  2:120297/1‑150 (MQ=255)
                                                                              aTTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATaa                                                                 <  2:104064/151‑1 (MQ=255)
                                                                                 tAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAac                                                               >  1:327693/1‑150 (MQ=255)
                                                                                   aCGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                                                             >  3:312784/1‑150 (MQ=255)
                                                                                   aCGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAacac                                                             >  2:713823/1‑150 (MQ=255)
                                                                                    cGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                                                           >  1:461161/1‑151 (MQ=255)
                                                                                     gCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAt                                                           <  4:114739/151‑1 (MQ=255)
                                                                                     gCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACAtt                                                          <  2:894068/151‑1 (MQ=255)
                                                                                         gATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGt                                                        <  4:965237/149‑1 (MQ=255)
                                                                                         gATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTa                                                       <  2:421722/150‑1 (MQ=255)
                                                                                           tttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAgg                                                    >  2:1574375/1‑151 (MQ=255)
                                                                                            ttGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAgg                                                    <  4:645312/150‑1 (MQ=255)
                                                                                             tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                  >  4:1682143/1‑151 (MQ=255)
                                                                                             tGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                  >  4:1454622/1‑151 (MQ=255)
                                                                                              gTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGAt                                                  <  2:390689/150‑1 (MQ=255)
                                                                                                 gTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACt                                              <  4:1186462/151‑1 (MQ=255)
                                                                                                   aCCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACtt                                             <  2:1177352/150‑1 (MQ=255)
                                                                                                   aCCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTa                                            <  4:1670529/151‑1 (MQ=255)
                                                                                                     cTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATg                                          >  3:117503/1‑151 (MQ=255)
                                                                                                         aGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACa                                      <  1:318030/151‑1 (MQ=255)
                                                                                                          gTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACa                                      >  3:229944/1‑150 (MQ=255)
                                                                                                          gTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACa                                      >  2:588177/1‑150 (MQ=255)
                                                                                                            cTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACa                                      >  1:473940/1‑148 (MQ=255)
                                                                                                             ttCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGAc                                  >  2:508284/1‑151 (MQ=255)
                                                                                                             ttCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGAc                                  >  2:1242645/1‑151 (MQ=255)
                                                                                                              tCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGa                                   <  2:1524956/149‑1 (MQ=255)
                                                                                                                ccTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTc                                >  4:586671/1‑150 (MQ=255)
                                                                                                                    tCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACt                            <  4:148257/150‑1 (MQ=255)
                                                                                                                    tCGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACt                            <  2:1351192/150‑1 (MQ=255)
                                                                                                                     cGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtg                           >  4:1289801/1‑150 (MQ=255)
                                                                                                                     cGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtg                           <  2:679553/150‑1 (MQ=255)
                                                                                                                     cGCCCCCCGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtg                           >  2:1602863/1‑150 (MQ=255)
                                                                                                                        cccccGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtgattg                       >  2:1600418/1‑151 (MQ=255)
                                                                                                                        cccccGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtgattg                       <  3:326933/151‑1 (MQ=255)
                                                                                                                         ccccGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACtgattga                      <  2:138856/151‑1 (MQ=255)
                                                                                                                            cGTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGt                   <  4:809300/151‑1 (MQ=255)
                                                                                                                             gTCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGt                   >  3:185648/1‑150 (MQ=255)
                                                                                                                              tCTGGTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTa                  >  3:149330/1‑150 (MQ=255)
                                                                                                                                 ggTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAAtt               <  1:24674/150‑1 (MQ=255)
                                                                                                                                 ggTCTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAAt                <  2:553581/149‑1 (MQ=255)
                                                                                                                                    cTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGc             >  2:1624714/1‑149 (MQ=255)
                                                                                                                                    cTACATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGc             >  2:897385/1‑149 (MQ=255)
                                                                                                                                      aCATTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCa            <  1:122096/148‑1 (MQ=255)
                                                                                                                                        aTTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGaa       <  2:803680/151‑1 (MQ=255)
                                                                                                                                        aTTTGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGaa       <  4:812571/151‑1 (MQ=255)
                                                                                                                                         tttGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGaa       >  4:304044/1‑150 (MQ=255)
                                                                                                                                         tttGGGGGGCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGaa       >  2:1043464/1‑150 (MQ=255)
                                                                                                                                            ggggggCG‑AAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGc    <  4:136283/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGGCCTTAGGGCAGACCTGCGTGAGGCGGGTGAGAGCAATATTGGTATAATTTTTCAGCAATAAGACCAGAAAACGTGATTTAACGCCTGATTTGTCGTACCTGGAGTCTTCCCTTTCGCCCCCCGTCTGGTCTACATTTGGGGGGCGAAAAAAAGTGGCTATCGGTGCGTGTATGCAGGAGAGTGCTATTCTGGCATTTCCGTCGCACTCGATGCTTAGCAAGCGATAAACACATTGTAAGGATAACTTATGAACAAGACTCAACTGATTGATGTAATTGCAGAGAAAGCAG  >  NC_000913/4200031‑4200323

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