Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,159 G→A 5.8% intergenic (+26/‑112) ygdQ → / → ygdR UPF0053 family inner membrane protein/DUF903 family verified lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1590GA 5.8% 350.4 / 7.5 120intergenic (+26/‑112)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (47/66);  new base A (4/3);  total (51/69)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 124 total aligned reads displayed.

TTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACT  >  NC_000913/2971011‑2971308
                                                                                                                                                    |                                                                                                                                                     
tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAggg                                                                                                                                                      <  1:1104159/150‑1 (MQ=255)
tttACACTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGAGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCTGGAggg                                                                                                                                                      >  2:1521433/1‑150 (MQ=255)
  tACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGTTTACATCTACTTCGCGATGTTCTTCTCTTTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCg                                                                                                                                                    <  1:453149/150‑1 (MQ=255)
    cccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                                                                  >  1:1653945/1‑150 (MQ=255)
     ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                                 >  1:1042306/1‑150 (MQ=255)
     ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                                 <  1:1156159/150‑1 (MQ=255)
     ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTaa                                                                                                                                                 >  1:855718/1‑149 (MQ=255)
     ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                                                                                 <  1:1436634/150‑1 (MQ=255)
      cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTact                                                                                                                                                <  1:6249/150‑1 (MQ=255)
      cTGATTCTGGAAAGTTGCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGGTGAAAGCCGCAACTTGATTCGCAACAAAAAGAATCCGCGCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacg                                                                                                                                                <  2:687586/150‑2 (MQ=255)
        gATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacgga                                                                                                                                              <  2:936472/150‑4 (MQ=255)
           tCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcactc                                                                                                                                           >  1:431145/1‑150 (MQ=255)
            cTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                                                                          >  2:1820027/1‑150 (MQ=255)
            cTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCATACTCACTcc                                                                                                                                          >  1:1587273/1‑150 (MQ=255)
            cTGGAAAGTTTCGACATCCACGGACCGAAAGGTTACATCTACTTCGCGCTGTTCTTCTCTATTGCGGTTGACAGCCTCAACTTGATTCGCACCAAACAGAATCCGCGCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                                                                          <  1:777091/150‑1 (MQ=255)
             tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCt                                                                                                                                         >  2:1678596/1‑150 (MQ=255)
             tGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCATGCGGGAGGGCGTACTCTCTCCt                                                                                                                                         >  2:288855/1‑150 (MQ=255)
                 aaGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTa                                                                                                                                     <  2:67546/150‑1 (MQ=255)
                  aGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                                                                                    <  1:1345729/150‑1 (MQ=255)
                       cGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTca                                                                                                                               >  2:772608/1‑150 (MQ=255)
                        gACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                                                                              >  2:904755/1‑150 (MQ=255)
                          cATCCACGTACCCAAAGGTTACATCTACTTCGCGATTTTCTTCTCTATTGCGGATGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTACGTACGCTCTCCGGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                                                                            <  1:56062/150‑1 (MQ=255)
                             ccACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGttt                                                                                                                         >  2:1475226/1‑150 (MQ=255)
                               aCGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                                       <  1:1739247/150‑1 (MQ=255)
                                cGTACCGAAAGGTTTCATCTACTTCGCGATGTTATTCTCTATTGCGGTTGAAAGCCGCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCGACGGCGGGAGGGCGGACTCACTCCTTTTTCCCTCACAGTTTaat                                                                                                                      >  2:821029/1‑149 (MQ=255)
                                cGTACAGAAAGGTTACATCTACTTAGCGATGTTCTTCTCTATTGCGGTAGAAAGCCTAAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTACGTACTCTATCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                                                                      <  1:771999/150‑1 (MQ=255)
                                      gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAGTTTc                                                                                                                >  1:1085130/1‑150 (MQ=255)
                                      gAAAGGTTACATCTACTTCGCGATGATCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                                <  1:62451/150‑1 (MQ=255)
                                       aaaGGTTACATCTACTACGCGATGTTCTTCTCTAATGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATTCTCCGTAAGCTCTCCGGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTcc                                                                                                               <  1:1017377/150‑1 (MQ=255)
                                        aaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                                                                              <  2:1874833/150‑1 (MQ=255)
                                         aGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGACTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                                                                             <  1:101370/150‑1 (MQ=255)
                                         aGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                                                                             <  2:1684903/150‑1 (MQ=255)
                                          ggTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTTCTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGc                                                                                                            <  2:1384293/150‑1 (MQ=255)
                                             tACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCttt                                                                                                         <  2:812810/150‑1 (MQ=255)
                                              aCATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTAAGCCCTCCTTTTACCCTCACAGTTTAAGAGTTTCCTGCTTTc                                                                                                        >  1:156793/1‑150 (MQ=255)
                                               cATCTACATCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAACAAGAATCCGCACTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCa                                                                                                       <  2:159642/150‑1 (MQ=255)
                                                 tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                                                                     <  1:210918/150‑1 (MQ=255)
                                                 tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAg                                                                                                     <  2:1947092/150‑1 (MQ=255)
                                                  cTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                                                                    >  1:1572104/1‑150 (MQ=255)
                                                  cTACTTCGCCTTGTTCTTCTCTATTGCGGTTGAAAGCCTCTACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTGCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                                                                    <  2:331622/150‑1 (MQ=255)
                                                   tACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                                                                   <  1:437646/150‑1 (MQ=255)
                                                    aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                                                                  <  2:809305/150‑1 (MQ=255)
                                                    aCTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                                                                  <  2:1505579/150‑1 (MQ=255)
                                                       tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatatat                                                                                               <  2:2000716/150‑1 (MQ=255)
                                                        cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                                                                              <  1:680852/150‑1 (MQ=255)
                                                        cgcgATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAtatata                                                                                              >  2:319115/1‑150 (MQ=255)
                                                           gATGTTCTTCTCTATTGCGGTTGAAAGCCCCAACTTGATTCGCAACAAAAAGACTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTACTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                                                                           <  2:381295/150‑1 (MQ=255)
                                                               ttcttcTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCgg                                                                                       <  1:1281904/150‑1 (MQ=255)
                                                                tcttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                                      <  2:238701/150‑1 (MQ=255)
                                                                 cttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                                     <  1:1362630/150‑1 (MQ=255)
                                                                  ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                                                                                    >  1:1890095/1‑150 (MQ=255)
                                                                  ttctCTATTGCGGTTGAAAGCCACAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGttt                                                                                    >  1:1002629/1‑150 (MQ=255)
                                                                    ctctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtat                                                                                  >  2:1683734/1‑150 (MQ=255)
                                                                         ttGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                                                                             <  1:772371/150‑1 (MQ=255)
                                                                           gCGGTTGAAAGCCTCAACTTGATTCGCAACAACAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGGGCGGGGGCCGGTTACCCTCACAGTTTAAGAGTTTGCTGCTTTCATGATTTATTACGTCGGGGTTTAAACAgc                                                                           >  2:141875/1‑149 (MQ=255)
                                                                           gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                           >  1:436788/1‑150 (MQ=255)
                                                                           gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACGCCTTTTACCCTCACAGTTTAAGATTTTCCTGCGTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                           >  1:1169269/1‑150 (MQ=255)
                                                                               ttGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                                                                       <  2:646111/150‑1 (MQ=255)
                                                                                  aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat                                                                    >  2:1203280/1‑150 (MQ=255)
                                                                                   aaGCCTCAACTTGATTCTCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATc                                                                   <  2:834132/150‑1 (MQ=255)
                                                                                   aaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTAGGTTTATAAACAGACTATTATc                                                                   >  1:1829667/1‑150 (MQ=255)
                                                                                    aGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAAATATATAACGTCGGTTTATAAACAGACTATTATCa                                                                  <  1:1678596/150‑1 (MQ=255)
                                                                                       cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtag                                                               <  2:1004081/150‑1 (MQ=255)
                                                                                            cTTGATTCGCAACAAAAAGAATCCTCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggtagg                                                          <  1:1767896/150‑1 (MQ=255)
                                                                                            cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCGCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggtagg                                                          <  2:1446909/150‑1 (MQ=255)
                                                                                             ttGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCCCTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGc                                                         >  2:1265409/1‑150 (MQ=255)
                                                                                              tGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGcc                                                        >  1:1019915/1‑150 (MQ=255)
                                                                                               gATTCGAAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCa                                                       <  1:1645386/150‑1 (MQ=255)
                                                                                                 ttCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCaaa                                                     <  2:1052874/150‑1 (MQ=255)
                                                                                                     caacatAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACATTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGac                                                 <  1:1421565/150‑2 (MQ=255)
                                                                                                     caacaaATAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTATTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAg                                                 <  1:1652022/150‑1 (MQ=255)
                                                                                                      aacaaAAAGAATCCTCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTACTGCTTTAAAGATATATAACGTCGGTTTATAAACGGACTATTATCATAGGTAGGCCAAATGAgg                                                <  1:944365/150‑1 (MQ=255)
                                                                                                      aacaaAAAGAATCCGCGCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAgg                                                <  2:642506/150‑1 (MQ=255)
                                                                                                        caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGAt                                              >  2:456591/1‑150 (MQ=255)
                                                                                                        caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGAt                                              >  2:1698102/1‑150 (MQ=255)
                                                                                                        caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGAt                                              >  1:1716633/1‑150 (MQ=255)
                                                                                                         aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATa                                             <  2:868578/150‑1 (MQ=255)
                                                                                                         aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATa                                             <  2:876216/150‑1 (MQ=255)
                                                                                                           aaaGAATCCGCTCTGATACTCCGTATGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAAc                                           >  2:1248950/1‑150 (MQ=255)
                                                                                                            aaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAAcc                                          <  2:1947526/150‑1 (MQ=255)
                                                                                                              gAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACATTTTAAGATTTTCCTGCTTTAAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCgg                                        <  1:772695/150‑1 (MQ=255)
                                                                                                                aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAAATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGAt                                      <  2:1106596/150‑1 (MQ=255)
                                                                                                                aTCCGATCAGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTATGATTATCCTGCTTTCAAGATATATAACGTCGGTTTATAACCAGACTATTATCATATGTAGGCCAAATGAGGATAACCGGAt                                      <  2:702904/150‑1 (MQ=255)
                                                                                                                     ctctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTCTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaa                                 <  1:904755/150‑1 (MQ=255)
                                                                                                                       ctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaaa                               >  2:931004/1‑150 (MQ=255)
                                                                                                                         gATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCTCAGTTTAAGATTTTCCTGCTTTCAATATATATAACGTCGGTTTATAAACAGACTATTATCATCGGTAGGCCAAATGAGGATAACCGGATGAAAAAATg                             <  1:969683/150‑1 (MQ=255)
                                                                                                                            aCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGc                          <  2:1007065/150‑1 (MQ=255)
                                                                                                                             cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCa                         >  2:1064434/1‑150 (MQ=255)
                                                                                                                              tgagTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAg                        >  2:905772/4‑150 (MQ=255)
                                                                                                                              tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAg                        >  1:221/1‑150 (MQ=255)
                                                                                                                                 gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTGCTGGTTTCAAGATATATAACGTCGGGTTATAAACATACTATTATCATAGGTAGGCGAATTGAGGTTAACCGGATGAAAAAATGGGCAGtaa                     >  1:1236516/1‑150 (MQ=255)
                                                                                                                                   aCGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTATCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGACAAATGAGGATAACCGGATGAAAAAATGGGCAGtaata                   >  1:1260614/1‑150 (MQ=255)
                                                                                                                                    cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAATTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaataa                  <  2:1550465/150‑1 (MQ=255)
                                                                                                                                    cGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTTTAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaataa                  >  2:510475/1‑150 (MQ=255)
                                                                                                                                      ctctCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACTGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt                >  2:508697/1‑150 (MQ=255)
                                                                                                                                           cTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGc           <  1:288855/150‑1 (MQ=255)
                                                                                                                                            tGCGGGAGGGCGTACTGACTCCTTTTACCCTCAGAGTTTAAGATATTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCATATGAGGATTACCGGATGAAAAAATGGGCAGTAATAATTTCTGCa          >  1:215065/1‑150 (MQ=255)
                                                                                                                                            tGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAACTGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCa          <  1:1458289/150‑1 (MQ=255)
                                                                                                                                              cGGGATGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATAAAACGTCGGTTTATAAACATACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGt        >  2:1052671/1‑150 (MQ=255)
                                                                                                                                                    ggCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAAGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACt  <  1:1475226/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTGCAGTCGGACT  >  NC_000913/2971011‑2971308

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