Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,151 T→C 5.6% intergenic (+18/‑120) 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,1510TC 5.6% 960.6 / 27.9 284intergenic (+18/‑120)ygdQ/ygdRUPF0053 family inner membrane protein/DUF903 family verified lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (137/131);  new base C (10/6);  total (147/137)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 284 total aligned reads displayed.

GGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTG  >  NC_000913/2971007‑2971298
                                                                                                                                                |                                                                                                                                                   
   cTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAg                                                                                                                                              >  4:47368/1‑149 (MQ=255)
    tttACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAgg                                                                                                                                             >  1:982284/1‑149 (MQ=255)
         ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                                                        >  2:442396/1‑149 (MQ=255)
         ccTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTa                                                                                                                                        >  4:122057/1‑149 (MQ=255)
          cTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                       >  1:1354239/1‑149 (MQ=255)
              ttCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcac                                                                                                                                   <  4:1510491/149‑1 (MQ=255)
               tCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTacggagt                                                                                                                                  <  4:720633/149‑6 (MQ=255)
                    aaaGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                                                             <  4:448221/149‑1 (MQ=255)
                     aaGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtttt                                                                                                                            <  3:177426/149‑1 (MQ=255)
                       gTTTCGACATCCACGTACCGAAAGGTTACATCTAGTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAGAATAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                                                                          <  2:1489231/149‑1 (MQ=255)
                       gTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTAc                                                                                                                          >  3:717939/1‑149 (MQ=255)
                             aCATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                                                                    >  3:1439700/1‑149 (MQ=255)
                               aTCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAg                                                                                                                  >  4:6149/1‑149 (MQ=255)
                                tCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGt                                                                                                                 >  3:1043811/1‑149 (MQ=255)
                                   aCGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTa                                                                                                              >  3:569721/1‑149 (MQ=255)
                                    cGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                             <  2:962202/149‑1 (MQ=255)
                                      tACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGa                                                                                                           <  4:487521/149‑1 (MQ=255)
                                         cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                        >  2:982054/1‑149 (MQ=255)
                                         cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                        >  1:939314/1‑149 (MQ=255)
                                         cGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                        >  4:1651427/1‑149 (MQ=255)
                                          gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                                                                >  2:368685/1‑124 (MQ=255)
                                          gAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTcc                                                                                                                                <  1:368685/124‑1 (MQ=255)
                                           aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                      <  2:1359636/149‑1 (MQ=255)
                                             aGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCt                                                                                                    <  3:1518410/149‑1 (MQ=255)
                                                     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                            <  1:1227332/149‑1 (MQ=255)
                                                     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                            <  2:766436/149‑1 (MQ=255)
                                                     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                            <  3:447750/149‑1 (MQ=255)
                                                         cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCATCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                                                        <  1:1384760/149‑1 (MQ=255)
                                                         cTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                                                        <  3:140035/149‑1 (MQ=255)
                                                         cTTCGCGATGTTCCTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGata                                                                                        <  4:1435354/149‑1 (MQ=255)
                                                          ttCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGatat                                                                                       >  2:706236/1‑149 (MQ=255)
                                                               gATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAAc                                                                                  >  3:110544/1‑149 (MQ=255)
                                                                 tGTTCTTCTCTATTGCGGTTGAAAGCATCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGt                                                                                >  4:1442090/1‑149 (MQ=255)
                                                                   ttcttcTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCg                                                                              >  3:488638/1‑149 (MQ=255)
                                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                             <  2:1180227/115‑1 (MQ=255)
                                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTaa                                                                                                             >  1:1180227/1‑115 (MQ=255)
                                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGa                                                                                          >  3:376554/1‑134 (MQ=255)
                                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                           >  3:466587/1‑149 (MQ=255)
                                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                           >  4:267794/1‑149 (MQ=255)
                                                                        ctctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTac                                                                                                                                       >  2:1384921/1‑87 (MQ=255)
                                                                         tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                              <  3:564903/127‑1 (MQ=255)
                                                                         tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTc                                                                                              >  4:564903/1‑127 (MQ=255)
                                                                         tctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtat                                                                        >  3:969106/1‑149 (MQ=255)
                                                                          ctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtata                                                                       >  2:1310222/1‑149 (MQ=255)
                                                                           tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGat                                                                                         >  3:1218445/1‑130 (MQ=255)
                                                                           tATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATaa                                                                      >  2:1678794/1‑149 (MQ=255)
                                                                            aTTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTcc                                                                                                     >  2:165995/1‑117 (MQ=255)
                                                                               gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAg                                                                  >  1:1303091/1‑149 (MQ=255)
                                                                                cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATGTATAACGTCGGTTTATAAACAGa                                                                 <  1:1055527/149‑1 (MQ=255)
                                                                                cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                 <  3:438197/149‑1 (MQ=255)
                                                                                cGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGa                                                                 <  2:483960/149‑1 (MQ=255)
                                                                                    tGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtat                                                             <  1:1705542/149‑1 (MQ=255)
                                                                                     gAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTtat                                                                        <  4:342056/137‑1 (MQ=255)
                                                                                      aaaGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtatta                                                           >  3:629566/1‑149 (MQ=255)
                                                                                         gCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCa                                                        <  1:1533543/149‑1 (MQ=255)
                                                                                          ccTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAt                                                       <  1:1500941/149‑1 (MQ=255)
                                                                                           cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGa                                                                                                           >  2:154337/1‑96 (MQ=255)
                                                                                           cTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta                                                      <  1:1660388/149‑1 (MQ=255)
                                                                                            tCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                            >  4:182674/1‑110 (MQ=255)
                                                                                                cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggtag                                                 >  2:1441420/1‑149 (MQ=255)
                                                                                                cTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAtaggtag                                                 >  3:1295804/1‑149 (MQ=255)
                                                                                                   gATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCg                                                                              >  1:1593691/1‑117 (MQ=255)
                                                                                                      tCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCaaa                                           <  4:502712/149‑1 (MQ=255)
                                                                                                        gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                            <  3:239485/114‑1 (MQ=255)
                                                                                                        gCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                            >  4:239485/1‑114 (MQ=255)
                                                                                                             aaaanGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGAt                                    >  3:1591800/1‑149 (MQ=255)
                                                                                                              aaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATa                                   >  1:800012/1‑149 (MQ=255)
                                                                                                                aaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAAc                                 <  2:243046/149‑1 (MQ=255)
                                                                                                                aaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTc                                                                                                      >  2:300952/1‑79 (MQ=255)
                                                                                                                 aGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCtt                                                                                                                              >  1:195167/1‑55 (MQ=255)
                                                                                                                 aGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAAcc                                <  4:1028369/149‑1 (MQ=255)
                                                                                                                   aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGGGGGCGTACTCACTCCTTTTACCCTGACAGTATCAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGAAGGCCAACTGAGGATACCCga                              >  4:846894/1‑148 (MQ=255)
                                                                                                                   aaTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCgg                              <  2:1215093/149‑1 (MQ=255)
                                                                                                                    aTCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGa                             <  1:1008846/149‑1 (MQ=255)
                                                                                                                     tCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGAt                            >  2:567013/1‑149 (MQ=255)
                                                                                                                        gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACtattat                                                          >  1:736287/1‑116 (MQ=255)
                                                                                                                        gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaa                         >  1:6284/1‑149 (MQ=255)
                                                                                                                        gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaa                         >  3:1266863/1‑149 (MQ=255)
                                                                                                                        gCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaa                         >  1:882266/1‑149 (MQ=255)
                                                                                                                          tctGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGa                                        <  2:205867/132‑1 (MQ=255)
                                                                                                                            tGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaaa                     >  3:299612/1‑149 (MQ=255)
                                                                                                                            tGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaaa                     >  3:118133/1‑149 (MQ=255)
                                                                                                                               tACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                            >  2:534358/1‑75 (MQ=255)
                                                                                                                                 cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGc                >  3:882690/1‑149 (MQ=255)
                                                                                                                                  tCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCa               >  2:1017678/1‑149 (MQ=255)
                                                                                                                                    agTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGt             >  4:1187272/2‑149 (MQ=255)
                                                                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                                                                   <  2:1156408/94‑1 (MQ=38)
                                                                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACa                                                                   >  1:1156408/1‑94 (MQ=38)
                                                                                                                                     gTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGta            >  1:288937/1‑149 (MQ=255)
                                                                                                                                      tACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaa           <  1:706236/149‑1 (MQ=255)
                                                                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGtt                                                                           <  2:1347590/85‑1 (MQ=38)
                                                                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACAGGATGAAAAAATGGGCAGtaa           >  3:750715/1‑149 (MQ=255)
                                                                                                                                        cGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGtaata         >  4:1346546/1‑149 (MQ=255)
                                                                                                                                         gccctcCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATg                   >  3:912263/4‑138 (MQ=255)
                                                                                                                                           tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGAt                                    <  4:892058/119‑1 (MQ=255)
                                                                                                                                           tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt      >  3:1007686/1‑149 (MQ=255)
                                                                                                                                           tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt      >  2:882518/1‑149 (MQ=255)
                                                                                                                                           tctcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAAtt      >  1:666850/1‑149 (MQ=255)
                                                                                                                                             tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc    <  4:432327/149‑1 (MQ=255)
                                                                                                                                             tcCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTc    >  4:1296998/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGGCTTTACCCTGATTCTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCTG  >  NC_000913/2971007‑2971298

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